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УДК 621.1 
И.В. Макаров, А.Н. Сизов, Ф.Л. Чубаров 


ВЛИЯНИЕ СОДЕРЖАНИЯ ВОЗДУХА 
В РАБОЧЕЙ ЖИДКОСТИ 

СИСТЕМЫ РЕГУЛИРОВАНИЯ ТУРБОПРИВОДА 
НА ЕГО ХАРАКТЕРИСТИКИ 


В работе исследуется влияние объема растворенного воздуха в 
рабочем масле системы регулирования паровой турбины на динами- 
ческую устойчивость ее работы. Исследование ведется на линеари- 
зованной математической модели турбины, учитывающей ее дина- 
мику, применительно к двум конкретным типам турбин ОАО «КТЗ», 
широко применявшихся в энергетике с 70-х гг. Для них определено 
расчетное максимально-допустимое объемное содержание воздуха в 
рабочей жидкости системы регулирования с т.з. потери ее устой- 
чивости, a также общие формулы для аналогичной оценки на других 
типах турбин схожей конструкции. 


Ключевые слова: паровая турбина, регулирование, рабочее масло, 
растворенный воздух, устойчивость, динамика, модель, исследова- 
ние, анализ. 


Введение. Паровые турбины малой мощности (до 25 МВт) произ- 
водства ОАО «КТЗ» давно известны на рынке энергомашиностроения 
как в РФ, так и за рубежом (работают в 42-х странах мира) [1]. Их отли- 
чает высокая надежность и неприхотливость к условиям эксплуатации, 
а также простота монтажа и настройки, ставшие результатом многолет- 
ней отработки конструкций. В процессе такой отработки неоднократно 
возникали и прорабатывались вопросы обеспечения более качественной 
и устойчивой работы систем регулирования и маслоснабжения турбин, 
достаточно полно освещенные в профильной литературе [2-4]. Одним 
из таких вопросов по-прежнему является влияние характеристик рабо- 
чих жидкостей в системах регулирования и маслоснабжения (у паровых 
турбин КТЗ они, как правило, совмещенные) на стабильность поддер- 
жания заданной частоты вращения в условиях действия разного рода 
возмущений и помех. Общая тенденция, выявленная на нескольких по- 
колениях турбин, заключается в том, что по мере эксплуатации характе- 
ристики рабочего масла неизбежно ухудшаются (меняется хим. состав, 
происходит насыщение воздухом и т. п.). В силу того, что системы ре- 
гулирования турбин выпуска до 90-х гг. конструктивно были, в основ- 
ном, гидродинамического типа [5], предполагающего обширное исполь- 
зование гидравлических узлов и элементов, их зависимостью от каче- 
ства масла пренебрегать недопустимо. Отмечая достаточно широкую 
проработку данных вопросов в литературе [8], остановимся более по- 
дробно на влиянии растворенного воздуха в рабочем масле системы ре- 
гулирования (Тп-22 по ГОСТ 9972-74) на устойчивость работы паровой 
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турбины применительно к двум типам турбин ОАО «КТЗ» (ОК-2 — для 
турбогенератора и ОР-18П — для приводной турбины питательного 
насоса [5]), широко применявшихся в энергетике с 70-х гг. 

Для исследования такого влияния воспользуемся математической 
моделью паровой турбины КТЗ с гидродинамической системой регули- 
рования [6, 9]. 

Полагая процесс сжатия (расширения) изотермическим и прене- 
брегая изменением степени растворения воздуха в рабочей жидкости в 
этом время, запишем уравнения движения системы в виде: 


p 29. Е 
dt 
dn ф 
Т — +1 = —; 
" dt " $ 


i (1) 


где Ф, N, с, и — относительное изменение угловой скорости вращения po- 
тора турбины, хода золотника трансформатора давления, хода отсечного 
золотника и хода сервомотора регулирующих клапанов; Г, Г, Tos Т, — 
динамические постоянные; д o - степень неравномерности. 


Для системы линейных дифференциальных уравнений (2) крите- 
рий устойчивости Гурвица будет иметь вид [7]: 


А T, 2 
Ta ха TT. E Ts +T >0. (2) 
Ф а 
Учитывая, что 
T, —T, (1+ Као); 


ж 
T, —T, (1 Ко), 
ЕЯ ЕЯ 
где Т - „Г, — динамические постоянные трансформатора давления и от- 


сечного золотника при отсутствии воздуха в рабочей жидкости; U- от- 
носительное объемное содержание нерастворимого воздуха в рабочей 
C, V, 


т н.т. 


2 
Е "о (Ро =E ^o) 
нерастворимого воздуха B импульсной линии OT насоса регулятора до 
СУ 


О П.и. 
2 
Е. (Р, = Р. ^o) 
ент, учитывающий влияние нерастворимого воздуха в проточной HM- 
пульсной линии; Уу. и = объемы импульсной линии от насоса-регуля- 


н.т.? 


жидкости; K, = коэффициент, учитывающий влияние 


золотника трансформатора давления; К. = — коэффици- 
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тора до золотника трансформатора давления и проточной импульсной ли- 
нии; С, Со - жесткость пружин трансформатора давления и отсечного 


золотника; F,,F, — площадь золотников трансформатора давления и от- 
сечного золотника; ( Po- Ву) - номинальный напор насоса регулятора; 
(А x: Ву) — номинальный перепад давления рабочей жидкости на отсеч- 


ном золотнике, и учитывая также, что для системы регулирования паровых 
турбин КТЗ всегда справедливы неравенства [5]: 


тека ко 
отво 


Приравняв левую часть неравенства (2) к нулю, найдем границу 
устойчивости системы (1): 


(3) 


где, 


Обычно в системах регулирования паровых турбин КТЗ всегда 

T, «T, [5]. Тогда, пренебрегая в выражении (2) членом T. как величи- 
ной второго порядка малости, вместо (2) получим 

8, T, (T, +1, )> 0. (4) 

Откуда найдем границу устойчивости по параметру о в более 

простом виде 
* ж 
E m) 


= + (5) 
гр * + 
T,-K,* T; Ке 


В rex случаях, когда D Ka < 0.1-T;. :K,, и учитывая, что для Na- 
ровых турбин КТЗ всегда имеем T. « T. „вместо (5) получим 
z б Тое 


с (6) 
i. ` Ks 


Tp 


Выражение (6) получим также, если в системе (1) положить T, =0 


и найти границу устойчивости системы регулирования по параметру U. 
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Подставив из [10] количественные характеристики рабочих пара- 
метров турбин ОК-18П и ОК-2, входящих в (3), получаем, что расчетное 
максимально-допустимое объемное содержание воздуха в рабочей жид- 
кости V, (20) для турбин типа ОР-18П составляет около 3 %, для турбин 


типа ОК-2 — около 9 %. 


В таблицах и на рисунках 1-4 приведен результаты расчета влия- 
ния объема основных гидравлических импульсных линий в блоках ре- 
гулирования турбин ОК-18П и ОК-2 на их динамическую устойчивость 


по частоте вращения. 


Таблица 1 
Турбина ОР-18П 
У. „см? 150000 | 120000 | 100000 | 80000 | 40000 | 25000 
Vap% 3,07 3,18 3,25 3,33 3,50 3,57 
Анализ устойчивости системы регулирования машины ОР-18П 
3.5 | ^ < | 
EN 
34| Si | 
£33 1 
32| ^ ч | 


2 4 6 


турбины ОР-18П на ее устойчивость 


8 


10 


VH.T.,CM3 
Рис. 1. Анализ влияния геометрического объема импульсной линии 
от насоса-регулятора до золотника трансформатора давления 


12 


Таблица 2 
Турбина ОК-2 
V, „см“ 50000 | 40000 | 30000 | 20000 | 15000 | 10000 
Ур» 2 9,64 9,9 10,17 | 10,46 | 10,61 | 10,76 
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Анализ устойчивости системы регулирования машины ОК-2 
0.8 


1 1.5 2 2.5 3 3.5 4 4.5 5 
VH.T.,CM3 < 104 
Рис. 2. Анализ влияния геометрического объема импульсной линии 
от насоса-регулятора до золотника трансформатора давления 
турбины ОК-2 на ее устойчивость 


Таблица 3 
Турбина ОР-18П 
V „см“ 8200 7400 6800 5100 4900 2200 
Vp% 3,07 3,34 3,58 4,49 4,63 7,88 


Анализ устойчивости системы регулирования машины ОР-18П 
8 


2000 _ 3000 4000 5000 600 700 _ 8000 9000 
Уп.и.,см3 
Рис. 3. Анализ влияния объема проточной импульсной линии блока 
регулирования турбинът ОР-18П на ее устойчивость 


1 1 1 1 L 
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Таблица 4 
Турбина ОК-2 
V, „см? 1600 | 1500 | 1400 | 1300 | 1200 | 1100 
Vp% 9,58 10,13 10,75 11,44 12,23 | 13,14 
is Анализ устойчивости системы регулирования машины ОК-2 
13 PN 
- 
12.5 а 
12 < 
^ 
E 11.5 EU 
> ES 
11 ~ = 
10.5 ia СИ 
ж 


9 5 de — — —- 1 -1. rann 21. —L 
1100 1150 1200 1250 1300 1350 1400 1450 1500 1550 1600 
Vn.M.,cM3 
Рис. 4. Анализ влияния объема проточной импульсной линии блока 
регулирования турбины ОК-2 на ее устойчивость 


Выводы 

1. Получены расчетные формулы определения максимально-допу- 
стимого объемного содержания воздуха в рабочей жидкости системы 
регулирования типовой паровой турбины КТЗ малой мощности, исходя 
из условия обеспечения ее устойчивой работы. 

2. Расчетное максимально-допустимое объемное содержание воз- 
духа в рабочей жидкости для штатной конструкции блока регулирова- 
ния турбин типа ОР-18П составило около 3%, для турбин типа ОК-2 — 
около 8% с тенденцией к увеличению по мере увеличения объемов ос- 
новных импульсных линий в блоках регулирования. 

3. Направлением дальнейшего исследования целесообразно вы- 
брать аналогичный анализ для систем регулирования нового типа — 
электрогидравлических, которыми в настоящее время комплектуются 
более 90% паровых турбин ОАО «КТЗ» по требованиям заказчиков. 


СПИСОК ЛИТЕРАТУРЫ 


[1] Референции / Информационный сайт ОАО «КТЗ. - [Электронный pe- 
сурс]. – URL: http://www.oaoktz.ru/customers/references/ (дата обращения 
22.01.2020) 

[2] Капелович Б.Э. Эксплуатация паротурбинных установок. - М.: Энергия, 
1975. — 288 c. 


Рубрика «Машиностроение» (ISSN 2413-6220) 12 


[3] Бойко Е.А. Тепловые электрические станции. Паротурбинные знергети- 
ческие установки ТЭС: справочное пособие / Е.А. Бойко, К.В. Баженов, 
П.А. Грачев. — Красноярск: ИПЦ КГТУ, 2006. - 152 с. 

[4] Веллер В.Н. Регулирование и защита паровых турбин. - М.: Знергоато- 
миздат, 1985. — 104 c. 

[5] Паровые турбины малой мощности КТЗ / И. Кирюхин, Н.М. Тараненко, 
Е.П. Огурцова и др. — М.: Энергоатомиздат, 1987. — 216 c. 

[6] Автоматическое регулирование энергоустановок: учеб. пособие для ву- 
зов / А.Е. Булкин. – М.: Издательский дом МЭИ, 2009. – 508 с.: ил. 

[7] Автоматическое регулирование паровых турбин и газотурбинных устано- 
вок: учеб. для вузов по спец. "Турбиностроение" / И. И. Кириллов. —2-е изд., 
перераб. и доп. — Л.: Машиностроение. Ленингр. отд-ние, 1988. — 445 с. 

[8] Гидравлика, гидромашины и гидроприводы: учебник для машинострои- 
тельных вузов / Т.М. Башта, С.С. Руднев, Б.Б. Некрасов и др. - 2-е изд., 
перераб. — М.: Машиностроение, 1982. — 423 с., ил. 

[9] Матричные методы расчета и проектирования сложных систем автома- 
тического управления для инженеров / Под ред. К.А. Пупкова и Н.Д. 
Егупова. — М.: Издательство МГТУ им. Н.Э. Баумана, 2007. — 664 с. 

[10] Информационные материалы ОАО «КТЗ» (технический отчет № 100-M-2570). 


СВЕДЕНИЯ ОБ АВТОРАХ 


Макаров Игорь Владиславович (МаКагоу Т.У.) - студент 
КФ МГТУ им. Н.Э. Баумана; таКат5883 @ gmail.com 


Сизов Александр Николаевич (5170у А.М.) — канд. техн. наук, до- 
цент кафедры «Тепловые двигатели и гидромашины» МГТУ 
им. Н.Э. Баумана; axelxd 9 yandex.ru 


Чубаров Федор Леонидович (Chubarov F.L.) — канд. техн. наук, зав. 
кафедрой «Механизации сельскохозяйственного производства» КФ ФГБОУ 
ВО "РГАУ – МСХА имени К.А. Тимирязева"; једог 1966сћир Ф yandex.ru 


AN INFLUENCE OF DISSOLVED AIR CONTAIN 
IN WORKING OIL OF A GOVERNING SYSTEM TURBINE 
DRIVE ON CHARACTERISTIC 


The paper investigates the influence of the volume of dissolved air in the 
working oil of a steam turbine control system on the dynamic stability of its 
operation. The study is conducted on a linearized mathematical model of a tur- 
bine, taking into account its dynamics, as applied to two specific types of tur- 
bines of KTZ OJSC, which have been widely used in power engineering since 
the 70s. For them, the calculated maximum permissible volumetric air content 
in the working fluid of the control system is determined with so on. loss of its 
stability, as well as general formulas for a similar assessment on other types of 
turbines of a similar design. 


Keywords: steam turbine, governing system, working oil, dissolved air, sta- 
bility, dynamics, model, research, analysis. 
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УДК 621 
В.М. Попков, Д.А. Тарасенков 


ИССЛЕДОВАНИЕ ВЛИЯНИЯ УСИЛИЯ НАТЯЖЕНИЯ 
ОБМОТОЧНОГО ПРОВОДА НА ЭЛЕКТРОМАГНИТНЫЕ 
ХАРАКТЕРИСТИКИ ПРИ ИЗГОТОВЛЕНИИ ОБМОТОК 
ЭЛЕКТРИЧЕСКИХ МАШИН 


В работе проводится исследование влияния усилия натяжения об- 
моточного провода на электромагнитные характеристики при изго- 
товлении обмоток электрических машин. Определена зависимость 
сопротивления обмотки и его индуктивности от величины усилия 
натяжения провода. Установлено влияние погрешности диаметра 
провода и погрешности удельного сопротивления на сопротивление 
обмоточного провода. 


Ключевые слова: электрическая машина, обмотка электрической 
машины, обмоточный провод, сопротивление провода. 


Введение. Обмотки являются важнейшими функциональными 
элементами электрических машин. От качества изготовления обмоток в 
немалой степени зависит качество изделий [1]. Дефекты обмоток отно- 
сятся к категории скрытых, а последствия их появления могут привести 
к внезапным и неконтролируемым отказам машин, поэтому проблема 
повышения качества обмоток имеет первоочередное значение, опреде- 
ляющее в конечном итоге над ежность функционирования изделия, а 
также производительность и себестоимость его изготовления [2]. 

Направления обеспечения бездефектности, повышения эффектив- 
ности процесса изготовления изделий известны: тщательная отработка 
технологии; детальное исследование причин и условий возникновения 
дефектов с целью их предотвращения; автоматизированное выполнение 
всех составляющих процесса «конструирование — производство» изде- 
лий, содержащих обмотки; совершенствование контрольно-испытатель- 
ных операций и их оснащения [3]. 

При разработке технологических процессов производства обмоток 
электрических машин учитывается влияние разнообразных факторов на 
точность получения их выходных характеристик. Для обоснованного и 
рационального назначения разброса электрических характеристик и их 
увязки с технологией изготовления обмоток целесообразно изучить воз- 
можные отклонения исходных параметров материалов и деталей, а также 
погрешности, вносимые в процессе намотки, и причины их появления. 

В процессе намотки провод испытывает различные воздействия, и 
его начальные характеристики претерпевают изменения. Таким обра- 
зом, в процессе исследования необходимо: установить фактические зна- 
чения исходных характеристик провода; определить величину измене- 
ния сопротивления провода при намотке. Первая задача сводится к ис- 
следованию производственных погрешностей обмоточного провода, а 
именно погрешности допуска на диаметр и погрешности удельного со- 
противления. 
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Производственные погрешности провода. При изготовлении 
обмоток электрических машин необходимо обеспечить получение со- 
противления обмотки в границах допустимых значений. На точность 
получения заданного параметра основное влияние оказывают режимы 
процесса намотки — натяжение провода, обеспечиваемое челноком 
намоточного станка, и скорость наматывания. Для исследования влия- 
ния указанных факторов (исходных отклонений и технологических по- 
грешностей), воспользуемся выражением для нахождения сопротивле- 
ния провода круглого сечения: 


41 
—;. (1) 


TD 
где p — удельное сопротивление провода, Ом iu M^; [— длина отрезка 
провода, M; D — диаметр провода, MM. 

Величина сопротивления изготовленной обмотки будет иметь рассеи- 
вание (погрешность параметра), обусловленное нестабильностью величин 
удельного сопротивления р, длины отрезка провода l и диаметра провода D. 

Чтобы установить зависимость погрешности параметра К от про- 
изводственных погрешностей величин выражения, воспользуемся пра- 
вилами дифференциального исчисления. Полный дифференциал урав- 
нения (1) имеет вид: 

ак = др+ 28: di zc лр, Q) 
пр TD лр 

Для относительного изменения сопротивления провода справед- 
ливо выражение: 


AR Ар Al 2AD 
гары (3) 
R p l D 
Рассмотрим влияние каждой составляющей суммарной погрешно- 


сти выражения (3). 
Погрешность сопротивления провода от нестабильности удель- 


А 
ного сопротивления АВ представим как да „а OT нестабильности диа- 


р 
метра провода как об, . Экспериментально установлено, что действи- 


тельное значение удельного сопротивления намоточного провода, как 
правило, отличается от стандартного (для медного провода ра = 


0,01754 Ом-мм/ м’). Тогда погрешность сопротивления обмотки, вы- 
званная отклонением действительной величины удельного сопротивле- 
ния провода ро от его стандартного значения Por: 


5R, = Pe Ро-100 %. (4) 
Per 
Погрешность сопротивления от колебания допуска на диаметр 
провода: 


2АР 
Rp = ———. 5 
Bos (5) 
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Следовательно, сопротивление провода изменяется не только в за- 
висимости от допуска на диаметр, но и от величины диаметра р, . 


Иными словами, чем меньше исходный диаметр провода, тем большее 
значение может принимать обусловленная им величина ÖR D: Допуски 


на диаметр провода стандартизированы. 

Зависимость погрешности сопротивления провода от его натя- 
жения при намотке. Наматывание провода на каркас производится с 
некоторым натяжением. Оптимальная величина натяжения для каждого 
вида намоточных проводов зависит от вида обмотки и предельно допу- 
стимой величины натяжения для данного провода. Чрезмерное натяже- 
ние влечет за собой недопустимое значение удлинения провода, изме- 
нение его сопротивления и повреждение изоляции. Недостаточное натя- 
жение вызовет погрешность формы, уменьшение коэффициента запол- 
нения обмоток и, как следствие, повлияет на качественные характери- 
стики намоточного изделия. 

Типовыми дефектами, вызванными выбором неоптимального уси- 
лия натяжения провода, являются хаотичное расположение провода (ри- 
сунок 1), «сбег» витка (рисунок 2) в случае недостаточного усилия и по- 
вреждение эмалированного покрытия провода (рисунок 3) в случае 
чрезмерного усилия [4]. 


4 


мА 
(LU MIT ШШ ЛА. 


Рис. 2. «Сбег» витка 
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Рис. 3. Повреждение змалированного покръпия провода 


На рисунке 4 представлена обмотка статора злектрической ма- 
шины, выполненная при оптимальном усилии натяжения обмоточного 


провода. 
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Рис. 4. Пример качественно уложенной обмотки статора 


В качестве справочных данных используются значения предельно до- 
пустимых натяжений конкретной марки и диаметра провода (таблица 1) [5]. 


Таблица 1 


Предельно допустимые натяжения медного провода ПЭТД-180 
ТУ16-705.264-82 при намотке 


Номинальный Номинальный 
Натяжение, г Натяжение, г 
диаметр, мм диаметр, мм 
0,03 7,5 0,12 120 
0,05 21 0,13 140 
0,07 42 0,14 165 
0,09 69 0,15 190 
0,10 85 0,16 220 
0,11 100 0,17 245 
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Изменение сопротивления провода от натяжения установим по 


формуле: 
Ю-В 
GR, = — — .100 96, (6) 


Ro 
где Rọ – первоначальное сопротивление провода при отсутствии раста- 
гивающих усилий, Ом; Ry - сопротивление провода, подвергнутого 


нагрузке в виде определ енного натяжения, Ом. 

Как было установлено, при приложении растягивающего усилия 
происходит уменьшение диаметра провода, что также влияет на элек- 
тромагнитные характеристики изделия. Это можно отразить следующей 
формулой нахождения индуктивности: 


= м 3 $4). (7) 


H 
где и — магнитная проницаемость воздуха, — ; | — начальная длина OT- 
м 


резка провода, MM; D — диаметр провода, мм; ly — длина провода, под- 
вергнутого нагрузке в виде определ енного натяжения, мм. 


Тогда относительное изменение индуктивности при изменении 
диаметра равно: 
АГ, _ р | 
L 4D-l-ly-4AD-l-ly -D 
Погрешности, определяющие точность сопротивления и индук- 
ТИВНОСТИ обмотки, подразделяются на систематические и случайные. 
Систематические погрешности необходимо учитывать при определении 
длины провода и числа витков обмотки. Случайные погрешности - 
следствие колебаний значений удельного сопротивления и диаметра 
провода. Случайные погрешности ведут к «рассеиванию» величины со- 
противления. 
Если представить допуск сопротивления обмотки как АЁ, то со- 
вокупность всех погрешностей, которые возникают в процессе намотки, 


(8) 


п 
должна быть меньше AR, то есть AR > > ÒR; . 
i=l 

Выводы. 

В ходе данной работы провед ен анализ факторов, влияющих на 
качество изготовления обмоток электрических машин. Установлено, 
что основным параметром, непосредственно влияющим на качество из- 
готовления обмоток электрических машин, является усилие натяжения 
провода при намотке. Параметр усилия натяжения провода является си- 
стематически изменяемым и требует уч ета при определении длины про- 
вода и числа витков обмотки. 
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INVESTIGATION OF THE INFLUENCE OF THE TENSION 
FORCE OF THE WINDING WIRE ON ELECTROMAGNETIC 
CHARACTERISTICS IN THE MANUFACTURE 

OF WINDINGS OF ELECTRIC MACHINE 


In this paper we study the influence of the tension force of the winding wire 
on the electromagnetic characteristics in the manufacture of windings of elec- 
trical machines. The dependence of the resistance of the winding and its induct- 
ance on the magnitude of the tension force of the wire is determined. The influ- 
ence of the error of the wire diameter and the error of the specific resistivity of 
the winding wire is obtained. 


Keywords: electrical machine, electrical machine winding, winding wire, 
conductor resistance. 
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УДК.678.019.13 
С.А. Бысов, И.Э. Гинсар 


ИССЛЕДОВАНИЕ ПРОЛИВАЕМОСТИ ФОРМЫ 
В ЗАВИСИМОСТИ ОТ ОПТИМИЗАЦИИ РАЗМЕРОВ 
ЗАГОТОВОК И ЛИТЕЙНОЙ СИСТЕМЫ 


В данной статье приведен имитационный эксперимент, направ- 
ленный на определение влияния диаметра центрального литникового 
канала и толщины стенки отливаемый детали на проливаемость 
пресс-формы. Эксперимент выполняется при использовании Solid- 
Works Plastics. Приводятся этапы проведения эксперимента. Пред- 
ставлена литниковая система, состоящая из пластмассовых дета- 
лей «Корпус» и «Крышка». Обоснована необходимость совместного 
изготовления деталей «Крышка» и «Корпус». Рассмотрены все соче- 
тания изменяемых размеров и их графическое представление. На ос- 
новании полученных результатов проводится регрессионный анализ 
для определения, какой из выбранных параметров наиболее влияет на 
проливаемость пресс-формы. Проводится имитационный процесс 
литья под давлением пластмассовых деталей с выбранными разме- 
рами. Приводятся результаты в виде следующих эпюр: «Давление по 
завершению процесса заполнения», «Центральная температура 
фронта процесса», «Давление на входе». 


Ключевые слова: литьё пластмасс под давлением, программное 
обеспечение, имитационное моделирование, центральный литнико- 
вый канал. 


Введение. Целью данного эксперимента является определение оп- 
тимальной толщины детали и диаметра центрального литникового ка- 
нала для обеспечения идеальной схемы проливаемости. Под идеальной 
проливаемостью понимается достижение полной заполняемости пресс- 
формы [1]. Расчетным методом было определено, проливаемость обес- 
печивается при диаметре центрального литника равным 1 мм. Но необ- 
ходимо определить размеры, при которых заполнение формы обеспечи- 
вается полностью с минимальным количеством дефектов. То есть экс- 
перимент направлен на сокращение количества дефектов литья в ре- 
зультате оптимизации размеров литых заготовок и литниковой системы. 

Так как деталь «Корпус» используется в автомобилестроении для 
монтажа реле различных типов, необходима хорошая собираемость с 
ответной деталью «Крышка» [2]. Поэтому детали отливаются в одной 


пресс-форме, объём которой равен 99 см?. Литниковая система с дета- 
лями представлена на рисунке 1. 


Рубрика «Машиностроение» (ISSN 2413-6220) 20 


Рис. 1. Литниковая система с деталями 


Эксперимент выполняется с использованием программного обес- 
печения SolidWorks Plastics. По предварительным экспериментам де- 
таль «Крышка» будет проливаться в любом случае, поэтому будет изме- 
няться толщина стенки детали «Корпус». Выбранные размеры для про- 
ведения эксперимента: 


1. Диаметр центрального литникового канала 4, =1мм и 
d, =Змм. 

2. Толщина стенки отливаемой детали b = 0,8 мм ир=12 мм. 

На основе полученных результатов проводится регрессионный 
анализ для определения 2 факторов. Также будет известно, какой 
именно размер в наибольшей степени влияет на проливаемость. 

Имитационный компьютерный эксперимент содержит 4 этапа, со- 
стоящих из различных сочетаний размеров. 

Та, 23MM и Б =0,8 мм; 

2. d, =1мм и Б=1,2 мм; 

3. d, =1мм и = 0,8 мм; 

4. d, =Змм и b =1,2 мм. 

При толщине детали 0,8 мм и диаметре литника 3 мм не достига- 
ется идеальной схемы проливаемости. Средний объём расплавленной 
пластмассы в этом случае равен 75,9 см“, который составляет 76 Фо от 
общего объема. Результат проведения 1 этапа представлен на рисунке 2. 
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Рис. 2. Проведение | зтапа 


При размерах d, =1мми b—1,2MM также расплавленная пласт- 


масса полностью не заполняет пресс-форму. При зтом по сравнению с 
предыдущим этапом проливаемость намного хуже. Средний объем со- 


ставляет 75,9 см? или 56 %. Результат представлен на рисунке 3. 


Рис. 3. Проведение 2 этапа 


При размерах d, =1мм и b —0,8 мм также не достигается запол- 
нение пресс-формы полностью (рисунок 4). Средний объем расплавлен- 
ной пластмассы равен 88,8 см“, что составляет 89 % от общего объема. 
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Рис. 4. Проведение 3 этапа 


И последний этап, в котором размеры диаметра центрального лит- 
никового канала и толщины стенки детали соответственно составляют 
3 мм и 1,2 мм. При таких размерах достигается идеальная схема проли- 
ваемости, при которой объем впрыска равен объему литниковой си- 
стемы с отливаемыми деталями. Результат проведения последнего этапа 
представлен на рисунке 5. 


Рис. 5. Проведение 4 этапа 


Именно этот случай будет использован для проведения регресси- 
онного анализа. В таблице 1 представлены исходные данные. 
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Таблица 1 
Исходные данные 
Факторы i ы i Интервал варьирования, Ах, 
+ E 
X, (4, ) 3 2 1 1 
Х„(Ь) 0,8 1 1,2 0,2 


Нормализованное значение уровней факторов связано с натураль- 
ными следующими зависимостями: 


а, – 2 
м = F ; (1) 
b-1 
X, == 2 
? 0,2 (2) 


План эксперимента в виде расширенной матрицы планирования 
представлен в таблице 2. 


Таблица 2 
Расширенная матрица планирования 
N Хо X, X, ХХ, 
1 + + + + 
2 + - + - 
3 + + - - 
4 + - - + 


План эксперимента с учётом результатов измерений объема расплав- 
ленной пластмассы, поступающей в пресс-форму, представлен в табл. 3. 


Таблица 3 
План эксперимента с учётом результатов измерений 
Результаты опыта 
1 2 3 
1 3 1,2 99,3 99,3 99,4 
2 3 0,8 71,3 71 71,1 
3 1 1,2 56,4 56,1 55,9 
4 1 0,8 75,6 75,4 75,4 


Обработка результатов проводится по следующим этапам: 
1. Для каждой строки матрицы определяется среднее значение от- 
клика по следующей формуле: 


г-н. (3) 
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2. Для каждого отклика вычисляется дисперсия по формуле: 


nrc 2 
wo 
ры (4) 
п-1 
3. Проверка гипотезы об однородности дисперсий по С-критерию 
проводится по следующей формуле: 


G, = ак. (5) 
2,5; 
i-l 


Затем расчетное значение критерия сравнивается с табличным, ко- 
торый определяется в зависимости от принятой доверительной вероят- 
ности и числа степеней свободы. 

В данном случае с доверительной вероятностью 95 % и числом 
степеней свободы равным 2 табличное значение С-критерия равняется 
0,7679. Так как С, < С а» гипотеза об однородности дисперсий прини- 


мается с доверительной вероятностью 95 %. 
4. Вычисляется дисперсия воспроизводимости по формуле: 


N 

Ys 
у EL. 6 
e (6) 


5. По результатам эксперимента вычисляются коэффициенты 
уравнения регрессионного анализа по формуле: 


N — 
>. 
p, = 1. 


1 


таб? 


(7) 
N 


При вычислении коэффициента B каждом опыте эксперимента 
необходимо приписать знаки в соответствии с вектор-столбцом. 

6. Проверка гипотезы о статистической значимости коэффициента 
уравнения регрессии. Вычисляется дисперсия для каждого коэффици- 
ента по формуле: 


52 
Soc 8 
в N-n (8) 
Вычисление доверительного интервала по формуле: 
А, = Гор Sp. (9) 


7. Сравнение каждого из ранее вычисленных коэффициентов с до- 
верительным интервалом. Если значение коэффициента больше значе- 
ния доверительного интервала, то он статистически значим с довери- 
тельной вероятностью 95 %. В противном случае коэффициент стати- 
стически незначим, и он может быть отброшен [3]. 
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После проведения всех зтапов регрессионного анализа было полу- 
чено следующее уравнение: 


y(V)=80,02+9,2- X, +30,3- X, +56,1.Х,-Х.. (10) 


Для применения уравнения в практических задачах уровни факто- 
ров необходимо представить в натуральном масштабе. Полученное 
уравнение имеет следующий вид: 


y(V) = 80,02 +9,2. du нюз (52 sea. m poner . (10) 
2 2 2 0,2 

Из полученного уравнения можно сделать следующие выводы: 

1. Влияние диаметра центрального литникового канала на про- 
ливаемость формы в 3 раза меньше по сравнению с толщиной отли- 
ваемой детали. 

2. При увеличении толщины стенки детали увеличивается и объем 
расплавленной пластмассы, поступающей в пресс-форму. Эффект взаимо- 
действия значителен по сравнению с влиянием факторов по отдельности. 

После определения оптимальных значений размеров проводится 
имитационный процесс литья под давлением. После проведения всех 
этапов компьютерного анализа получаем детали, не имеющие каких- 
либо дефектов [4]. Результаты проведения имитационного процесса ли- 
тья могут быть представлены в виде эпюры «Давление по завершению 
процесса» (рисунок 6), эпюры «Центральная температура фронта про- 
цесса» (рисунок 7), диаграммы «Давление на входе» (рисунок 8). 


E) 


' | 
Рис. 6. Эпюра «Давление по завершению процесса» 


е? 


Рис. 7. Зпюра «Центральная температура фронта процесса» 
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[ES HMP 0 ,8 По умолчанию : Давление на входе] 
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Рис. 8. Давление на входе 
Выводы. 


1. Проведен имитационный эксперимент при использовании про- 


граммного обеспечения. 


2. Проведен регрессионный анализ и представлены все его этапы. 
3. Определены степени влияния диаметра центрального литнико- 


вого канала и толщины стенки отливаемой детали «Корпус» на проли- 
ваемость формы. 


4. При определении оптимальных размеров проведен имитацион- 


ный процесс литья под давлением для подтверждения отсутствия каких- 
либо дефектов. 


[1] 


[2] 


[3] 
[4] 
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THE RESEARCH OF MOLD SPILLAGE DEPENDING 
ON THE OPTIMIZATION OF THE SIZES 
OF WORKPIECES AND GATED PATTERN 


In the article we conduct the simulation experiment aimed at determining 
the influence of the diameter of the central gating channel and the thickness of 
the wall of the cast part on the mold spillage. The experiment is performed using 
SolidWorks Plastics. The stages of the experiment are given. The gating system 
consisting of the parts «House» and «Lid» is presented. The necessity of simul- 
taneous production of the parts. All combinations of resizable sizes and their 
graphical representation are considered. Based on the results, a regression 
analysis is performed to determine which of the selected parameters most af- 
fects the spillage of the mold. The simulation process of injection molding of 
plastic parts with the selected dimensions is carried out. The results are pre- 
sented in the form of the following diagrams: "Pressure at the end of the filling 
process", "Central temperature of the process front", "Inlet pressure". 


Key words: injection molding of plastic, software, simulation modeling, cen- 
tral gating channel. 
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УДК 621.9 
Н.В. Зенкин, К.А. Шулепов 


МОДЕЛИРОВАНИЕ ЗАТРАТ ВРЕМЕНИ 
НА НАЛАДКУ И СМЕНУ ИНСТРУМЕНТА 


Статья посвящена исследованию затрат времени на наладку и 
смену инструмента, как составляющей штучного времени. Исследо- 
вана зависимость времени обслуживания рабочего места, в частно- 
сти замены режущих пластин во фрезах, установки осевого инстру- 
мента в патрон и настройки инструмента на размер на фрезерном 
станке с ЧПУ, от количества повторений данных переходов, выпол- 
няемых рабочим. Получены математические модели данных процес- 
сов в виде уравнений регрессии. Построены линии тренда и козффи- 
циенты достоверности аппроксимации с помощью пакета приклад- 
ных программ МАТГАВ. 


Ключевые слова: способ наладки, время на смену инструментов, 
нормирование. 


Интерес к вопросам оптимизации штучного времени достаточно 
высок и постоянно продолжает расти. Существуют различные пути ре- 
шения данного вопроса [1-6]. Однако следует отметить, что и время об- 
служивания рабочего места может занимать немалую часть штучного 
времени [7]. С целью исследования и моделирования процессов затрат 
времени на смену и наладку инструмента, в рамках выполнения научно- 
исследовательской работы, был проведен эксперимент. В данном экспе- 
рименте моделировался процесс смены и наладки инструмента на раз- 
мер на станке Mini Mill 450. 

Измерению подлежали: процесс замены твердосплавных пла- 
стин на фрезах R390-020A20-11M (рис. 1, рис. 2), ISCAR HCE 012- 
A-L120-C12 (рис. 3), GARANT 215 157 32/3 (рис. 4), процесс замены 
осевого инструмента в патроне, в частности: сверла ВТ-40-ЕВ 16-100 
(рис. 5), фрезы BT-40-ER32-200 (рис. 6) и фрезы BT-40-SLA20-90 
(рис. 7), а также процесс наладки инструмента на размер на станке 
Mini Mill 450 (рис. 8). 


wee | 
.1 20А2О-ЛТМ 
PR - 775 


Рис. 1. Фреза R390-020A20-11M 
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Рис. 3. Фреза ISCAR НСЕ 012-A-L120-C12 


Рис. 4. Фреза САКАМТ 215 157 32/3 


Рис. 5. Сверло ВТ-40-ЕК16-100 в патроне (6) ив разобранном виде (а) 
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Рис. 6. Фреза ВТ-40-ЕКЗ2-200 в патроне (6) ив разобранном виде (а) 


Рис. 7. Фреза BT-40-SLA20-90 в патроне (б) ив разобранном виде (a) 
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Рис. 8. Процесс наладки инструмента на размер на станке Mini Mill 450 


Для построения математической модели (уравнения регрессии) — 
был использован МНК (метод наименьших квадратов). Метод наимень- 
ших квадратов дает состоятельные, эффективные и несмещенные 
оценки параметров уравнения регрессии [8]. 

Формально критерий МНК представляет собой следующее выра- 
жение (1): 


"nv у 
5=> (у-у) — —) min . (1) 
Система уравнений для определения коэффициентов регрессии (2, 3): 
bnth: M х=у; (2) 


2 
by х+5- У х =} yu. (3) 
С учетом полученных данных система уравнений примет следую- 
щий вид (4, 5): 


25-b) +325-Ь = 694,5; (4) 
325. b, 4 5525. b, = 8195,1. (5) 
Решив систему уравнений, найдем (6, 7): 
1300-5, =-833,4. (6) 
Откуда 
b, =-0, 6411. (7) 


Для нахождения коэффициента Бо подставим полученное значение 
в систему уравнений (8-11): 


25-b) 4 325: b, = 694,5; (8) 
25. b, +325- (—0,6411) = 694,5 ; (9) 
25. b, —902,85; (10) 


by 236,114. (11) 
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Эмпирические коэффициенты регрессии: b, =-0, 6411, 
Бу = 36,114. С учетом коэффициентов уравнение регрессии будет иметь 
следующий вид (12): 
у=-0, 6411.х+36,114. (12) 
Для остальных экспериментов математические модели имеют вид: 
* Замена пластин на фрезе ISCAR HCE 012-A-L120-C12 (13): 


у=-0,3185.х+ 28, 784. (13) 
* Замена пластин на фрезе GARANT 215 157 32/3 (14): 

у=-0,3992.х+32,409. (14) 
* Замена осевого инструмента в патроне (15-17): 

у=-1,241-х+ 61,837 ; (15) 

у=-0,3672-х+52,218; (16) 

у=-0,8231.х+48,724, (17) 
* Процесс наладки инструмента на размер на станке (18): 

у=-1,2769.х+80,48. (18) 


Линии тренда (рис. 9-15) и коэффициент достоверности аппрокси- 
мации получены с помощью пакета прикладных программ МАТГ АВ. 


[12345678910 11 12 13 14 15 16 17.. 

18 19 20 21 22 23 24 25]; 

у = [49.9 31.3 32,3 30.2 34.5 29.7 24.3 28.7 2842. 
35.4 25.60 25.9 24 24.5 23.8 25.8 22.5 33.9 23.1. 
22.4 23.4 21.1 24.3 21.1 22.60]; 


grid on 
plot(x, y, '*r'); 
xi = min(x):1:max(x); 
N = 4; 
coeffl = polyfit(x, у, N); 
у2 = 1; 
for k-0:N 
y2 = y2 + coeffl(N-k41) * xi.^k; 
end 
hold on; plot(í(xi, y2; 'r'); grid on; 
М = 21; 


coeff2 = polyfit(x, у, N); 

std (у- (соеЕЕТ (1) *х+соеЕЕТ (2))); 

ста (у- (соеЕЕ2 (1) *х.^2+соеЕЕ? (2) *х+соеЕЕ? (3))); 
Узк=зим (у) ./25; 

Взаг2-1- (зам ( (у4-у2) . ^2) ./sum((y4-Ysr) .^2)) 
у4=-0.6411.*х+36.114; 

figure (2); ріої (х,у, 'm'); grid on; 
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Рис. 9. Замена пластин на фрезе R390—020A20-11M 
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Рис. 15. Процесс наладки инструмента на размер на станке 
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В результате эксперимента были получены математические мо- 


дели процесса смены и наладки инструмента на размер, в виде уравне- 
ний регрессии, по ним построены линии тренда, которые позволяют 
проанализировать характер поведения кривых времени при дальнейшем 
увеличении повторений процесса. 
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SIMULATION OF TIME FOR ADJUSTMENT 
AND CHANGE OF TOOLS 


The article is devoted to the study of the time spent on setting up and chang- 
ing a tool as a component of piece time. The dependence of the service time of 
the workplace, in particular the replacement of cutting inserts in milling cutters, 
the installation of an axial tool in a chuck and the tool setting to size on a CNC 
milling machine, on the number of repetitions of these transitions performed by 
the worker, is investigated. Mathematical models of these processes are ob- 
tained in the form of regression equations. Trend lines and approximation con- 
fidence coefficients are constructed using the MATLAB application package. 


Keywords: setup method, tool change time, rationing. 
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УДК 62-1 
И.А. Зенкина 


НАХОЖДЕНИЕ ГЛАВНОГО МОМЕНТА 

СИЛ СОПРОТИВЛЕНИЯ 

В ГАЗОДИНАМИЧЕСКОМ ПОДШИПНИКЕ, 
ПРОФИЛИРОВАННОМ СПИРАЛЬНЫМИ КАНАВКАМИ 


Статья посвящена процессам, протекающим в газодинамическом 
подшипнике со спиральными канавками, у которого одна из поверх- 
ностей профилирована спиральными микроканавками, а другая явля- 
ется гладкой. Рассматривается косоугольная система координат, 
которая представлена логарифмическими спиралями и окружно- 
стями. В этой системе координат записываются уравнения, опреде- 
ляющие изменение давления в активной зоне. Для каждой пары пере- 
мычка-канавка найдены силы вязкого трения и их главный момент 
сопротивления относительно оси подшипника. Получено уравнение 
для главного момента сил сопротивления в активной и гладкой зоне. 
Итогом работы является алгоритм, для вычисления главного мо- 
мента сил вязкого трения, развивающегося в слое смазки газодина- 
мического подшипника, профилированного спиральными канавками. 


Ключевые слова: газодинамический подшипник, спиральные ка- 
навки, косоугольные координаты, момент сопротивления. 


Введение. Газодинамические подшипники активно используются 
в технике. К их неоспоримым достоинствам можно отнести экологич- 
ность, уверенную работу при различных температурах, увеличение ско- 
рости ротора в результате низкой вязкости газов. Рассматривается газо- 
динамический подшипник, у которого рабочая поверхность состоит из 
двух близко расположенных поверхностей, и одна из них профилиро- 
вана спиральными микроканавками. Целью работы является вычисле- 
ние момента сил вязкого трения, которые развиваются в слое смазки га- 
зодинамического подшипника со спиральными канавками. 

Косоугольная система координат. Рассмотрим плоский газодина- 
мический подшипник, состоящий из двух поверхностей, гладкая по- 
верхность неподвижна, а поверхность с нанесенными спиральными ка- 
навками вращается [1] (рис. 1). 


Рис. 1. Газодинамический подшипник со спиральными канавками 
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Рабочая зона парциального подшипника представлена парами ка- 
навка-перемычка, и распределение давлений в каждой паре одинаково. 
Одна пара изображена на рис. 2, где О — некоторая точка на фронтальной 
кромке, которой канавка отделяется от перемычки. 


Рис. 2. Пара канавка-перемычка 


Введем ортогональную систему координат у, 9 


x-(r-n)/R: 9-9-фо, (1) 
здесь ђ и фу — полярные координаты точки О. Обозначим 
рег Е. (2) 


Рассмотрим еще криволинейную косоугольную систему коорди- 
нат Е, п причем линии 6 = const являются окружностями с центром в 


О, ‚а ЛИНИИ у = const являются логарифмическими спиралями с полю- 


COM в О, и углом атаки ҷу [3]. 


А А 
= i =п а 4 94 Хову ЕЕ 1 Хиву, (3) 
фу +p 2n 2m p 2n p 
здесь п — число спиральных канавок. 


Значения 1] на границах пары можно определить через параметры 
«__% а 

+o фу +ф, 

Вычисление сил вязкого трения. Рассмотрим полосу смазочного 


слоя (рис. 2) шириной dr , расположенную вдоль линии г = const. Вы- 
делим элементарный участок 45 этой площадки 


n 


„к+а=1. (4) 


2 
ds = == К?рарап, (5) 
п 


Если площадка ds расположена на высоте = ћ от плоскости пере- 
городок нижней поверхности, имеющей канавки, то проекция силы вязкого 
трения, приложенной к этой злементарной площадке, на ось 9 равна 


ак, = (Ze ds, (6) 
dz h 
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здесь и — динамический коэффициент вязкости, а У, - проекция скоро- 


сти газовых частиц на ось 9. Исследуем силы трения в окрестности ка- 
навки. На высоте z от нижней стенки смазочного слоя проекция скоро- 
сти на ось 9 определяется выражением 


ска (Весна) op 
c+h Фик(р+х) 09. 


Здесь p - давление, с — глубина спиральных канавок, h — рабочий 


Vs = oR(p- x) 


(7) 


зазор, < — угловая скорость подвижной детали. 
Продифференцировав равенство (7), получим 
ду, oR(p*x) h-c-2z др 
Oz һ+с 2иК (р+х) 09. 
Полоса содержит точку О, где у =0 (рис. 2). Тогда из выражения 


(8) 


(8), учитывая, что у=0, z = 1, получим 


ЕЗ .9Rp, htc др. (9) 
Oz Јр h«c 2uRp 29 
Подставив (9) и (5) в (6), найдем 
210 р др 
ак -TR +(h+ dpd 10 
9 Е- me К?р? +(һ+с) 2| pdn. (10) 


Данное уравнение получено для канавки, позтому вместо перемен- 
ной р логично взять p, . Для перемычки надо принять с = 0 изаменить р 


на ро. Таким образом, для канавки и перемычки соответственно получим: 


dRg, = л п 222 р + (n+) арат 
h+c 09 (11) 
2р0 Һ Z) 
dR =- R ^R dpd 
52 Ex y Ree Hen 
C учетом (3), получим 
Oy Hs ОП ара ОП, Сва Не. GOES. о) 


09 2P 09 4nP да 09 2P 09 4nP Om 
здесь П, и IL, — квадраты безразмерного давления B области канавки и пе- 
ремычки, p, — атмосферное давление, Р = р/р, — безразмерное давление. 


Квадраты безразмерного давления можно аппроксимировать дву- 
мерными сплайнами — кубическими по тү и линейными по Е [2,4,5]: 


Tt 
IT, =П+- (аб +am+am? +a"); 

á (13) 
IL, -II4 nee от + Бул? +). 


Здесь П - квадрат а давления в точке О. 
Продифференцировав равенства (13) по п, получим: 
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ап, 


On 
Выражения (13) и (14) позволяют привести соотношения (11) к 
следующему виду: 


b *2b*3ba). (14) 


= MC +2ауп+ Зап.) : E - | 


T 2uoR? 2 hcc 2 
аб р, | И p. & —— (a, + дар + Зал) |dp dy; 
ae д6 АР ( 1 21 3 ) pen 
(15) 
üh cR ВИК з (bi 25m + Зл: ) арап 
92 a Pah 4P 1 2 3 у 


Расчет главного момента сопротивления в активном слое. Найдем 
моменты сил трения относительно оси подшипника, для чего надо умно- 
жить dRo И аКу на г = Rp 

2uoR? у һ+с 
Ра (^ t c) i 4P 


dm, "fen, [eset ehe?) papam 
n 


(16) 
2uoR? 2 h 


+ 
Рай Е 4Р 


Пусть dm - момент сопротивления полоски шириной dr = Кар. 


dm, -Ertz А Барат 
п 


Она пересекает пару канавка-перемычка по линии р = const (рис. 2). 
1-0 п=а 
dm= | ат (п) + | ат, (п). (17) 
п=-к 1-0 
Обозначим за dM, момент сил сопротивления, приложенных к 
кольцевому фрагменту непрофилированной детали, шириной dr = Кар, 


относительно оси подшипника. Так как dM, в п раз больше ат ‚то: 


п-0 n=a 
dM, =n; | dm(n)+ | tnl. (18) 
T--K т=0 

Выполним вспомогательные расчеты: 
0 0 0 
[а= к; 2 | ndn- —к°; 3f па = кз; 
-к -к -к (19) 
а а а 
[4n =q; 2fn ат = а?; Зл ал = а3. 
0 0 0 
Полученные соотношения (19) и (16) позволяют привести выраже- 
ния (18) квиду 


2uoR? h+ 
am, = akèn] = е | као ка} 
Ра +с 


(20) 
һ 
ric + ab, + ab, Ji» dp. 
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Обозначим за Л число сжимаемости и за у безразмерный пара- 
метр 


и брок? т һ 
pak? | hec. 


Учитывая эти выражения, запишем 


2uoR? 2uoR? 
ES | È нар. ES нак 


А 


(21) 


је = 1 Аһ(о.+ку}ўр?; 


p, “А+с h рай? һ+с 3 
һ+с h 
тр (ка = 5a) + ка, )+ (nbi e b, наь) = 
h | к(А+с 
-+ ш Да ssa enn eb наь). (22) 


h K 2 2 
= Б] (а – ка, +K аз) + a (b, + ар, +a s)|; 
Используя соотношения (22) dM, можно привести K виду 
А 1 
dM, = IC + ку)р? + б -= Ka, + E 
+ам (b, +ab, + «^b, р dp. 
Выражение для главного момента M, сил трения, приложенных K 


непрофилированной детали со стороны смазочного слоя, относительно 
оси подшипника можно получить, что интегрируя dM, (23) oT p =p; до 


p-l 


1 
М, = черан) [к(а – ка +к?аз)+ 
1 


(24) 
+ам (b, +ab, + ob, јр ар. 
Выражение (24) преобразуется следующим образом 
ж А QU dd 2 
M,-nR' ph — (o ку)(1 ри) + | к(а ка, + к аз )+ 


+ ОУ (b, t ab, + ob, JJa}. 


Введем в рассмотрение безразмерный момент сопротивления М К 
2 * 

M, =тК рМ). (26) 

Обозначим для безразмерный зазор 6, =ӧ&; и нормированное 


число сжимаемости A, (где  — номинальный зазор) 


и=1+6; = = ; А, =Ли?. (27) 
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Вьражение М | для активной зоны подшипника примет вид 


1 
M, = “(а +ку)(1-р!)+-= f [s (s -ка, +к2аз)+ 
Pi (28) 
+ам (b + ор, + ab, )р ар. 


Коэффициенты сплайнов а;, 5, (13), выражаются через Р и 4Р/др, 


исходя из гидродинамической эквивалентности всех пар канавка-пере- 
мычка и из условия неразрывности локальных массовых расходов газа 
[1]. В результате преобразований получим следующее уравнение для 
доминирующего давления [4, 6, 7]: 


vo =1–ур а 21-1; y = yo/(Le vC); v 21- yi B у"; 
u - 1G A = Аи; Л, = (xA /nP)p? sin? y; e = 5, sin2w; 
2(a + кВ) 
В(1+ кд) 
9, =а(3+од); оу = 4в(а? —^.); с, =4(1-В); по = 35 * B9; 


лите ар ЕЕ 


9, . 
& = (9,6; cos^ y - Bo; )/n, 17 ; 2y; 
0 


E3 = 06, + 4e; £, = 06; A = 3(1- кб); 8 = 20/A; 6, = 864; 
E6 = беј: А = A(k ap); A, =к(1-В); (29) 
W = А +2. (-61 + ква - к26,); 
Фу = 2Луу?д (Е - кб, + ко јр /W o; - 28 и WpP)); 


—— =-Ф,-Ф„О. 
др 


с 
Здесь уа = n — нормированная глубина канавок. 
+c 


0 
* ~ 
Окончательная формула М, для активной зоны подшипника при- 


мет вид 


* ГА, 4 ку Ло | 3 
M, = —— (а + кУ )|1— +—4— yv | фр ар+ 
1715 2(9+ку)( Bs cec pup 
' (30) 
1 
и ар 2 
+— | 0, —p°dp тезу; 
| ias р (18 


Pi 


éo = k(1- кв); Ө, = 563 – 51; 0; = 6, - 6054. (31) 
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Для нахождения безразмерного момента сопротивления активной 
зоны надо решить уравнение (29), используя численные методы, потом 


* 

определить вспомогательные переменные (31), и далее вычислить М, 
по формуле (30). 

Расчет главного момента сопротивления для гладкой зоны. Непрофи- 

лированная зона имеет форму кольца с радиусами R, и R, толщиной A. 


Внизу она ограничена неподвижным диском, а верхний диск вращается с 
угловой скоростью о. Безразмерный момент сил сопротивления гладкой 


зоны М Х можно получить, сделав допущения отсутствия канавок. 
в Пе x 
M =. 20р. -р2). (32) 


Фактический момент сил сопротивления для гладкой зоны связан 
с безразмерным таким же образом, как и для активной зоны (26) 


2 ж 
Расчет главного момента сил сопротивления в смазочном слое. 


Окончательно безразмерный момент сил сопротивления М смазоч- 
ного слоя получается в результате сложения безразмерных моментов со- 


противления активной М 1 (30) и гладкой (32) М 2 30H. 


« ТА 
M = — = (a-«xv)(1-pt)* pt -p3 |+ 
(36) 
Ky | Ло р 3 и | ар > 
+— < — үу | фр'ар+—– | Ө, ——р др го. 
Iz Ка А 


Фактический момент сопротивления М определяется равенством 
2: ж 

М = лК pM . (37) 
Выводы. В данной статье было получен алгоритм расчета глав- 
ного момента сил вязкого трения, создающегося в смазочном слое газо- 
динамического подшипника, одна из двух поверхностей которого про- 
филирована спиральными канавками. Для нахождения момента была 
применена косоугольная система координат, из окружностей и логариф- 

мических спиралей. 
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FINDING OF MAIN MOMENT OF RESISTANCE FORCES 
IN GAS-DYNAMIC BEARING PROFILED 
WITH SPIRAL GROOVES 


The paper focuses on processes taking place in a gasdynamic bearing with 
spiral grooves, in which one of the surfaces is profiled by spiral microchannels 
and the other is smooth. The oblique coordinate system is considered, which is 
represented by logarithmic spirals and circles. This coordinate system records 
equations that determine the change in pressure in the core. Viscous friction 
forces and their main moment of resistance relative to bearing axis are found 
for each pair of bridge-groove. The equation for the main moment of the re- 
sistance forces in the active and smooth zone is obtained. The result of the work 
is an algorithm for calculating the main moment of viscous friction forces de- 
veloping in the lubricant layer of a gasdynamic bearing profiled by spiral 
grooves. 


Keywords: gas-dynamic bearing, spiral grooves, oblique coordinates, re- 
sistance moment. 
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УДК 062 
Е.Н. Ильдияров, В.В. Коротков 


ОБЗОР И АНАЛИЗ СУЩЕСТВУЮЩИХ МЕТОДОВ БОРЬБЫ 
С ВИХРЕОБРАЗОВАНИЕМ В ЦЕНТРОБЕЖНЫХ НАСОСАХ 


Проводится анализ источников возникновения гидродинамических 
вибраций в центробежных насосах. Особый интерес в данной работе 
представляет возникновение обратных токов, вихреобразование и 
неравномерное обтекание элементов проточной части. Представ- 
лена картина течения жидкости в центробежном насосе при нали- 
чии на входе обратных токов. В настоящей работе приведены суще- 
ствующие методы уменьшения негативного влияния обратных то- 
ков на вибрационные характеристики центробежного насоса. 


Ключевые слова: гидродинамическая вибрация, пульсация давле- 
ния, центробежный насос, обратные токи, повышение эффективно- 
сти насоса. 


Введение. При работе осевой ступени насоса на подаче меньше 
50 % от оптимальной, во входной части осевого колеса, образуются 
вихри (вихревые течения) жидкости, усиливающиеся с уменьшением 
подачи. Вихри приводят к повышенной вибрации, неустойчивой работе, 
генерации низкочастотных пульсаций давления, имеющие характер ав- 
токолебаний. Возникшие автоколебания достигают большой ампли- 
туды и могут привести к выходу насоса из строя [1]. 

Рассмотрение вопросов с генерацией пульсаций давления и вибра- 
ций в насосах следует приступать с общей классификации источников 
возникновения вибраций. Выделяют две основные группы источников 
вибрации: 

ф гидродинамические; 

+ механические. 

К источникам гидродинамической вибрации относятся: 

1. Кавитационные процессы. 

2. Возникновение обратных токов, вихреобразование и неравно- 
мерное обтекание элементов проточной части насоса. 

3. Пульсации давления после рабочего колеса насоса. 

К механическим источникам вибрации относятся: 

1. Дисбаланс ротора. 

2. Некомпенсированные центробежные силы. 

3. Работа привода насосного агрегата. 

4. Влияние подшипниковых опор. 

Ведя речь о самих источниках вибрации, анализ научно-исследо- 
вательских работ по исследованию их влияния на вибрацию в насосах 
необходимо проводить отдельно для каждого источника. Такой метод 
дает возможность оценить имеющийся способ снижения вибрационной 
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нагрузки путем снижения воздействия данного источника вибрации или 
полной его нейтрализации. 

Особый интерес представляет источник гидродинамической виб- 
рации - возникновение обратных токов, вихреобразование и неравно- 
мерное обтекание злементов проточной части. 

В ходе проектирования и эксплуатации насосных агрегатов часто 
фиксировалось появление в их проточной части сложных вихревых те- 
чений. Природа таких течений, как и причины, их появления, могут 
быть неодинаковыми. 

Опираясь на анализ литературы следует выделить основные при- 
чины возникновения вихревых течений: 

+ отрывы потока жидкости при обтекании элементов проточной 

части насоса; 

* перетекание жидкости через зазоры из зоны высокого давления 

в зону низкого давления; 

* взаимодействие ротора насоса с ero статором; 

* закрутка потока на входе в рабочее колесо насоса; 

Ф работа насосов на режимах сильно отличающихся от номи- 

нального. 

Нужно отметить, что появление вихревых структур течения обычно 
сопровождается высокими потерями полного давления в рабочей жидко- 
сти, высокими уровнями пульсаций давления, локальным снижением ста- 
тического давления и вибрацией частей насосного агрегата [2]. 

Именно поэтому появилась необходимость в проведении новых 
исследований направленных на определение путей снижения неодно- 
родности в потоке гидравлической жидкости. Эти исследования прово- 
дятся как для насосов общей промышленности, так и для насосов си- 
стемы подачи авиационных двигателей и ЖРД. 

Имеющиеся методики борьбы со сложными вихревыми структу- 
рами течения в большинстве случаев состоят из нескольких этапов: 

+ проведение натурного или численного эксперимента для обна- 

ружения неоднородностей в рабочей жидкости; 

* выявление источников или причин возникновения вихревых те- 

чений в жидкости; 

Ф проработка изменений, которые следует внести в конструкцию 

насосного агрегата с целью достижения однородности потока; 

* проведение дополнительных экспериментальных или числен- 

ных исследований для оценки эффективности проведенных 
конструктивных изменений. 

Данный подход позволяет выработать рекомендации по выравни- 
ванию потока рабочей жидкости в проточной части насоса, а именно, 
предотвращению возникновения в нем вихревых структур и неоднород- 
ностей полей скоростей и давлений [3]. 

В работе [6] с помощью гибких нитей исследовалось относитель- 
ное движение жидкости непосредственно перед шнековым преднасосом 
и в самом шнеке. Исследования проводились на реконструированном 
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авиационном насосе ЗЦН-23. Наиболее ярко все особенности располо- 
жения гибких нитей в каналах шнекового преднасоса проявились на рас- 
ходах, близких к нулевым. На основании анализа расположения гибких 
нитей во входной области шнека была построена картина течения жид- 
кости в шнековом преднасосе при наличии на входе обратных токов, 
представленная на рис. 1. 


Рис. 1. Картина течения во входной части шнекового преднасоса: 
1 — точка торможения потока; 2 — источник; 3, 5 — циркуляционные потоки; 
4 – обратные токи на входе; 6 – шнек; 7 — активный поток; 

8 — обратные токи на выходе: 

a — периферийное сечение; б — меридиональное сечение 


Как видно из рисунка обратные токи имеют вихревой характер. 
Нестационарный отрыв потока с поверхностей межлопастного канала 
приводит к появлению значительных пульсаций давления на его стен- 
ках, что приводит к появлению вибрации шнека, снижая его эксплуата- 
ционную надежность [5]. Пульсации давления, возникающие в зоне от- 
рыва, передаются как против потока, так и по потоку к центробежному 
колесу. Эти пульсации давления, являясь одним из источников шума, 
вызывают дополнительные нестационарные силы, действующие на эле- 
менты насоса, что отрицательно сказывается на работе всего насоса. 

В работе [4] рассматривается пример возникновения сложных вих- 
ревых структур на входе в осевую ступень (шнек) насоса системы по- 
дачи авиационного двигателя. Исследования проводились по выявле- 
нию причин появления неоднородности потока на входе в насос. В ре- 
зультате чего источник был найден. В конструкциях насосов, подобных 
рассматриваемой в работе [4], на входе в насос устанавливается направ- 
ляющий аппарат, с целью снижения закрутки потока перед шнековым 
колесом. Однако, вместо выполнения своей основной функции, направ- 
ляющий аппарат насосного агрегата, образовывал условия для появле- 
ния сложных вихревых структур. Вихревые структуры данного типа 
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были источниками пульсаций давления, являлись причиной появления 
больших уровней вибрации. 

Для снижения пульсаций давления предложили изменить кон- 
струкцию лопаточной системы на входе в насос и таким образом оста- 
новить возникновение неоднородностей в потоке, и таким образом сни- 
зить уровень пульсаций давления. 

«Модернизированная» лопаточная система представляла собой 
«экран», выполненный из участка перфорированной трубы с фланцем, 
в который помещена лопаточная решетка с количеством лопаток боль- 
шим, чем в направляющем аппарате первоначальной конструкции. На 
рис. 2 представлен модернизированный участок трубопровода. 


Рис. 2. Участок трубопровода с перфорацией и направляющим аппаратом 


Как можно судить по результатам экспериментальных исследова- 
ний, вопреки применению данной конструкции, удалось снизить общую 
неоднородность потока на входе в насос и обеспечить незначительное 
снижение величин пульсаций давления и вибрации. 

В целом, выполненный обзор показывает, что не существует еди- 
ной методики снижения уровней неоднородности потока и способов 
предотвращения возникновения вихревых структур в перекачиваемой 
насосом жидкости. 

Однако, при детальном изучении конструкции насоса, а также 
структуры течения рабочей жидкости, была выдвинута идея по модер- 
низации насоса, позволяющая уменьшить интенсивность вихревого 
течения и тем самым снизить его влияние на основной поток. 

Реализация выдвинутой идеи, доработка а также результаты 
данной модернизации будут представлены в следующей работе. 

Выводы. 

1. Проанализированны причины возникновения вибрации в 
насосах. 

2. Перечислены последствия негативного влияния обратных токов 
на работу насоса. 

3. Проанализированны рассмотренные методы уменьшения нега- 
тивного влияния обратных токов на вибрационные характеристики цен- 
тробежного насоса. 
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[4] 
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REVIEW AND ANALYSIS OF THE EXISTING METHODS 
OF FIGHTING VORTEX FORMATION IN CENTRIFUGAL PUMPS 


The sources of hydrodynamic vibrations in centrifugal pumps are analyzed. 
Of particular interest in this work is the occurrence of reverse currents, eddy 
formation and uneven flow around the elements of the flow part. A picture of 
the flow of liquid in a centrifugal pump in the presence of reverse currents at 
the inlet is presented. This paper presents existing methods for reducing the 
negative impact of reverse currents on the vibration characteristics of a centrif- 
ugal pump. 


Keywords: hydrodynamic vibration, pressure ripple, centrifugal pump, re- 
verse currents, increasing the efficiency of the pump. 
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УДК 624.408.64 
В.В. Калмыков, Е.А. Медведева, Р.А. Конозобко 


СПОСОБЫ ОЦЕНКИ КРАЕВОГО УГЛА 
МЕТОДОМ ЛЕЖАЩЕЙ КАПЛИ 


Статья посвящена изучению способов измерения краевого угла 
смачивания одним из наиболее распространенных и доступных ста- 
тических методов – методом лежащей капли. Все из рассмотренных 
способов позволяют получить достаточно точные результаты из- 
мерений в зависимости от точности фиксации контура капли или 
пятна контакта со смачиваемой поверхностью. Исследуемые капли 
должны быть определенного размера и равного объема, чего можно 
добиться при помощи специальных дозаторов фирм «Ленпипет» с 
ценой деления 2-20 мкл и «Колор» с ценой деления 0,5-10 мкл. При 
этом нет необходимости применять специальное дорогостоящее 
оборудование, что играет немаловажную роль при планировании и 
проведении экспериментов в имеющихся лабораторных условиях. 


Ключевые слова: краевой угол смачивания, метод лежащей капли, 
дозаторы, смачиваемость. 


Многие химические и химико-технологические процессы и явле- 
ния связаны с растеканием жидкостей по твердым или жидким подлож- 
кам. Смачиванием называется совокупность явлений на границе сопри- 
косновения трёх фаз: твёрдое тело, жидкость и газ. Интенсивность сма- 
чивания характеризуется величиной краевого угла 0. В настоящее 
время большое количество работ посвящается вопросам смачивания, 
т. к. результаты теоретического и практического исследования статики 
и динамики капли играют важную роль при решении научных и техно- 
логических задач в ряде промышленных областей. Например, в тек- 
стильной промышленности при изучении смачиваемости одежды, при 
типографских работах, в ряде отраслей машиностроения при изучении 
смазывающей способности и т. д. [1-3]. 

Существует множество методов различной сложности по определе- 
нию статического краевого угла смачивания в лабораторных условиях — 
метод лежащей капли, метод диаграммы отраженного света, гониомет- 
рический метод и др. Данная работа посвящена анализу способов 
оценки краевого угла методомлежащей (неподвижной) капли, т. к. он 
является наиболее распространенным и относительно надежным из ста- 
тических методов и не во всех случаях требует использования сложного, 
дорогостоящего оборудования. 

При проведении экспериментов размер капли не меняется в тече- 
ние всего измерения, однако угол контакта изменяется под воздей- 
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ствием внешних факторов. Влияние седиментации, испарения и различ- 
ных химических или физических взаимодействий приводит к самопро- 
извольному изменению краевого угла со временем. 

Капля жидкости при помощи дозатора помещается на твердую по- 
верхность (рис. 1). Важным параметром является размер исследуемой 
капли. В случае очень малых капелек будет велико влияние поверхност- 
ного натяжения жидкости, а в случае больших капель преобладают силы 
гравитации, вследствие чего капля быстро растекается по поверхности. 
Наиболееоптимальный диаметр капли составляет 2-5 мм; при этих зна- 
чениях краевой угол не будет зависеть от диаметра [4,5]. 


Рис. 1. Дозаторы 


Соотношение сил межфазного и поверхностного натяжения в 
точке контакта трех фаз описывается на основании уравнения Юнга 


(рис. 2). 


Рис. 2. Определение краевого угла методом «лежащей капли»: 
1 — воздух; 2 — жидкость; 3 — твердое тело 


Предполагая, что капля жидкости находится в термодинамиче- 


ском равновесии, величина краевого угла смачивания: 
бу 3-0 
13702,3 
cos 0 = ————, (1) 
918 


где 0 — краевой угол смачивания, сүз - поверхностное натяжение на 
границе твердое тело-воздух; O54 — поверхностное натяжение на rpa- 
нице твердое тело-жидкость, O] — поверхностное натяжение на границе 


жидкость-воздух. 
Наиболее распространены следующие способы оценки формы ле- 
жащей капли: 
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1. Способ Юнга-Лапласа 

Способ основан на решении дифференциальных уравнений Юнга- 
Лапласа. При построении контура капли учитывается, что на форму 
капли оказывают влияние не только межфазные взаимодействия, но и 
собственный вес жидкости. Эта модель предполагает использование 
симметричности формы капли, поэтому её нельзя использоваться для 
исследования динамических краевых углов. Значительным недостатком 
способа является получение относительно точных результатов. Геомет- 
рические размеры капли определяются путем обмера ее фотографиче- 
ского изображения при помощи микроскопа или более простой реали- 
зации — фиксирования профиля капли при помощи цифрового фотоап- 
парата. Угол измеряется между касательной, проведенной к исследуе- 
мой жидкости и смачиваемой поверхностью твердого тела. Касатель- 
ную проводят через точку соприкосновения трех фаз: твердой фазы, 
жидкости, газа или воздуха. 

(412) =- 
cos0 = -a (2) 
(4/ 2) +h 
где 4 — диаметр капли, h — высота капли. 

Недостаток способа заключается в том, что процесс обмера доста- 
точно трудоемок, а его результаты содержат погрешность, связанную с 
индивидуальными особенностями наблюдателя [6, 7]. Исключить такие 
погрешности позволяют специальные приборы, например, DSAHT, 
DSA100, или портативный MobileDrop. 

2. Полуугловой способ 

Подходит для измерения капель малого размера. Каплю исследуе- 
мой жидкости можно представить как часть сферы (рис.3). Определив 
по полученному изображению высоту капли и её радиус, можно при- 


ближенно рассчитать угол контакта 6. 


Рис. 3. Измерение краевого угла полуугловым методом 


бане бе. (3) 


F 
где h — высота капли, r — половина ширины базовой линии. 
Тогда в случае простой геометрии угол краевой угол смачивания 
можно определить следующим образом: 
Ө = 20. 
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3. Круговой способ 

Способ также подходит для исследования капель малых размеров 
и является одним из наиболее распространенных и основывается на зна- 
нии простейших законов геометрии. 

По профилю полученного изображения капли строится окруж- 
ность радиуса г с центром в т. О (рис. 4). Базовая линия обозначается 
АВ, от которой будет измеряться краевой угол. На кривой указываются 
характерные точки А, В, С и точка пересечения базовой линии и верти- 
калью — т. D. 


Е 


Рис. 4. Определение краевого угла смачивания круговым способом 


В треугольнике АОС: 
ОА = OC = г, ХАСО = (САО =ß +y. АЕ- касательная к профилю 
капли. 


ZEAO =а+В+у=90°. (4) 
Из треугольника АРС: 
И САБ = /АСр= В+ В +уү= 90°. (5) 


Приравнивая выражения (4) и (5) следует: а = В. Тогда краевой 

угол смачивания будет определяться следующим образом: 
0 = агссоѕ(20). 

4. Тангенциальный способ 

Полученная форма капли представляется как часть контура пред- 
полагаемой окружности (рис. 5). Затем определяется центр предполага- 
емой окружности (т. М), а угол контакта представляет собой угол между 
касательной и построенной окружностью. 
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Рис. 5. Пояснение к определению краевого угла смачивания 
тангенциальным способом 


На предполагаемой окружности отмечаются три точки L1, L2 и L3 
для определения левого угласмачивания и аналогично Rl, R2, и R3 с 
противоположной стороны. Угол 0 определяется между касательной к 


профилю капли т и базовой линией. 

Рассмотренные способы обладают общим свойством, подразуме- 
вающем, что точность оценки краевого угла зависит от точности фикса- 
ции контура капли или пятна контакта со смачиваемой поверхностью. 
Прежде, такие способы, считались неточными и отдавалось предпочте- 
ние применению дорогостоящего оборудования. Сегодня, при исполь- 
зовании специальной оснастки, позволяющей формировать капли с точ- 
ностью до 0,5 мкл и четкого изображения контура капли, возможно оце- 
нивать краевой угол с точностью до минут. Следовательно, рассматри- 
ваемые способы методом оценки краевого угла могут быть использо- 
ваны для точных измерений. 
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METHODS FOR ASSESSING THE CONTACT ANGLE 
BY THE METHOD OF LYING DROP 


The article is devoted to the study of methods for measuring the contact an- 
gle by one of the most common and available of static methods - the method of 
а lying drop. All of the methods considered allow obtaining fairly accurate re- 
sults of measurement depending on the accuracy of fixing the contour of the 
drop or the contact spot with the wetted surface. The studied drops should be of 
a certain size and equal volume, which can be achieved with the help of special 
dispensers from «Lenpipet» firm with a division value 2-20 uL and «Color» 


firm with a division value 0,5-10 uL. In this case, there is no need to use special 


expensive equipment, which plays an important role in planning and conducting 
experiments in the existing laboratory conditions. 


Key words: contact angle, method of lying drop, dispensers, wettability. 
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УДК 621.65.03 
А.Н. Маклачков, В.В. Коротков 


СРАВНЕНИЕ ОЖИДАЕМЫХ И РАСЧЕТНЫХ 
ЭНЕРГЕТИЧЕСКИХ ХАРАКТЕРИСТИК СКВАЖИННОГО 
ПОГРУЖНОГО НАСОСА 


Проводится сравнительный анализ ожидаемых и расчетных 
энергетических характеристик скважинного погружного насоса. 
Расчетные характеристики получены при помощи автоматизиро- 
ванных средств, результаты представлены графически. Проведена 
оценка возможности применения методики расчета ожидаемых 
энергетических характеристик для многоступенчатого центробеж- 
ного насоса. 


Ключевые слова: ожидаемые энергетические характеристики 
насоса, максимально достижимое КПД насоса. 


Введение. 

Сложность рабочего процесса лопастного насоса затрудняет рас- 
чет его характеристики. Между тем уже при проектировании часто 
необходимо иметь эту характеристику, чтобы установить эксплуатаци- 
онные свойства насоса, оценить его эффективность и конкурентоспо- 
собность с имеющимися аналогами. 

В работе [1] приведена оригинальная методика, по которой проек- 
тант может самостоятельно рассчитать ожидаемые характеристики цен- 
тробежного насоса, исходя из его номинальных параметров и частоты 
вращения вала насоса. 

Целью данной работы является оценка возможности применения 
методики [1] на примере расчета ожидаемых характеристик многосту- 
пенчатого скважинного погружного насоса, основные параметры кото- 
рого приведены в таблице 1. 


Таблица 1 
Параметры проектируемого насоса 
Наименование параметра, размерность Значение 
Наружный диаметр корпуса насоса, не более, мм 94 
Подача номинальная, м/ч 4,0 
Напор номинальный, M 200 
Номинальная частота вращения, об/мин 2850 


По предварительной оценке данный насос имеет 50 ступеней. Ме- 
тодика определения количества ступеней скважинного погружного 
насоса подробно описана в работе [2]. 

Достоверность ожидаемых энергетических характеристик проекти- 
руемого насоса, будем оценивать, сравнивая их с характеристиками, по- 
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лученными по программе «Автоматизированное проектирование много- 
ступенчатого центробежного насоса» [3], которые условимся называть 
расчетными характеристиками. Программа [3] используется в научных и 
учебных целях в МГТУ имени Н.Э. Баумана (Калужский филиал) для 
проектирования многоступенчатых центробежных насосов с коэффици- 
ентом быстроходности n, = 40...300. Для проектирования проточных по- 


лостей используются методы расчета, основанные на одно- или двумер- 
ных моделях течения, которые создавались на протяжении десятков лет 
в МВТУ им. Н.Э. Баумана, МЭИ, ВНИИГидромаше, Ленинградском по- 
литехническом институте, ВНИИАЗН, МАИ и др. организациях. 
Определение максимально достижимого КПД насоса npa- 


Номинальные параметры насосов — подача и напор в номинальной 
точке работы насоса определяются в процессе проектирования без боль- 
ших ошибок. Хуже с оценкой коэффициента полезного действия (КПД). 
Поэтому широко используются статистические методы оценки дости- 
жимого КПД в зависимости от типа насоса, определяемого коэффици- 
ентом быстроходности п, . 


Рассмотрим наиболее известные способы определения максималь- 
ного КПД. 

Первый способ является наиболее популярным и подробно описан 
в работе [4]. Он используется инженерами и студентами технических 
ВУЗов для расчета КПД центробежных насосов. Значение максималь- 
ного КПД представляет собой произведение гидравлического, объем- 
ного и механического КПД: 


TN max 1,1,1, › (1) 
где n, – гидравлический КПД, n, - объемный КПД, Nye — механиче- 
ский КПД. 

Для определения гидравлического КПД: 
42 
е ш | (2) 


(ер, - 0.172) 


где О — приведенный диаметр горловины колеса, мм: 


О = 4000? 2 | (3) 
п 
где О – подача, м/с; n — частота вращения, об/мин. 
Для объемного КПД: 
1 
ов = 1+0, 685, C^ 9 , (4) 
где коэффициент быстроходности 
n, =3,65 "No (5) 


H^" 2, 
где Н – напор, м. 
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Значение механического КПД, не совершая большой ошибки, при- 
мем постоянным для всех значений n, = Па = 0,96 для прогнозных Ie- 


лей зто вполне допустимо. 

Второй способ заключается в определении максимально достижи- 
мого КПД с использованием режимных параметров работы насоса. В 
работе [5] приведены обобщенные данные по 132 высокооборотным 
насосам п, =16...290 и предлагается определять максимальный КПД 


как произведение двух сомножителей, один зависит от их, второй зави- 
сит от режимных факторов: 


ах - 1,1 ? (6) 
где 
1 

Ns = 1 , (7) 

+0,91+0,00054(n, —3) 

n — 
0,2 
пр -1-exp 2062. : (8) 
n 


Третий способ представляет собой, предложенное автором [1], ма- 
тематическое выражение зависимости КПД от п, и подачи 


Nmax = f (s> Ом), приведенной в графическом виде в [6]: 


Nmax = (0,042 1 О 


ном 


+0,56) (0,011 (1а n,) –0,23(1ап,) +1,5 ти, – 2,12). 
(9) 


Для оценки ожидаемого максимального КПД проектируемого 
насоса используем среднее значение, вычисленное по трем вышеопи- 
санным методикам. Результаты расчета представлены в таблице 2. 


Таблица 2 
Определение максимально достижимого КПД 
Коэффициент Максимально достижимый КПД Nmax 
быстроходности 
п, 1 способ | 2 способ | 3 способ | Среднее значение 
123 0,683 0,711 0,617 0,67 


Для построения напорной и мощностной характеристик нужно 
знать не только пъ» НО и зависимость 1 = f (О). Как известно, это па- 


рабола, однако ее форма также зависит от n, и О 


HOM * 


B работе [6] приве- 
день: графики парабол относительного значения т/т,„„, для нескольких 


значений п, (рисунок 1). 
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700046, Yo от номинальной 
Рис. 1. Кривые КПД для насосов с различными значениями п,. 


Таким образом, зная среднее значение максимально достижимого 
КИД и подобрав форму кривой по n, , была построена ожидаемая КПД xa- 
рактеристика проектируемого скважинного погружного насоса (рисунок 2). 


0,8 
0,7 


0,6 
0,5 


= 0,4 
А Расчетная хар-ка 


0,3 
9 Ожидаемая xap-ka 


О, «3/4 
Рис. 2. Расчетная и ожидаемая КПД характеристики 


Как видно из рисунка 2 ожидаемая КПД характеристика проекти- 
руемого насоса, построенная на основе описанных выше соображений, 
хорошо согласуется с расчетной характеристикой, полученной по про- 
грамме [3]. На номинальном режиме ожидаемое Npa Ha 3 % ниже 3Ha- 
чения, рассчитанного по программе [3]. 

Расчет ожидаемой мощностной характеристики насоса. Характе- 
ристика потребляемой насосом мощности № = f(Q) прогнозируется 


проще и достовернее. Для большинства насосов в широком диапазоне 
n., она представляет собой прямую линию, описываемую уравнением: 


МЕКО +К,, (10) 
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где К, = N, — потребляемая насосом мощность при О —0, а К, – опре- 


деляет наклон мощностной характеристики в рабочем диапазоне расхо- 
дов. Для построения зависимости достаточно определить две точки: №, 


при Q,, и №, при О=0. Потребляемая насосом мощность на номи- 


нальном режиме определяется по номинальным параметрам насоса: 


Н 
М, — Мом — ре ы С өм 5 (11) 
нах 


где р - плотность рабочей жидкости, d ; = - ускорение свободного 


падения VA А 


Для построения прогнозируемой зависимости N = f (О) необходима 
вторая точка - значение мощности насоса №, при О = 0. Она определяется 
исходя из накопленных статистических данных для насосов, близких по П 


В работе [6] приведены графики, из которых можно определить относитель- 


ную потребляемую насосом мощность a при Q=0 для n, = 47...650. 


Полученная зависимость хорошо аппроксимируется уравнением: 


N 
? = 0,0000097, + 0,00002л, --0,254 . (12) 


HOM 


Таким образом, определив N ou и No, зависимость N= f (О) 


определена во всем диапазоне подач (рисунок 3). 


= 
= 2 А Расчетная хар-ка 


Z 15 • Ожидаемая хар-ка 


о ~ 


0 2 4 6 
О, «3/4 
Рис. 3. Расчетная и ожидаемая мощностные характеристики 


Как видно на рисунке 3, ожидаемая мощностная характеристика про- 
ектируемого насоса хорошо согласуется с расчетной характеристикой, по- 
лученной по программе [3]. На номинальном режиме ожидаемая Nou 


на 2 % превышает расчетное значение. При малых подачах ожидаемые зна- 
чения мощности несколько ниже рассчитанных по программе [3]. 


Рубрика «Машиностроение» (ISSN 2413-6220) 62 


Определение напорной характеристики насоса в рабочем диапа- 
зоне подач. Зная зависимость потребляемой насосом мощности 
N = f (О) и КПД насоса в рабочем диапазоне подач, можно построить 


ожидаемую напорную характеристику (рисунок 4), по формуле: 


Ма 
н(0)---, 
(0) pgQ (13) 


n- f (0). 


А Расчетная хар-ка 


9 Ожидаемая хар-ка 


0 ud T T 1 
0 2 4 6 


а, м3/ч 


Рис. 4. Расчетная и ожидаемая напорные характеристики 


Как видно на рисунке 4, ожидаемая напорная характеристика про- 
ектируемого насоса в рабочем диапазоне подач хорошо согласуется с 
расчетной характеристикой. Исключение составляют малые подачи, на 
которых ожидаемый напор превышает расчетные значения. Это объяс- 
няется наличием обратных токов на входе в центробежное колесо при 
работе насоса с малыми подачами. Для упрощения расчета в данной ра- 
боте обратные токи не учитывались. Следует отметить, что на номи- 
нальном режиме значения ожидаемого и рассчитанного по программе 
[3] напора имеют полное совпадение. 

Выводы. 

1. Методика [1] расчета ожидаемых энергетических характеристик 
обеспечивает хорошую точность для многоступенчатого скважинного 
насоса. На номинальном режиме значения ожидаемого КПД и мощно- 
сти отличаются от рассчитанных по программе [3] на 3 % 2 Фо соответ- 
ственно. Значения ожидаемого и расчетного напора на номинальном ре- 
жиме полностью совпали. 

2. Для получения более точных ожидаемых значений мощности и 
напора проектируемого насоса на малых подачах, необходимо в расчете 
учитывать обратные токи. 

3. Методика [1] пригодна для получения ожидаемых энергетиче- 
ских характеристик проектируемых многоступенчатых центробежных 
насосов в рабочем диапазоне подач. 
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COMPARISON OF EXPECTED AND CALCULATED 
ENERGY CHARACTERISTICS 
OF A DOWNHOLE SUBMERSIBLE PUMP 


A comparative analysis of the expected and calculated energy characteris- 
tics of downhole submersible pump is performed. The calculated characteristics 
are obtained using automated tools, and the results are presented graphically. 
The expected energy characteristics were evaluated. The possibility of applying 
the method of calculating the expected energy characteristics for a multistage 
centrifugal pump is evaluated. 
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УДК 621.65.03 
А.А. Анкудинов, Н.С. Цыганов 


ЭНЕРГЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ 
ВИХРЕВОЙ ГИДРОМАШИНЫ С ЦИЛИНДРИЧЕСКИМ 
И КОНИЧЕСКИМ РОТОРОМ 


В статье проводятся результаты исследования новой вихревой 
гидромашины состоящей из ротора, в виде цилиндрического и кони- 
ческого барабана с глухими отверстиями и расположенного на его 
периферии винтового канала. Получены экспериментальные знерге- 
тические характеристики, определено влияние основных геометри- 
ческих параметров на характеристики, выбраны оптимальные зна- 
чения. Зависимости напора от подачи устойчивы во всем диапазоне, 
работа гидромашины не сопровождается низкочастотными пульса- 
циями давления. Гидромашина может найти широкое применение в 
различных областях техники. 


Ключевые слова: вихревая гидромашина, цилиндрический бара- 
бан, конический барабан, винтовая решетка, экспериментальные 
напорные характеристики, вихревое течение. 


Введение. Вихревые насосы и компрессоры с рабочим процессом, 
основанным на взаимодействии жидкости в рабочем колесе и боковом 
кольцевом канале, широко применяются в технике [1]. Работа вихревого 
насоса (рисунок 1) основана на передаче энергии лопатками рабочего 
колеса потоку жидкости в канале. Такие вихревые машины имеют вы- 
сокий коэффициент напора в их проточной части происходит интенсив- 
ный обмен количества движения жидкости, причем давление в канале 
возрастает линейно по мере движения потока. Вихревые гидравличе- 
ские машины, содержат рабочее колесо 1 (рисунок 1) с радиальными ло- 
патками и боковой канал 2. Имеются патрубки нагнетания 3 и всасыва- 
ния 5, которые разделены перемычкой 4. В боковом канале образуется 
вихревое течение, повышающее напор линейно от входа к выходу [2]. 


Рис.1. Схема вихревого насоса: 
1 — рабочее колесо; 2 — боковой канал; 3 — патрубок нагнетания; 4 — перемычка; 
5 — патрубок всасывания 
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Вихревъте насосы обычно применяются для создания достаточно 
большого напора при малой подаче. Насосы такого типа обладают свой- 
ствами самовсасывания, имеют крутопадающую напорную характери- 
стику и поэтому малочувствительны к колебаниям давления в системе. 

В вихревой машине, имеющей ротор в виде цилиндрического или 
конического барабана с рабочим процессом, основанном на вихревом 
взаимодействии жидкости в каналах проточной части, происходит по- 
добный процесс с образованием вихрей в глухих отверстиях и их рас- 
пространением в неподвижные каналы. 

Разработана вихревая реверсивная гидромашина, имеющая ротор 
с глухими сферическими отверстиями и течением основного потока по 
неподвижным винтовым каналам (рисунок 2). Такие устройства иссле- 
довались в широком диапазоне подачи и геометрических параметров. 
Насос состоит из ротора в виде цилиндрического или конического бара- 
бана с глухими отверстиями на наружной поверхности и охватывающей 
его неподвижной винтовой решетки |3, 4]. При вращении ротора под 
действием центробежных сил и воздействия винтовой решетки на поток 
в глухих отверстиях ротора возникают вихри, имеющие жесткую струк- 
туру ираспространяющиеся в каналы решетки (рисунок). Вихри распро- 
страняются в каналы решетки и перемещают находящуюся там жид- 
кость от входа к выходу. Нетрадиционная схема гидромашины приво- 
дит к качественно новым результатам. Эксперимент показал, что напор 
такого устройства прямопропорционален осевой длине и массе жидко- 
сти, участвующей в осевом течении. 


A-A 
разбертка 


Рис. 2. Схема вихревой гидромашины с цилиндрическим ротором 
и неподвижными винтовыми каналами: 
1 – цилиндрический барабан; 2 — винтовые каналы 
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Рассматривая схему работы (рисунок 3) примем допущения: тече- 
ние жидкости установившееся; отсутствует момент от трения жидкости 
о стенки неподвижного канала и поверхность барабана; давление возрас- 
тает по длине винтового канала линейно; картина течения в области всех 
глухих отверстий, включая входные и выходные участки, одинакова. 


Рис. 3. Схема течения жидкости в глухих отверстиях барабана 
и винтовом канале 


Выделим сектор, а на секторе угол Ô , соответствующий диаметру 
глухого отверстия в плоскости перпендикулярной оси барабана. При 
вращении ротора под воздействием центробежных сил и лопастей вин- 
тового канала жидкость в глухих отверстиях перемещается с окружной 
составляющей скорости у. и нормальной у,. Обозначим: Е - площадь 


сечения глухого отверстия; г — количество глухих отверстий; v, — 
окружная составляющая скорости жидкости на площадке Е; v,— нор- 


мальная составляющая скорости жидкости на площадке Е. 

Влияние геометрических параметров на энергетические характе- 
ристики определялось экспериментальным путем. [5, 6]. Исследовались 
проточные части с наружными диаметрами барабанов D, =126; 140; 


189; 470 мм и различным числом глухих отверстий. На наружных по- 
верхностях барабанов выполнялись отверстия различного диаметра, 
формы, глубины. Винтовые решетки имели шаг одного витка 
5 = 64; 112; 160;192 мм. Высота лопастей составляла h, =9,5;12;18,5; 

24; 31 мм. Число каналов винтовой решетки выбиралось <: =3;4;5; 6. 
Толщина лопастей неподвижной винтовой решетки составляла 
Ь=2...3мм. Число глухих отверстий в барабане достигало z,- 270. 
Всего было испытано более 30 вариантов вихревой машины. Испытания 
проводились на рабочей жидкости - воде с температурой в диапазоне 
t = 20...80 °С. Частота вращения барабана n = 800...2000 об/мин. Диапа- 
зон расходов при определении характеристик составлял 
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О =0...300 м/ч. Приведены напорные характеристики вихревой Ma- 
шины H = f(Q) отличающейся только отношением глубины глухого 


отверстия к его диаметру i^ (рисунок 4). Из приведенных характери- 
0 


стик следует, что выполнение на барабане глухих отверстий приводит к 
значительному увеличению напора. 


, n - 1500-22. 


| —— 
3 


м 
в == 
ч 


Рис. 4. Напорные характеристики с различным соотношением глубины 
к диаметру глухих отверстий в барабане: 
1— d) =15 мм, № =4 мм; 2- d, < 10 мм, № =8 мм; 3 — 4у=15 мм, № 212 мм; 
4— dg =16 MM, № =12 мм; 5 – dg =18 MM, № =15 мм 


Испытания на различном числе оборотов моделей показали, что 
для исследованного диапазона скоростей выполняются законы динами- 
ческого подобия и справедливы соотношения для модельного пересчета 
энергетических характеристик по известным формулам. В приведенном 
диапазоне геометрических параметров, при выполнении их оптималь- 
ных значений, получена эмпирическая формула расчета напорной ха- 
рактеристики: 

2 2 2 
9 Di | 0,1+0,93sin a, сова, № |-(1-1,674+ 0,674), (1) 
8 D, (D; -D, ) 


где а = 0/0,,„; О - текущее значение подачи; О, – максимальная IO- 


Н = 


дача, при которой напор равен нулю. Максимальный расход рассчиты- 
вается как половина расхода через винтовую решетку при условии без- 
ударного входа О = 0,5F u/cosa, . [7]. 

Вихревую гидромашину имеющую ротор в виде конического ба- 
рабана можно испытывать в качестве компрессора [3, 4]. На рисунках 5 


тах 
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и 6 представлен вариант вихревого компрессора, когда область всасы- 
вания находится на большем диаметре ротора, хотя возможен и вариант 
машины, когда патрубок всасывания находится на меньшем диаметре. 


Рис. 6. Схема неподвижной винтовой решетки 


При вращении ротора под действием центробежных сил и воздей- 
ствия винтовой решетки на поток в глухих отверстиях ротора возникают 
вихри. Вихри распространяются в каналы решетки и перемещают нахо- 
дящуюся там среду от входа к выходу. Основной поток перемещается 
по неподвижным каналам. Упрощенную схему работы компрессора 
можно представить как работу объемного пластинчатого компрессора, 
где роль пластин выполняют вихри, возникающие в глухих отверстиях, 
а роль отвода винтовые каналы. Приведены напорные характеристики 


вихревой машины Н = f (Q, ) (рисунок 7). 
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Рис. 7. Напорные характеристики вихревой машины 


Напорные характеристики определялись при постоянной частоте 
вращения, которая поддерживалась постоянной. Работа компрессора 
устойчива во всем диапазоне изменения подачи. 

Выводы 

1.При снижении давления на входе происходит плавное снижение 
напора, низкочастотные пульсации давления отсутствуют. 

2.Работа машины устойчива во всем диапазоне подачи. 

3.Полученные данные открывают перспективу использования раз- 
работанного устройства в качестве насоса, компрессора, нагрузочного 
устройства, так как ротор представляет собой барабан с большим запа- 
сом прочности. 

4.На базе представленной вихревой гидромашины возможно со- 
здание герметичных насосов для химической, нефтеперерабатывающей, 
и газовой промышленности, авиации и судостроения. 
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ENERGY CHARACTERISTICS OF VORTEX HYDRAULIC 
MACHINE WITH CYLINDRICAL AND CONIC ROTOR 


The article presents the results of a study of a new vortex hydraulic machine 
consisting of a rotor in the form of a cylindrical and conical drum with blind 
holes and a helical channel located on its periphery. Experimental energy char- 
acteristics are obtained, the influence of the main geometric parameters on the 
characteristics is determined, the optimal values are selected. Dependences of 
the pressure on the flow are stable over the entire range, the operation of the 
hydraulic machine is not accompanied by low-frequency pressure pulsations. 
The hydraulic machine can be widely used in various fields of technology. 


Keywords: vortex hydraulic machine, cylindrical drum, conical drum, heli- 
cal lattice, experimental pressure characteristics, vortex flow. 
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УДК 004.921 
Ю.Ю. Болденков, Ю.С. Белов 


МЕТОД ЛОКАЛИЗАЦИИ ПОЛЬЗОВАТЕЛЯ 
В ПОМЕЩЕНИИ НА ОСНОВЕ МАРКЕРОВ 
ДЛЯ MICROSOFT HOLOLENS 


Дополненная реальность (Augmented Reality, AR) как популярный 
метод взаимодействия человека с компьютером (Нитап-Сотршег 
Interactions, НСІ), касается взаимодействия между пользователем и 
компьютером, особенно в области объединения виртуальных и ре- 
альных элементов. Мобильные устройства дополненной реальности 
(AR) для внутренних помещений, такие как Microsoft HoloLens, обла- 
дают потенциалом для визуализации данных о модели здания на ме- 
сте. Хотя HoloLens обладает достаточной способностью отсле- 
живания в режиме реального времени, чтобы обеспечить простран- 
ственно-правильную и стабильную визуализацию виртуального кон- 
тента относительно его окружения, размещение виртуальных объ- 
ектов обычно должно выполияться пользователем. Рассматрива- 
ется метод локализации на основе маркеров для HoloLens, который 
достаточен для наложения внутренней среды на виртуальные дан- 
ные модели в масштабе комнаты с пространственной точностью в 
несколько сантиметров. 


Ключевые слова: дополненная реальность, HoloLens, локализация, 
визуализация. 


Введение. Как подсфера смешанной реальности (Mixed Reality, 
МВ), АК добавляет виртуальные элементы к воспринимаемой реально- 
сти и позволяет взаимодействовать в реальной среде, в то время как 
пользователь получает дополнительную визуальную компьютерную 
или смоделированную информацию для поддержки поставленной за- 
дачи. В прошлом АК-средьг применялись главным образом в научной 
визуализации и игровых развлечениях. В последние годы они были изу- 
чаются для применения в образовательных и визуализационных целях в 
АЕС (Architecture, Engineering, and Construction) и для взаимодействия 
при совместном проектировании [1]. 

С помощью HoloLens (версии Microsoft 2018а) корпорация 
Microsoft выпустила мобильную АК-платформу с головной консолью, 
которая способна обеспечить высококачественное отслеживание из- 
нутри и пространственно стабильную визуализацию виртуального кон- 
тента (так называемые «голограммы») в помещениях [12]. Исходя из 
технической оценки этого АВ-устройства можно сделать вывод, что оно 
способно отслеживать объекты с точностью до нескольких сантиметров 
и с точностью до миллиметра в случае медленных движений [1, 2]. 

Описание и характеристики HoloLens. HoloLens это дисплей- 
ный блок, монтируемый на голове, с расположенными перед глазами 
тонированными линзами с волнообразной призматической структурой. 
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В верхней части расположены 2 пары кнопок: для управления яркостью 
и громкостью. Динамики расположены у нижнего края устройства; они 
позволяют слышать, как звуки виртуальной реальности, так и звуки, ис- 
ходящие извне [3]. 

HoloLens автономны и не требуют подключения к ПК или смарт- 
фону, в отличие от большинства других устройств виртуальной, допол- 
ненной или смешанной реальности. 

В основе HoloLens лежит 64-разрядный 4-ядерный процессор Intel 
Atom x5-Z8100 с частотой 1,04 ГГц. В дополнение HoloLens имеет голо- 
графический процессор holographic processing unit (hpu), разработанный 
Microsoft специально для HoloLens. 

Встроенное хранилище данных имеет объем 64 Гб, однако пользо- 
вателю остаются доступны 54,09 Гб так как около 10 Гб занимает опе- 
рационная система. Объем оперативной памяти - 2 Гб. 

Приложения для HoloLens не могут использовать больше 900 Мб па- 
мяти. При превышении этого лимита работа приложения прерывается. 

Для корректной работы со смешанной реальностью HoloLens 
оснащены 4 камерами (по 2 с каждой стороны) для сканирования окру- 
жения и ориентации в пространстве, 4 микрофонами, гиростабилизато- 
ром, датчиком глубины, 2МР видеокамерой, сенсором окружающего 
освещения. 

Оптика HoloLens устроена очень сложно, что обусловлено необхо- 
димостью не просто выводить изображение на экран, но еще и пра- 
вильно совмещать его с объектами реального мира. Жидкокристалличе- 
ские проекторы с разрешением сторон 16:9, создают изображение, ко- 
торое затем проходит через визуализационную оптику, волновод, 
combiner (устройство, совмещающее проекцию и изображение реаль- 
ного мира) и дифракционные решетки. Линзы имеют 3 слоя - для си- 
него, зеленого и красного цветов. Каждый слой имеет свои дифракци- 
онные свойства [4]. 

В контексте большинства существующих приложений HoloLens 
отдельные виртуальные объекты размещаются в среде пользователя 
устройства HoloLens. Обычно это происходит либо активным «ручным 
размещением» подобных виртуальных объектов пользователями с по- 
мощью жестов, либо, альтернативно, автоматическим размещением. 
Автоматическое размещение в большинстве случаев ограничивается 
размещением виртуальных объектов на открытых пространствах на 
полу, на столах или на стенах таким образом, чтобы естественным 06- 
разом интегрировать их в среду пользователя [5]. 

При использовании устройства дополненной реальности для визу- 
ализации на месте компоненты модели здания, такие как, например, 
стены, двери или кабели внутри стен можно представить, как «голо- 
граммы». 

В этом случае размещение этих компонентов виртуальной ком- 
наты должно выполняться таким образом, чтобы гарантировать, что 
виртуальные компоненты здания точно перекрывают свои материаль- 
ные объекты в окружении пользователя [6]. 
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Решение задачи локализации пользователя в помещении при 
помощи маркеров. Для достижения пространственно корректного до- 
полнения сцены внутри помещения с помощью содержимого модели 
виртуального здания на HoloLens необходимо знать пространственную 
взаимосвязь между средой, в которой работает устройство, и его соот- 
ветствующей моделью [7]. 


Это сводится к определению позиции га. модели здания в си- 


стеме координат приложения HoloLens. Велична ТА обътчно обозначает 
позу объекта А в координатной сетке В и состоит из матрицы 4х4 в виде: 


B ‚В 
Ка t 
А 
0 1 
Здесь R? является ортогональной матрицей 3х3, которая описы- 
вает ориентацию соответствующего объекта А относительно системы 


координат В то время как вектор 1; описывает свою позицию. 


Для наложения реальных объектов с данными модели виртуаль- 
ного здания в виде голограмм необходимо определить положение 
устройства HoloLens относительно модели здания. Это равносильно 


определению положения ТАРР модели здания в коо динатной системе 
Model 


приложения HoloLens, что приводит к правильному выравниванию 
между виртуальной и физической геометрией здания. 

Система координат приложения HoloLens — это система коорди- 
нат, в которой описывается положение устройства HoloLens, определя- 
емое его системой отслеживания внутри помещения. Определение ко- 
ординатной рамки приложения получено из положения устройства в мо- 
мент запуска соответствующего приложения HoloLens. 


Арр ~ 
Текущая позиция Ту, „л Устройства HoloLens в этом этапе pa- 


боты приложения может запрашиваться приложениями HoloLens через 
HoloLens SDK. 
Кроме того, произвольный виртуальный объект («голограмма») 


А 
ТАРР 


может быть расположен путем установки его позы Тир ит 


в этом кадре 


приложения через SDK. 
Таким образом, построение данных модели может располагаться в 
кадре приложения в виде голограмм, если можно определить правиль- 


ную позицию Т?Р модели здания в кадре приложения. 
Model 


Самый простой способ добиться этого - поместить маркер B среду 
здания, которая должна быть дополнена данными модели. Положение 
этого маркера в системе координат соответствующей модели здания 
должно быть известно. Если позицция этого маркера также может быть 
определена в кадре приложения HoloLens, то в результате положение 
модели определяется как: 


ТАРР = ТАРР T Маг ker App Model ! 


Model Marker* Model T * Магкег! Marker * 
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Положение маркера в кадре приложения можно определить, 


наблюдая за ним с помощью камеры «Live Capture» HoloLens и, таким 


T HoloLens 


образом, определяя ero положение Тул. В локальной системе KOOP- 


динат этой камеры. Положение камеры в кадре приложения в момент 
захвата изображения маркера может быть запрошено через HoloLens 


Арр : 
SDK. Положение маркера Туу, в кадре приложения приводит к: 


ТАРР np App т HoloLens 
Marker ^. ^ HoloLens” Marker 


Схематический обзор зтого метода Hà основе маркеров для лока- 
лизации устройства HoloLens внутри моделей зданий, соответствующих 
окружающей его внутренней среде, показан на рисунке 1 [8-10]. 

Комната может быть дополнена соответствующими данными мо- 
дели с использованием Microsoft HoloLens путем определения положе- 


Арр ~ ~ 
НИЯ Тае модели виртуальнои комнаты внутри координатнои рамки 


приложения HoloLens. Положение n. Lens Самого устройства HoloLens 


относительно фрейма приложения может быть запрошено через 
HoloLens SDK [9, 11]. 


Комната 


Рис. 1. Схема модели 


Заключение. Надетые на голову АК-устройства могут обеспечи- 
вать средства захвата изображений, которые сами включают в себя ви- 
зуализацию виртуальных объектов, дополняющих сцену, как в случае с 
HoloLens и его камерой «Live Capture» [13]. Однако ни в коем случае 
нельзя быть уверенным в том, что видимое положение виртуальных 
объектов, изображенных на подобных изображениях, соответствует 
расширенной сцене, которую он испытывает при просмотре изображе- 
ния на дисплее установленного на голове АВ-устройства [14]. 

В случае Microsoft HoloLens видимое положение виртуальных 
объектов B дополненных изображениях, снятых камерой «Live Capture», 
может заметно отличаться от их видимого положения, которое наблю- 
дает пользователь, носящий устройство. Поэтому для адекватной и объ- 
ективной оценки точности позиционирования виртуальных объектов 
крайне важно всегда измерять видимое положение голограммы непо- 
средственно во время ее восприятия пользователем, а не с помощью до- 
полнительного механизм захвата, как у камеры HoloLens «Live Capture». 
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MARKER BASED USER LOCALIZATION METHOD 
FOR MICROSOFT HOLOLENS 


Augmented Reality (AR) as a popular technique in Human-Computer Inter- 
action (HCI), which deals with the interaction between the user and the com- 
puter, especially in the field of combining virtual and real-world elements. Mo- 
bile devices which support augmented reality (AR) for indoor environments, 
such as the Microsoft HoloLens, have the potential to visualize building model 
data in place. While the HoloLens has enough real-time tracking capability to 
provide spatially-correct and stable visualization of virtual content relative to 
its environment, the placement of virtual objects should be usually done by the 
user. Consider a simple marker-based localization method for the HoloLens, 
which is sufficient to impose the indoor environment on the virtual data of the 
model at the scale of a room with a spatial accuracy of several centimeters. 


Keywords: augmented reality, HoloLens, localization, visualization. 
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УДК 004.93.11 
В.Ю. Кириллов, А.А. Кузнецова 


ПОСТРОЕНИЕ МОДЕЛЕЙ ДЛЯ АНАЛИЗА КОМПОЗИЦИИ 
РЕТИНАЛЬНЫХ СНИМКОВ 


В работе проводится исследование способов анализа набора 
снимков ретинальной камеры. Дан обзор существующих плоскост- 
ных и пространственных способов композиции. Предложены два ва- 
рианта расширения технологии сверточных нейронных сетей на 
сферическую поверхность для альтернативных способов ее описания 
и рассмотрены правила обхода. Выявлены потенциальные трудно- 
сти обработки каждого случая. 


Ключевые слова: анализ ретинальных снимков, адаптивное ядро 
сверточной нейронной сети, квазиравномерная сетка на сфере. 


Для построения скрининговой системы, диагностирующей заболева- 
ния по набору ретинальных кадров, требуется инструмент выделения опре- 
деленных паттернов и закономерностей (области патологий, сосудистая 
сетка и т.п). Графическая форма представления и большая информацион- 
ная насыщенность каждого кадра, а также сложность и нелинейность рас- 
сматриваемых ключевых объектов стимулируют применять сверточные 
нейронные сети, архитектуры, предназначенные для распознавания обра- 
зов. Проблематичность прямого их использования обуславливается тем, 
что поверхность глазного дна не является плоскостью и при съемке возни- 
кают нелинейные искажения, значительно усиливающиеся при съемке пе- 
риферийных областей с отклонением от главной оптической оси. 

С этой проблемой не сталкиваются диагностические системы, рас- 
сматривающие только фронтальный кадр, в центре которого располо- 
жены макула и диск зрительного нерва. В данном случае достаточно 
применить классический подход к сверточным нейронным сетям, изоб- 
раженный на рисунке 1. 


„во Бов 


Рис. 1. Схема работы классической сверточной сети. 
Область соседства определяется по прямоугольной сетке 


Зффективность применяемой методологии обуславливается срав- 
нительной простотой машинных вычислений при обходе изучаемой ре- 
гулярной прямоугольной сетки. 

К сожалению, для диагностирования таких заболеваний как рети- 
нопатия недоношенных требуется полная информация о состоянии пе- 
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риферии, в которой наблюдаются области фиброзной ткани, аваскуляр- 
ные области и пролиферативный вал. Современные ретинальные ка- 
меры имеют угол обзора 120°-130°, что мало сопоставимо с полным те- 
лесным углом, поэтому необходимо оценивать пространственную ком- 
бинацию полученных снимков. При этом нелинейный характер искаже- 
ний не позволяет совместить снимки на плоскости и прямо транслиро- 
вать существующий матричный подход. Необходимость приближения 
пространственно-масштабных характеристик модели реальному глазу, 
приводит к потребности исследовать не плоскую комбинацию снимков 
(которая будет ухудшаться, например, при необходимости стыковки 
трех и более кадров, т.к. сосуды неодинаковы на них во взаимоотдале- 
нии и ширине), а наложение комбинации на некую криволинейную по- 
верхность: сферу [1], полусферу, анатомически приближенную форму. 
Таким образом, предметом исследования является построение модели 
глазного дна, получаемой из совокупности ретинальных снимков, кото- 
рая пригодна для применения методов машинного обучения, в частно- 
сти, сверточных сетей, сегментирующих сетей (U-Net) и т. п. Модель 
должна обладать регулярной структурой, для которой можно опреде- 
лить ядро свертки, правила обхода и субдискретизации (pooling). 

Существует ряд решений, экстраполирующих механизм свертки с 
плоскости на трехмерное пространство, в качестве ядра свертки в таком 
случае выступает кубическая матрица. Для рассматриваемой задачи 
этот подход является слишком прямолинейным, так как приводит к 
огромной размерности образца. Известны и способы улучшения каче- 
ства распознавания, например, мультипроекционный подход, позволя- 
ющий снизить влияние ориентации объекта на конечный результат [2]. 
Однако на основании имеющихся кадров почти невозможно судить о 
форме глазного дна, а решения, рассматривающие именно цвето-тоно- 
вые характеристики неплоской поверхности сравнительно редки [3, 4]. 

Необходимо помнить, что сверточные сети устойчивы лишь к аф- 
финным преобразованиям, в то время как искажения, возникающие при 
съемке в сложной оптической системе, не являются таковыми. Для уста- 
новления имеющихся зависимостей (в частности, формы сосудистой 
сетки) может потребоваться отход от классической формы ядра свертки, 
что усложнит описание процесса, но не обесценит решение вообще. Так, 
например, для случая единственного кадра в работе [5] показано приме- 
нение адаптивной сверточной сети с деформируемым ядром (которое, 
тем не менее, располагается в регулярной структуре). 

Для простого перенесения операции свертки с прямоугольной 
сетки на сферу необходимо задать на сфере некоторую равномерную 
сетку (критерий равномерности, связан с чувствительностью к пово- 
роту), однако решения данной проблемы не существует. Полигональ- 
ные сетки, которыми сфера приближается в компьютерной графике, не- 
регулярны за счет присутствия пятигранных элементов в некоторых 
точках (рис. 2). Можно описать дополнительные правила соседства для 
областей пятигранников, но это придаст дополнительную сложность 
правилам обхода, особенно учитывая, что данные правила не эквива- 
лентны для различных размеров ядра сверки и страйда. 
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Рис. 2. Проблема триангуляции сферы 


Другим способом задания регулярной системы точек на поверхно- 
сти сферы является редуцированная Гауссовская сетка, которая обла- 
дает свойством равенства расстояний между соседними точками. Такая 
сетка не описывается матрицей, число точек на параллелях, а соответ- 
ственно и количество соседей, уменьшается к полюсам, что делает ее 
неприменимой для задачи свертки — нельзя описать правила обхода и 
невозможно выделить одинаковое ядро. 

В качестве решения проблемы работы свертки на сфере предлагается 
отказ от использования точек на ее поверхности. В процессе обхода сферы 
ядро можно перемещать по плоскостям, касательным к центральной точке 
области, к которой будет применяться ядро. Рассмотрим плоскость АВСР 
на рисунке 3, здесь точка 7'— центр малого квадрата ABCD на касательной 
плоскости, О — центр сферы AJ, B1, C1, DI - проекции соответствующих 
точек на центр сферы. Тогда применяя свертку с размером ядра 2х 2 с цен- 
тром в точке T можно составить матрицу из значений проекций этих точек 
на плоскость (А ВС“ — для вершин квадрата). Таким образом, сохраня- 
ется работа с квадратной сеткой, удобной для машинной обработки, и оста- 
ется вариативность в правилах обхода поверхности. 


Рис. 3. Метод касательных плоскостей 
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Данный способ применим для любой выпуклой поверхности, что 
ценно в случае использования уточненной формы глаза. Другим его до- 
стоинством является то, что поверхность не обязана быть описана как 
набор точек, поскольку по итогам склейки ретинальных кадров, такое 
описание можно получить лишь со значительной долей приближения. 
Классические способы соединения изображений основаны на поиске де- 
скрипторных точек и сопоставлении их на соседних кадрах. Они дают 
плохой результат для снимков сетчатки, поскольку исходные изображе- 
ния неодинаково искажены и из-за особенностей предметной области 
дескрипторы разных точек похожи [6]. 

Сложность представляет совмещение снимков на неплоской по- 
верхности глаза. В работе [1] был предложен способ решения этой за- 
дачи вручную. В настоящей работе модель глазного дна представляет 
собой совокупность опорных плоскостей 7, T5, ..., Ту (соответствую- 
щих набору исходных ретинальных снимков) касательных к сфере в не- 
которых точках (рис. 4) и правило расчета значений на плоскости, каса- 
тельной к сфере в произвольной точке Т, как линейную комбинацию 


проекций плоскостей с весовыми коэффициентами, характеризующими 
степень близости к точке касания. 


N 


NN ра 


Рис. 4. Метод касательных плоскостей 


Так, область сферы аппроксимируемая плоскостью } (Т) ‚ описы- 
вается за счет комбинации проекций 1-й, 2-й и 4-й плоскостей 4 (Т, ‚Т,) 
в каждую ее точку с весовыми коэффициентами с, — зависящими от 
расстояния Ту до точки касания кадра (1). 

f (15) = € *d (19,7; - 9» *d (T9, T) - o, “а (11) (1) 

Коэффициент о; может представлять собой и функцию, отражаю- 


щую помимо близости соседства качество конкретного снимка. Реаль- 
ные наборы ретинальных снимков содержат области засветки, бликов, 
размытия, а также слабо различающиеся кадры. К сумме весовых коэф- 
фициентов для каждой точки предъявляется требование нормировки — 
равенства единице, что нетривиально для случая зависимости от рассто- 
яния и качества снимка. Пусть ) порядковый номер соседа начиная с 0, 
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в множестве, упорядоченном по угловому расстоянию, тогда компо- 
нента близости, зависящая от постоянной 0: 
=; -a(1- a) (2) 
А собственно козффициент зависит or m; — Beca качества кадра, 
задающего плоскость /, и компонента близости и нормируется относи- 
тельно всех попавших в выборку соседей: 


т же, 


O. 7 7 (3) 


ит 
тр те, 
Я 


Проблема определения качества снимка сложна еще и тем, что ар- 
тефакты, возникшие при съемке, могут быть выявлены именно при со- 
поставлении нескольких кадров, а этот процесс трудноосуществим 
стандартными способами из-за недостатка и малой вариативности де- 
скрипторных точек на фотографиях сетчатки [6]. 

Другой вариант расширения технологии свертки на сферическую 
поверхность — аппроксимация сферы проекциями на грани усеченного 
икосаэдра — описан в работах [7, 8]. Преимуществами этого подхода яв- 
ляются адаптивность к разрешению исходных снимков и возможность 
применения для анализа посредством плоских сверточных преобразова- 
ний (возможность использования имеющихся библиотек и трансфер- 
ного обучения). Подход основывается на использовании гномонических 
проекций и для рассматриваемой предметной области прямо не приме- 
ним — ретинальные снимки делаются с поверхности глаза и обладают 
сильной неоднородностью по цветопередаче и артефактам. 

Модель, предлагаемая в настоящей работе, позволяет однозначно 
определить гномонические проекции по совокупности снимков с воз- 
можной коррекцией по качеству снимков. 

Таким образом, для использования технологии свертки в диагно- 
стической системе, работающей с наборами ретинальных снимков, 
предложены вариант с касательными плоскостями, получаемыми есте- 
ственным образом при совмещении кадров, или работа с их проекцией 
на описанный многогранник с большим числом граней. Хотя незнание 
оптических параметров системы не дает делать предположения об ис- 
тинном виде искомой комбинации, оба предложенных способа характе- 
ризуются минимальным искажающим влиянием на внутренние связи 
внутри кадров, что говорит о релевантности анализа. 


СПИСОК ЛИТЕРАТУРЫ 


[1] Терещенко А.В. Компьютерный анализ сетчатки и ретинальных сосудов 
при ретинопатии недоношенных / А.В. Терещенко, Ю.А. Белый, 
М.С. Терещенкова, И.Г. Трифаненкова, А.А. Кузнецов, Ю.А. Юдина // 
Офтальмохирургия. – 2009. – № 5. С. 48—51. 

[2] ОС. Volumetric and Multi-View СММ for Object Classification on 3D Data 
/ C. Qi, Н. Su, M. Niefiner // IEEE Conference on Computer Vision and Pat- 
tern Recognition. — 2016. — P. 5648—5656. 


Электронный журнал: наука, техника и образование (ISSN 2413-6220) 83 


[3] Su Y.-C. Learning Spherical Convolution for Fast Features from 360? Im- 
agery / Y.-C. Su, K. Grauman // Clin. Perinatol. — 2018. — Vol. 40. — No. 2. — 
P. 271—296. 

[4] > Cohen T. Spherical CNNs / T. Cohen, M. Geiger, J. Kohler // Proceedings of 
the International Conference on Learning Representations. — 2018. — P. 15. 

[5]  Qiangguo J. DUNet: A deformable network for retinal vessel segmentation. / 
J. Qiangguo, M. Zhaopeng, Ph. Tuan // Journal of Class Files — 2015. Vol. 14. 
№ 8. — P. 1-12. 

[6] > Кузнецова А.А. Исследование способов автоматизации диагностирова- 
ния ретинопатии недоношенных / АА. Кузнецова, Е.В. Вершинин // 
Наукоемкие технологии. - 2019. — T. 3. – № 10. – C. 146-149. 

[7] Eder M. Convolutions on spherical images /М. Eder, J.-M. Frahm. // IEEE 
Conference on Computer Vision and Pattern Recognition Workshops. — 2019. 
— Р. 5. 

[8] Eder M. Tangent Images for Mitigating Spherical Distortion. / M. Eder, 
M. Shvets, J. Lim. — [Электронный ресурс]. - URL: https://arxiv.org/abs/ 
1912.09390 (дата обращения: 16.01.20). 


СВЕДЕНИЯ ОБ АВТОРАХ 


Кириллов Владимир Юрьевич (КиШоу V.Yu.) — канд. физ.-мат. 
наук, доцент кафедры «Системы обработки информации» КФ МГТУ 
им. Н.Э. Баумана; sapr91 @ mail.ru 


Кузнецова Анна Андореевна (Kuznetsova А.А.) - студент 
КФ МГТУ им. Н.Э. Баумана; 521annakuznetsova @ gmail.com 


MODEL BUILDING FOR FUNDUS IMAGES 
COMPOSITION ANALYSIS 


In this paper retinal camera image set analysis methods as well as existing 
approaches to combine sets on the plane and in three dimensions were studied. 
We introduce two ideas of extending convolutional neural network to spherical 
surface which can be described either as explicit points or as a set of tangent 
planes. Both cases are supplemented with surface bypass rules and potential 
edge difficulties. 
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Е.В. Вершинин, И.В. Лаковщиков, А.С. Никулин 


ПРИМЕНЕНИЕ РЕКУРРЕНТНЫХ НЕЙРОННЫХ СЕТЕЙ 
ДЛЯ ОПРЕДЕЛЕНИЯ ЭМОЦИОНАЛЬНОЙ ОКРАСКИ ТЕКСТА 


В работе рассматривается применение рекуррентных нейронных 
сетей для определения эмоциональной окраски текста. Представлен 
набор данных для обучения и описан способ их представления в циф- 
ровом виде. Рассмотрен процесс построения модели нейронной сети, 
с использованием архитектуры Г5ТМ. Описано обучение полученной 
модели и сделаны выводы по результатам работы на проверочных 
данных. 


Ключевые слова: рекуррентные нейронные сети, LSTM, змоцио- 
нальная окраска текста, плотное векторное представление, Python, 
Keras. 


Введение. На сегодняшний день задача по анализу змоциональной 
окраски текста имеет большую популярность. Это обусловлено тем, что 
количество информации в текстовом представлении постоянно растет. 
Маркетологи, лингвисты, государственные служащие и ряд других про- 
фессий постоянно сталкиваются с задачами по обработке текста. Авто- 
матизация данного процесса может привести к увеличению качества и 
скорости выполнения работы. 

Данная задача входит в раздел компьютерной лингвистики, 
ставящей своей целью представление естественных языков в виде ма- 
тематической модели. Благодаря развитию технологии искусствен- 
ного интеллекта решение поставленной проблемы не требует углуб- 
ленного изучения в сфере машинного обучения и искусственных 
нейронных сетей. Современные инструменты программирования o6- 
ладают уровнем абстракции, позволяющим быстро вникнуть в про- 
цесс разработки. 

Постановка задачи. В рамках данной работы рассматривается 
определение полярности текста — позитивный или негативный. Эмоци- 
онально нейтральные сообщения не учитываются для упрощения про- 
цесса обучения. Вследствие чего можно определить исходную задачу 
как задачу классификации. В работе использовались только русскоязыч- 
ные тексты. 

В качестве входных параметров для обучения и тестирования 
нейронной сети будут служить 2 набора данных: 

1. RuTweetCorp. Это набор, состоящий из отзывов на платформе 
Twitter. В него входят около 115000 положительных, 112000 отрица- 
тельных и 17,5 миллиона неразмеченных текстовых сообщений [1]. 
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2. ЕКТС. Набор размеченных данных, включающий в себя около 
14000 отрицательных отзывов с платформ 2ch.hk и pikabu.ru. Он отли- 
чается сильной эмоциональной окраской сообщений, с использованием 
ненормативной лексики [2]. 

Инструменты для разработки. В качестве языка для разработки 
был использован Python, так как он обладает большим количеством го- 
товых решений, применимых к нашей задачи. Для построения модели и 
обучения нейронной сети использовалась платформа Google 
Colaboratory. Этот облачный сервис предоставляет возможность бес- 
платно воспользоваться средой разработки, подготовленной для работы 
с Руфоп. Она обладает современными и производительными аппарат- 
ными средствами, что позволяет выполнять работу по подготовке дан- 
ных и обучению с высокой производительностью. 

Для работы с нейронной сетью использовалась библиотека Keras. 
При подготовке данных используются инструменты из библиотек Рап- 
das, Numpy и Sklearn. 

Подготовка входных данных. Для более удобного представления 
данные были модифицированы следующим образом: 

* удалены все знаки пунктуации в текстах; 

+ заменены упоминания пользователей Twitter на токен «User»; 

+ заменены ссылки, которые могут содержаться в тексте, на токен 

«URL»; 

* заменена буква «6» на «e», для уменьшения количества различ- 

ных вариаций одного и того же слова); 

ф все тексты приведены к нижнему регистру. 

Далее было проведено форматирование данных для представления 
текста в цифровом виде. Токенизация происходила на уровне слов и ис- 
пользовала алгоритм, согласно которому слова в предложениях заменя- 
лись на их индексы в словаре. После этого происходит процесс векто- 
ризации в формате плотного векторного представления (embedding). На 
рисунке 1 представлено отображение слов на плоскости в данном фор- 
мате. Если вместо выбранного вектора использовать унарный тип век- 
торов (one hot encoding), то получатся большие разряженные вектора. 
Современные СРО и СРО работают с ними недостаточно производи- 
тельно, а также данный тип занимает большее пространство в памяти, 
по сравнению с разряженными векторами [3]. 
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Рис. 1. Относительное расположение слов 
в формате плотного векторного представления 


Построение модели нейронной сети. Для данной задачи наибо- 
лее подходит тип рекуррентных нейронных сетей. С его помощью 
можно анализировать текст как последовательность токенов, в которой 
будет учтен порядок слов в предложении. Сети LSTM обладают невос- 
приимчивостью к длительности временных разрывов, благодаря чему 
имеют преимущество при решении поставленной проблемы [4]. 

В качестве основы была создана последовательная модель буду- 
щей сети. Первым добавлен слой плотных векторных представлений, 
Embedding. Длина вектора для представления слова выбрана 64. После 
него в модель добавляется слой с долгой краткосрочной памятью 
(LSTM) с размерностью выходного пространства 128. Рекуррентные 
сети склонны к переобучению, вследствие чего далее идет слой для ре- 
гуляризации — Dropout. С его помощью случайно выбранные нейроны 
игнорируются во время обучения. Вероятность исключения нейрона на 
каждом цикле обновления веса установлена в 10 %. В конце добавляется 
полносвязанный слой Dense с одним нейроном, который будет выпол- 
нять функцию классификации. Функция активации - сигмоида. 

Итоговая модель была сконфигурирована с бинарной кросс-энтро- 
пией в качестве функции ошибки и функцией оптимизации Адат. 

Обучение нейронной сети. С учетом того, что данные представ- 
лены в формате плотного векторного представления, на первом этапе 
элементы вектора инициализируются случайным образом. Далее ис- 
пользуется обучение с учителем, метод обратного распространения 
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ошибки, в процессе которого значения векторов изменяются итераци- 
онно. Обучение происходило в течение 10 зпох. Размер вариационной 
выборки составляет 10 %. Группы примеров, используемых для обуче- 
ния, состоят из 256 элементов. В качестве метрики качества обучения 
был выбран параметр точности на вылидационных данных [5]. Процесс 
обучения отображ ен на рисунке 2. 
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Рис. 2. График обучения на протяжении 10 зпох 


Судя по графику, можно сделать вывод о том, что полностью 
предотвратить переобучение не получилось. Доля верных ответов на 
обучающем наборе стабильно растет, в то время как доля верных отве- 
тов на проверочном наборе уменьшается после третьей зпохи. 

Выводы. В итоге была разработана нейронная сеть, с вероятно- 
стью верных ответов равной 77 Фо. Для достижения лучших результатов 
могут быть рассмотрены различные методы предотвращения переобу- 
чения рекуррентных сетей. Полученная модель может использоваться в 
качестве предобученной, сокращая время и требования к количеству 
размеченных данных при обучении новой нейронной сети. 
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APPLICATION OF RECURRENT NEURAL NETWORKS 
TO DETERMINE THE EMOTIONAL COLORING OF THE TEXT 


The paper considers the use of recurrent neural networks to determine the 
emotional coloring of a text. A set of data for training is presented and a method 
for their presentation in digital form is described. The process of constructing 
a neural network model using the LSTM architecture is considered. The train- 
ing of the obtained model is described and conclusions are drawn from the re- 
sults of work on the verification data. 
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ПРОЕКТИРОВАНИЕ БАЗЫ ДАННЫХ 
ДЛЯ МЕДИЦИНСКОЙ СИСТЕМЫ 


ПОДДЕРЖКИ ПРИНЯТИЯ РЕШЕНИЙ 


В работе проводится исследование проблемы проектирования 
базы данных в системе поддержки принятия решений, предназначен- 
ной для консультирования в сфере медицины. Приведена типовая 
схема взаимодействия между базой данных и пользователем. Разра- 
ботана объектная модель предметной области, позволяющая учи- 
тывать всю необходимую для принятия решений информацию. 
Сформулирована совокупность требований к медицинской базе дан- 
ных и приведен пример ее реализации. Разработана политика управ- 


ления доступом на основе ролевой модели. 


Ключевые слова: медицинская база данных, объектная модель 
предметной области, ролевая модель управления доступом. 


Введение. Проектированию медицинской базы данных предше- 
ствует моделирование предметной области, для чего необходимо пред- 
варительно провести анализ с учетом современных потребностей в ин- 
формации при принятии врачебных решений и выделить совокупность 
задач системы поддержки принятия решений [1]: 


ф 
ф 


ф 


описание диагноза; 
классификация в соответствии с современными стандартами и 
методиками; 

описание существующих схем лечения, а также выбор наиболее 
оптимальной из них; 

анализ совместимости различных методов лечения; 

сбор информации о потенциально возможных побочных дей- 
ствиях на пациента; 

поддержка принятия решения при выборе метода лечения. 


Следовательно, принятие врачебного решения тесно связано с 
процедурами поиска и перебора вариантов лечения. Если обобщить и 
формализовать процедуру принятия врачебного решения, то можно 
представить это в виде следующих этапов [2]: 


ф 


описание состояния пациента (проведение анализов, поста- 
новка диагноза); 

анализ ограничений (наличие аллергии, нарушение функций 
подсистем организма и т. п.); 

выбор метода лечения/лекарственного препарата; 

прогноз вариантов результата лечения (оценка безопасности, воз- 
можности возникновения нежелательных побочных реакций). 
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Согласно теории принятия решений, необходимо идентифициро- 
вать совокупность критериев выбора, множества альтернатив и алго- 
ритма выбора альтернатив согласно заданным критериям. В поставлен- 
ной задаче можно выделить следующие критерии [3]: 

* диагноз пациента (согласно международной классификации бо- 

лезней); 

+ схема лечения (множество, состоящее из элементов — методов 

лечения); 

+ наличие аллергических реакций у пациента. 

Рассмотренные принципы построения и особенности реализации 
системы поддержки принятия решений позволили представить следую- 


щую структуру взаимодействия между базой данных и пользователем 
(рисунок 1). 
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Рис. 1. Схема взаимодействия между базой данных и пользователем 
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Объектная модель предметной области. Система состоит из 3 
основных слоев [4]: 

Ф слой хранения данных, в котором хранится рассмотренная 

выше информация об объектах предметной области; 

Ф слой основной бизнес-логики, в котором хранятся правила ло- 

гического вывода, специфические для каждой подсистемы; 

Ф слой представления данных, отображающий пользовательский 

интерфейс - набор необходимых элементов управления. 

На основе выявленных особенностей предметной области выпол- 
нено ее моделирование. При помощи унифицированного языка модели- 
рования UML разработана объектная модель, позволяющая учитывать 
всю необходимую для принятия решения информацию (рисунок 2). 
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Рис. 2. Объектная модель предметной области 


Требования к реляционной базе данных. При переносе разрабо- 


танной объектной модели в реляционную СУБД необходимо обеспе- 
чить соответствие следующим требованиям к схеме данных [5]: 


Ф схема данных должна давать полное представление о предмет- 
ной области; 

* в схему данных должны быть включены все необходимые для 
выполнения запросов отношения (таблицы) и их атрибуты 
(столбцы) соответственно; 

Ф названия отношений должны быть уникальными; 

+ названия атрибутов в пределах одного отношения должны быть 
уникальными; 

Ф необходимо гарантировать однозначную трактовку схемы данных; 

* каждое отношение должно содержать первичный ключ (иден- 
тифицирующий набор атрибутов); 

Ф схема данных должна быть гибкой, при появлении необходимо- 
сти в новых запросах должна дополняться без масштабных пре- 
образований существующей схемы данных; 

Ф список отношений (таблиц) должен быть минимальным; 

Ф отношение используется тогда и только тогда, когда без него 
невозможно выполнить запрос; 
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список атрибутов должен быть минимальным; 
+ атрибут включается в отношение тогда и только тогда, когда без 
него описание предметной области не будет исчерпывающим; 
* первичный ключ отношения должен быть минимальным, невоз- 
можно исключить ни один атрибут из идентифицирующего 
набора, не нарушив при этом однозначной идентификации. 
Перечень запросов, выполнение которых необходимо обеспечить, 
включает в себя следующие [6]: 
* идентификация пользователя; 
Ф подключение к архивным данным (история болезни); 
* формирование перечня пациентов; 
* формирование исчерпывающей информации о конкретном па- 
циенте; 
добавление результатов анализа/осмотра; 
запись/изменение медицинской информации о пациенте; 
формирование шаблонов; 
подключение к вспомогательной информации (лабораторные 
данные, лекарственные препараты, международная классифи- 
кация болезней, графические и видеофайлы). 
Фрагмент схемы данных, обеспечивающей соответствие всем 
представленным выше требованиям, представлен на рисунке 3. 
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Рис. 3. Фрагмент схемът данных 


Политика управления доступом. Наконец, для обеспечения безопас- 
ной работы с данными необходимо разработать политику управления досту- 
пом к данным в соответствии с одной из перечисленных моделей [7-9]. 

Избирательное управление доступом — это управление доступом 
субъектов к объектам на основе списков управления доступом или мат- 
рицы доступа. Для каждой пары «субъект — объект» должно быть задано 
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явное и недвусмысленное перечисление допустимых типов доступа (чте- 
ние, запись и т. д.), то есть тех типов доступа, которые являются санкцио- 
нированными для данного субъекта (пользователя) к данному объекту (ре- 
сурсу базы данных). Поскольку при функционировании базы данных ко- 
личество пользователей и ресурсов будет исчисляться в пределах порядка 
от 10? до 103, перечисление допустимых типов доступа для каждой пары 
«субъект — объект» является неоправданно громоздкой задачей. 

Мандатное управление доступом — это разграничение доступа 
субъектов к объектам, основанное на назначении метки конфиденциаль- 
ности для информации, содержащейся в объектах, и выдаче официаль- 
ных разрешений субъектам на обращение к информации такого уровня 
конфиденциальности. Такие метки доступа являются ключевым отли- 
чием систем защиты государственной тайны и не подходят для меди- 
цинских систем, где необходима более гибкая настройка. 

Управление доступом на основе ролей - это развитие политики из- 
бирательного управления доступом, при этом права доступа субъектов 
системы на объекты группируются с учетом специфики их применения, 
образуя роли. Таким образом, отсутствует необходимость задания прав 
доступа для каждой пары «субъект — объект», характерная для избира- 
тельного управления доступом, поскольку пользователей можно сгруп- 
пировать по ролям в соответствии с должностью в медицинском учре- 
ждении. Также отсутствует жесткая привязка к меткам доступа, харак- 
терная для мандатного управления доступом, в связи с чем ролевое 
управление доступом обеспечивает более гибкую настройку. 

В политике управления доступом на основе ролевой модели опре- 
делено, к каким объектам будет иметь доступ каждый субъект и какие 
действия может выполнять (таблица 1): 

1) К— чтение; 

2) № – запись/модификация; 

3) С – создание; 

4) D – удаление. 


Таблица 1 
Политика управления доступом 
Роль кн ы Данные Данные | Системные 
персонала | пациентов| данные 
учреждения 
Администратор CRWD CRWD 
Главный врач RW CRWD CRW 
Заместитель главного R CRWD CRW 
врача 
Заведующий отделением R RW CRW 
Дежурный врач К В CRW 
Bpau R CRW 
Медсестра R CRW 
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Выводы. 

1. Обеспечение гибкого взаимодействия между базой данных и 
пользователем является основной задачей проектировщика медицин- 
ской системы поддержки принятия решений. 

2. Грамотность проектирования медицинской базы данных зависит 
от полноты объектной модели предметной области и соблюдения совокуп- 
ности требований при ее преобразовании в реляционную модель. 

3. Информационная безопасность медицинской базы данных до- 
стигается с помощью разработки политики управления доступом на ос- 
нове ролевой модели, которая в данной области обладает преимуще- 
ствами по сравнению с избирательной и мандатной моделями. 
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A DATABASE DEVELOPMENT 
FOR A MEDICAL DECISION SUPPORT SYSTEM 


The paper investigates the problem of a database development in a decision 
support system designed for consulting in the field of medicine. A typical scheme 
of interaction between the database and the user is given. An object model of 
the domain is developed, allowing to take into account all the information nec- 
essary for decision-making. A set of requirements for a medical database is 
formulated and an example of its implementation is given. An access control 
policy based on a role model is developed. 
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УДК 004.624 
А.Н. Молчанов, Д.С. Серпинская 


СПОСОБЫ ПРОВЕРКИ ПОДЛИННОСТИ 
ЭЛЕКТРОННОЙ ПОДПИСИ 


Отмечается актуальность необходимости защитных атрибу- 
тов для юридически важных документов, передающихся в электрон- 
ном виде. Раскрывается сущность электронной подписи. Дается 
определение понятиям, относящимся к электронной подписи. Пере- 
числяются существующие виды электронной подписи и их отличия 
друг от друга. Обращается внимание на то, что сертифицирован- 
ные электронные подписи выдаются владельцам в сертифицирован- 
ных Удостоверяющих центрах. Рассматриваются различные сред- 
ства проверки подлинности злектронной подписи. Приводится по- 
дробное описание каждого из них. При этом подчеркивается, что 
выбор конкретного метода зависит от возможностей и предпочте- 
ний участника электронного документооборота. 


Ключевые слова: электронная подпись, открытый ключ, закры- 
тый ключ, сертификат, удостоверяющий центр, криптопровайдер. 


Введение. В эпоху глобализации мировой экономики активное 
развитие получила сфера информационных технологий. Одним из 
наиболее важных изобретений этой эпохи, которое навсегда изменило 
жизнь человечества — это создание распределенных компьютерных се- 
тей (Интернет). Данный факт поспособствовал оптимизации взаимодей- 
ствия людей с окружающим миром. Это дало возможность передавать 
различные документы в электронном виде. Однако для обеспечения без- 
опасного юридически-значимого документооборота появилась необхо- 
димость обеспечения защиты прав владельца передаваемых докумен- 
тов, а также защиты непосредственно передаваемых документов от мо- 
дификации [1]. 

Одним из таких средств защиты является электронная подпись (ЭП). 

Понятия электронной подписи. Электронная подпись позволяет 
установить отсутствие искажения информации в электронном доку- 
менте с момента формирования ЭП и подтвердить принадлежность ЭП 
ее владельцу. Принцип работы ЭП основан на применении алгоритма 
асимметричного шифрования. Данный вид шифрования подразумевает 
использование связи пары ключей: открытого и закрытого. Они пред- 
ставляют собой последовательность символов, генерация которых про- 
исходит в специальных шифровальных программах [2]. 

Закрытый ключ хранится у владельца подписи в закрытом до- 
ступе, как правило, на цифровом носителе: токене или смарт-карте. Для 
подписываемого документа формируется хеш-значение с помощью 
криптографических модулей и данный ключ его шифрует. Открытый 
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ключ, наоборот расшифровывает подписанный соответствующей злек- 
тронной подписью документ. Данный ключ передается в открытом виде 
всем участникам электронного документооборота. 

При получении ЭП в Удостоверяющем центре ее владельцу выда- 
ется сертификат ЭП. ОН содержит в себе информацию об удостоверяю- 
щем центре, который его выдал, о дате создания и окончания срока дей- 
ствия подписи, а также данные о ее владельце. Данную часть подписи 
используют для ее верифицирования [3]. 

Согласно Федеральному закону от 06.04.2011 М 63-ФЗ (ред. от 
23.06.2016) "Об электронной подписи" существуют несколько видов 
электронных подписей [4]: 

ф простая; 

* усиленная. 

Простой можно называть такую электронную подпись, которая 
подтверждает факт формирования электронной подписи определенным 
лицом посредством использования паролей, кодов или иных программ- 
ных или технических средств. 

Усиленная электронная подпись в свою очередь делится на неква- 
лифицированную и квалифицированную. 

Усиленной неквалифицированной электронной подписью счита- 
ется та, которая обладает следующими признаками [5]: 

Ф является результатом криптографического преобразования ин- 

формации с применением ключа электронной подписи; 

Ф позволяет идентифицировать личность человека, который под- 

писал электронный документ; 

* позволяет обнаружить факт наличия модификации в электрон- 

ном документе с момента его подписания; 

Ф создается с помощью средств электронной подписи. 

Для данного вида подписи не обязателен факт наличия сертифи- 
ката ключа в случае, если соответствие электронной подписи признакам 
неквалифицированной электронной подписи, установленным настоя- 
щим Федеральным законом, может быть обеспечено без использования 
сертификата ключа проверки электронной подписи. 

Усиленная квалифицированная электронная подпись - это элек- 
тронная подпись, которая соответствует всем признакам неквалифици- 
рованной электронной подписи, но при этом имеет дополнительные 
признаки: 

* вквалифицированном сертификате указан ключ проверки элек- 

тронной подписи указан; 

Ф для создания и проверки электронной подписи используются 

такие средства электронной подписи, которые имеют подтвер- 
ждение соответствия требованиям, установленным согласно с 
настоящим Федеральным законом. 

Сертифицированные электронные подписи выдаются владельцам 
в сертифицированных Удостоверяющих центрах. 
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В настоящее время уполномоченным органом Российской Федера- 
ции, который осуществляет аккредитацию удостоверяющих центров, 
проводит проверки соблюдения ими требований, установленных насто- 
ящим Федеральным законом и иными принимаемыми в соответствии с 
ним нормативными правовыми актами, является Министерство цифро- 
вого развития, связи и массовых коммуникаций Российской Федерации 
(МИНКОМСВЯЗЬ) [6]. 

Порядок верификации подписи, используемой для подписанного до- 
кумента, состоит в проверке данных сертификата используемой подписи: 

1. Уникального номера квалифицированного сертификата, а также 
даты начала и окончания срока его действия; 

2. Личных данных, указанных в сертификате: 

* для физического лица, не являющегося индивидуальным пред- 
принимателем, производится проверка фамилии, имени и отче- 
ства (если имеется) владельца квалифицированного сертификата. 

* для физического лица, являющегося индивидуальным предприни- 
мателем — фамилия, имя, отчество (если имеется) и основной гос- 
ударственный регистрационный номер индивидуального пред- 
принимателя — владельца квалифицированного сертификата; 

* для российского юридического лица - наименование, место 
нахождения и основной государственный регистрационный но- 
мер владельца квалифицированного сертификата; 

* для иностранной организации (в том числе филиалов, предста- 
вительств и иных обособленных подразделений иностранной 
организации) — наименование, место нахождения владельца 
квалифицированного сертификата, а также идентификацион- 
ный номер налогоплательщика (при наличии). 

3. Такие юридические данные как: 

* для физического лица — страховой номер индивидуального ли- 
цевого счета и идентификационный номер налогоплательщика 
владельца квалифицированного сертификата; 

* для юридического лица — идентификационный номер налого- 
плательщика владельца квалифицированного сертификата. 

4. Уникальный ключ проверки электронной подписи; 

5. Наименования средств электронной подписи и средств аккреди- 
тованного удостоверяющего центра, которые были использованы для 
создания ключа электронной подписи, ключа проверки электронной 
подписи, квалифицированного сертификата, а также реквизиты доку- 
мента, который подтверждает соответствие указанных средств тем тре- 
бованиям, которые установлены в соответствии с настоящим Федераль- 
ным законом; 

6. Полное наименование и место нахождения аккредитованного 
удостоверяющего центра, который выдал квалифицированный сертифи- 
кат и номер квалифицированного сертификата удостоверяющего центра; 

7. Ограничения использования квалифицированного сертификата 
(если такие ограничения устанавливаются) [7]. 
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Средства для верификации электронной подписи. На данный 
момент существуют различные средства для проверки подлинности 
злектронной: 

* на официальном портале государственных услуг; 

Ф с помощью специального программного обеспечения (напри- 
мер, программы «КриптоПро» и аналогичных криптографиче- 
ских комплексов); 

+ средствами программного комплекса Microsoft Word / PDF или 
по значениям хеш-функций (однако они требуют наличия на 
компьютере пользователя специализированных утилит); 

* автоматическая проверка документов внутри сервисов, которые 
непосредственно связаны с документооборотом. Например, 
сервисы операторов фискальных данных (ОФД) и электронного 
документооборота (ЭДО). 

Самым простым методом проверки такой подписи является ис- 
пользование данной услуги на официальном портале государственных 
услуг. Данный способ подходит тем пользователям, которые не имеют 
возможности или желания применять стороннее программное обеспече- 
ние (ПО) или сложные методы проверки, однако хотят произвести вери- 
фикацию: 

* отсоединенной электронной подписи в формате РКС$#7; 

* сертификата электронной подписи; 

* сертификата открытого ключа. 

Для такого вида проверки пользователю необходимо посетить 
официальный сайт государственных услуг, в главном окне данного 
электронного ресурса выбрать документы для уточнения, загрузить 
нужный сертификат и нажать на соответствующую кнопку для начала 
проверки (рисунок 1). 

Однако, такой метод верификации ЭП позволяет установить под- 
линность только сертифицированной электронной подписи, сертификат 
которой содержится в общем реестре открытых ключей, действующих 
в России электронных подписей. 


Подтверждение подлинности ЭП: 


— сертификата 


Подтверждение подлинности ЭП Справо-нал индрормашия по 
сертификата подписа (РО. ЗОВ) 


Данный сервис предназначен для подтверждения подлинности ЭП сертификатов Аккредитованные 
в форматах Х 509 и BASE64. удостоверяющие центры 


Можно подтвердить подлинность ЭП сертификата, изданного удостоверяющим 
центром, входящим в список аккредитованных удостоверяющих центров 
Министерства связи и массовых коммуникаций 


В поле «Выберите сертификат для проверки» следует выбрать сертификат, подлинность ЭП которого следует подтвердить, и 
нажать кнопку «Проверить». На экране будет отображена информация о результатах проверки сертификата 


Предоставляемая услуга носит информационный характер и не может быть использована в качестве доказательств в 
судах различных инстанций. 


" выберите сертификат для проверки: 
загрузить файл 
Рис. 1. Проверка злектронной подписи 
с помощью портала государственных услуг России 
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В случае, если пользователь обладает средствами для приобрете- 
ния специализированного криптографического ПО (криптовровайдера), 
то выбор конкретного уже состоит из выделяемых средств на его при- 
обретение и типа операционной системы, на которую будут установлена 
криптографическая программа. Программы для шифрования данных 
для своей работы требуют наличие на компьютере пользователя крип- 
топровайдера, и в таком случае дают возможность не только осуществ- 
лять проверку электронных подписей, но и подписывать документы 
имеющимися у пользователя собственными электронными подписями, 
а также осуществлять шифрование различных файлов или даже диско- 
вых пространств. Такие лицензионные продукты можно приобрести на 
официальных сайтах их разработчиков. 

На данный момент популярными программами для шифрования 
являются: 

+ УегаСгурб 

* АхСгурб 

* Folder Lock; 

* CryptoExpert 8; 

* CertainSafe. 

А в качестве производителей криптопровайдеров часто предпочи- 
тают выбирать: 

* ООО «КРИПТО-ПРО»; 

* ООО «Лисси»; 

Ф ОАО «ИнфоТеКС»; 
* ЗАО «Сигнал-КОМ»; 
* ОКБ САПР; 

* AVEST CSP; 

* КриптЭК; 

+ Криптолайн. 

При установленном на компьютере криптопровайдере y пользова- 
телей есть возможность произвести проверку злектронной подписи с 
помощью программы MicrosoftWord/PDF. Для этого необходимо: нали- 
чие компьютера с лицензионной операционной системой Windows 7 и 
выше, точно установленные время и дата на компьютере, настроенная 
на работу с ЭП лицензионная версия Microsoft Word 2007 и выше, 
настроенный лицензионный криптографический продукт, например, 
КриптоПро CSP, установленный модуль КриптоПро Office Signature, а 
также доступ в Интернет. 

В случае, если все необходимые требования учтены, пользователь 
может открыть подписанный документ в одной из программ 
MicrosoftWord или PDF. Над отображением данных в документе будет 
выведено сообщение о том, что документ содержит электронную под- 
пись и кнопку для ее открытия информации о сертификате открытого 
ключа (СОК). На данном этапе подпись может имееть статус недействи- 
тельной. В данном случае у пользователя есть возможность произвести 
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проверку подлинности данной подписи двумя способами: с помощью 
специального ресурса на официальном сайте КриптоПро или «вруч- 
ную». Для обоих случаев первоначально необходимо сохранить СОК 
как отдельный документ. Далее первый способ требует перейти на озву- 
ченный ранее сайт, открыть ссылку «проверить ЗЦП сертификата», за- 
грузить нужный сертификат с компьютера и нажать проверить. Если 
при этом электронная подпись была подтверждена, то полученную стра- 
ницу надо распечатать и приложить к котировочной заявке. 

Для «ручной» проверки ЭП пользователь должен открыть СОК 
встроенными средствами операционной системы — двойным щелчком 
на сохраненном сертификате. В появившемся окне найти наименование 
удостоверяющего центра (УЦ), выпустившего СОК, точку распростра- 
нения списка отозванных сертификатов. Для проверки действительно- 
сти корневого сертификата (КС) УЦ необходимо перейти на сайт под- 
системы удостоверяющих центров и, используя данные из СОК через 
диалог поиска, по имени УЦ найти необходимую электронную подпись, 
сохранить все файлы КС и установить на компьютер только те, которые 
ранее отсутствовали в системе. Далее скачать и установить список ото- 
званных сертификатов, перейдя по адресу ОКТ, ресурса, который указан 
в сертификате. На следующем шаге пользователю необходимо с помо- 
щью MS Office открыть подписанный ЭП документ и проверить статус 
подписи. На зтом зтапе у подписанного ЗП документа будет отобра- 
жаться текущее состояние использованной для отправки данных злек- 
тронной подписи [8]. 

Однако все перечисленные выше способы верификации ЗП под- 
ходят для людей, которые редко взаимодействуют с электронным доку- 
ментооборотом. Пользователи, которые постоянно взаимодействуют с 
получением и передачей электронных документов, обычно пользуются 
услугами специализированных сервисов, например, операторы фис- 
кальных данных (ОФД) или сервисы электронного документооборота 
(ЭДО). Данные службы имеют встроенные и специально разработанные 
модули для подписания и проверки проходящих через них документов. 
В случае необходимости подписать передаваемые документы данные 
модули автоматически считывают с компьютера клиента данные мо- 
дули и незаметно для пользователя производят проверку поступающих 
документов. 

В сервисе Астрал.ОФД используется следующий алгоритм про- 
верки подписи: 

1. При поступлении в сервис документа, подписанного электрон- 
ной подписью, криптографический модуль обращается полученному 
файлу и извлекает из него электронную подпись. 

2. Для взаимодействия с электронной подписью специальный мо- 
дуль обращается к криптопровайдеру, с его помощью извлекает из нее 
ту часть, которая содержит ее сертификат. 
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3. Далее модуль обращается к хранилищу, где содержатся все вы- 
данные удостоверяющими центрами действующие сертификаты, поме- 
щает в него открытую часть сертификата ЭП, которую нужно прове- 
рить, и сопоставляет его с содержащимися в реестре сертификатами. 

Сначала происходит поиск соответствующего сертификата из 
списка по значению объектного идентификатора области использования 
ключа. Затем сопоставляются значения полей сертификатов согласно 
структуре стандарта Х509 для открытого ключа. 

В случае нахождения соответствия сертификатов происходит 
сравнение хеш-значения данных в подписях с помощью российского ал- 
горитма шифрования ГОСТ Р 34.10-2012. 

4. Если данные совпали происходит проверка документа на целост- 
ность данных. Программа снова обращается к полученному документу, 
вычисляет для него значение хеш-функции, с помощью открытого ключа 
электронной подписи зашифровывает его и сравнивает со хеш-значением 
документа, которое содержится в прикрепленной к нему ЭП. 

Если все проверяемые данные соответствуют друг другу, то ЭП 
считается действительной. Данная проверка происходит «налету», а 
пользователь сервиса получает лишь результат проверки ЭП документа, 
который содержит данные о ее владельце и статус актуальности. 

Заключение. Таким образом, существует достаточно способов 
для получения данных об электронных подписях, использованных для 
передаваемых документов, и проверки их подлинности. Выбор конкрет- 
ного способа для этого зависит от возможностей самого пользователя и 
его личных предпочтений. Для людей, постоянно взаимодействующих 
с юридически значимым электронным документооборотом, одним из 
наиболее удобных способов обеспечения безопасного электронного до- 
кументооборота является использование услуг специализированных 
сервисов. 
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The urgency of the need for protective attributes for legally important doc- 
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УЛУЧШЕНИЕ КАЧЕСТВА МОДЕЛЕЙ 

МАШИННОГО ОБУЧЕНИЯ В ЗАДАЧАХ 
КЛАССИФИКАЦИИ ИЗОБРАЖЕНИЙ 

НА ОСНОВЕ ПОДХОДОВ ИЗВЛЕЧЕНИЯ ПРИЗНАКОВ 
И ТОЧНОЙ НАСТРОЙКИ МОДЕЛИ 


В данной статье указывается основная цель и особенности под- 
хода передачи знаний, а также приводятся основные стратегии, ха- 
рактерные для данного подхода. Кроме того, раскрывается сущность 
таких методов, как точная настройка модели и извлечение признаков 
и приводится аргументация того, почему для извлечения признаков 
нужно использовать свёрточные нейронные сети, учитывая их архи- 
тектурные особенности. Также приводится обоснование необходи- 
мости в точной настойке моделей глубокого обучения. Помимо этого, 
в данной статье рассматриваются основные техники, используемые 
для точной настройки моделей машинного обучения. В заключении де- 
лается вывод о перспективности использования приемов извлечения 
признаков и точной настройки модели, а также выделяются ключе- 
вые момеиты дальнейшего развития данных подходов. 


Ключевые слова: передача знаний, извлечение признаков, точная 
настройка, классификация изображений. 


Введение. Как известно, человек при решении различных задач накап- 
ливает опыт и стремится перенести его на решение новых схожих задач. 
Иными словами, люди получают и применяют соответствующие знания из 
предыдущего опыта обучения, когда сталкиваются с новыми задачами [1]. 
Чем больше новая задача связана с предыдущим опытом, тем легче им спра- 
виться с этим. В свою очередь стандартные алгоритмы машинного обуче- 
ния, напротив, традиционно решают отдельные задачи [2, 3]. 

Такой подход, как передача знаний пытается исправить это путем 
разработки методов передачи знаний (признаков), полученных в одной 
или нескольких исходных задачах, и их использования для повышения 
качества обучения в связанной целевой задаче. Методы, позволяющие 
передавать знания, несут прогресс в том, чтобы сделать машинное обу- 
чение столь же эффективным, как и человеческое обучение. Методы пе- 
редачи знаний, как правило, сильно зависят от алгоритмов машинного 
обучения, используемых для изучения задач, и часто могут просто рас- 
сматриваться как расширение этих алгоритмов. Некоторая работа по пе- 
редаче знаний проводится в контексте индуктивного обучения и вклю- 
чает расширение хорошо известных алгоритмов классификации и логи- 
ческого вывода, таких как нейронные сети, байесовские сети и марков- 
ская логика [4]. 

Целью передачи знаний является улучшение качества обучения в 
целевой задаче путем использования знаний из исходной задачи [5]. 
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Данная техника может улучшить обучение по трём наиболее распро- 
страненным метрикам. Во-первых, это начальная эффективность, дости- 
жимая в целевой задаче с использованием только переданного знания, 
прежде чем будет выполнено какое-либо дальнейшее обучение, по срав- 
нению с начальной эффективностью обучаемой модели без этих знаний. 
Во-вторых, количество времени, необходимое для полного изучения це- 
левой задачи с учетом переданного знания по сравнению с количеством 
времени, чтобы изучить его с нуля. В-третьих, это конечный уровень 
эффективности, достижимый в целевой задаче, по сравнению с конеч- 
ным уровнем без передачи [6]. 

Передача знаний. Существует две основных и наиболее популяр- 
ных стратегии передачи знаний: 

1. Использование обученной модели для извлечения признаков. В 
задачах классификации изображений одним из наиболее применяемых 
для этой цели методов являются сверточные нейронные сети (CNN). 
СММ относятся к специфичной категории нейронных сетей [7]. СММ мо- 
жет достаточно точно изучать изображения. К тому же она оказалась 
намного лучше многих традиционных методов обработки изображений. 
Стандартная сверточная нейронная сеть состоит из повторяющихся 
групп сверточных слоев, каждый из которых обладает одинаковым чис- 
лом фильтров; слоев пулинга и полносвязных слоёв. Последние функ- 
ционируют в качестве классификатора. Сверточные нейронные сети об- 
ладают тремя архитектурыми концепциями, которые позволяют им 
быть достаточно эффективными: локальные рецептивные поля, распре- 
деление весов и операции субдискретизации [8]. Первое означает, что 
каждый нейрон принимает входные данные только лишь от небольшой 
части нейронов предыдущего слоя. Локальные рецептивные поля ис- 
пользуются в сверточных слоях и слоях пулинга. Распределение весов 
используется в сверточных слоях для уменьшения сложности модели. 
Наконец, нелинейная субдискретизация (используемая в слоях пулинга) 
уменьшает пространственный размер изображения и количество сво- 
бодных параметров. В качестве примера можно взять исходную модель 
СОМУМЕТ, обученную на датасете ImageNet, удалить последний пол- 
носвязный слой (выходные данные этого слоя представляют собой про- 
гнозы для 1000 классов изображений), азатем использовать оставшуюся 
часть для получения признаков из нового набора данных (рис. 1). После 
извлечения признаков для всех изображений на этих признаках можно 
обучить собственный классификатор (например, SVM или Softmax 
классификатор). 

2. Настройка готовой модели. Вторая стратегия заключается не 
только в замене и переобучении на новом наборе данных классифика- 
тора из верхней части модели, но также и в точной настройке весов 
предварительно обученной нейронной сети путем продолжения обрат- 
ного распространения ошибки (рис. 2). Можно выполнить точную 
настройку всех слоёв сети или зафиксировать некоторые из более ран- 
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них слоёв (из-за проблем с переобучением) и настроить только высоко- 
уровневую часть сети. Это объясняется тем, что начальные слои содер- 
жат более общие признаки, которые могут быть полезны для многих за- 
дач. В свою очередь нижние слои модели постепенно становятся все бо- 
лее специфичными для деталей конкретных классов, содержащихся в 
исходном наборе данных. Например, в случае ImageNet, который содер- 
жит множество пород собак, значительная часть репрезентативных воз- 
можностей моделей, обученных на данном датасете, может обнаружи- 
вать особенности, которые характерны для различия между этими поро- 
дами. Более подробно основные техники точной настройки модели бу- 
дут рассмотрены далее в этой статье. 


Вход 
| Извлечение признаков 


Признаки 
— 


Входное изображение 


1-й Скрытый слой N 


2-й Скрытый слой 
і-й Скрытый слой 
п-Й Скрытый слой 


Рис. 1. Общая схема извлечения признаков 


Исходная Целевая 
модель модель 
Е Инициализация 
Выходной ч Ж Обучение 
слой произвольными 
весами с нуля 
Предварительно 
— копирование — 
обученная E 
модель Настройка 
А — — копирование — + р 
модели 


Исходные 
данные 


Рис. 2. Общая схема точной настройки модели 


Целевые 
данные 


Точная настройка модели. Перед рассмотрением методов точной 
настройки необходимо ответить на вопрос: почему вообще нужно 
настраивать модели? Когда на практике возникает задача глубокого 
обучения, например, задача, включающая в себя обучение сверточной 
нейронной сети (СММ) на некотором наборе изображений (с целью их 
дальнейшей классификации или распознавания), одним из самых пер- 
вых и очевидных способов является обучение какой-либо готовой мо- 
дели (с грамотно подобранной и проверенной архитектурой) с нуля. В 
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теории данный способ, несомненно, позволит получить на выходе хо- 
рошо обученную модель, ориентированную на работу с данными из рас- 
сматриваемой предметной области. Однако на практике существует ряд 
ограничений. Прежде всего стоит отметить, что глубокие нейронные 
сети, такие как СОУМЕТ или VGG-16, состоят из большого количества 
слоёв и имеют огромное число параметров [9]. Часто их количество до- 
стигает нескольких миллионов. Обучение подобных моделей на неболь- 
шом наборе данных может сильно повлиять на их способность обоб- 
щать и извлекать из данных только важные признаки, что часто приво- 
дит к переобучению и достаточно низкой точности на выходе. Кроме 
того, стоит учесть и тот факт, что на промежуточных слоях эти модели 
оперируют тензорами очень больших размерностей. Соответственно по- 
требуется наличие колоссальных объемов памяти (чаще всего речь идет 
о видеопамяти, поскольку обучение нейронных сетей на практике про- 
исходит именно на дорогостоящих графических картах). Помимо фи- 
нансовых затрат следует также принять во внимание и временной фак- 
тор. Ведь несмотря на использование высокопроизводительных графи- 
ческих процессоров время обучения «глубоких» моделей с нуля может 
занимать несколько дней или даже недель. 

Поэтому на практике возникает потребность в настройке суще- 
ствующих моделей, которые заранее были обучены на достаточно боль- 
шом наборе данных, таком как ImageNet (включает в себя порядка 1.2 
млн. помеченных изображений, относящихся к 1000 классам) с продол- 
жением их обучения (частичном) на меньшем наборе данных, который 
имеется для конкретной предметной области. При условии, что целевые 
данные не сильно отличаются от исходного набора данных (например, 
от того же ImageNet), предварительно обученная модель уже изучит 
признаки, которые имеют отношение к целевой задаче (распознавания 
или классификации). 

Рассмотрим основные техники точной настройки модели [10]: 

1. Обычной практикой является удаление последнего слоя (слоя 
softmax) предварительно обученной сети и замена его новым слоем 
softmax, который соответствует целевой задаче. Например, предвари- 
тельно обученная Ha датасете ImageNet сеть имеет на выходе слой 
softmax c 1000 категориями. Если целевой задачей, например, является 
классификация изображений по 10 категориям, новый слой softmax дол- 
жен состоять из 10 классов вместо 1000. Затем необходимо запустить 
обучение сети на целевых данных для точной настройки весов. При этом 
рекомендуется выполнять перекрестную проверку, чтобы сеть могла хо- 
рошо обобщать входные данные. 

2. Использование меньшей скорости обучения сети. Поскольку 
ожидается, что предварительно обученные веса уже будут достаточно хо- 
рошими по сравнению со случайно инициализированными весами, нет 
необходимости изменять их слишком быстро и слишком сильно. Распро- 
страненный подход - сделать начальную скорость обучения в 10 раз 
меньше, чем та, которая использовалась для тренировки модели с нуля. 
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3. Еще одним достаточно популярным подходом является «замо- 
раживание» весов на первых нескольких слоях предварительно обучен- 
ной модели. Это связано с тем, что первые несколько слоев, как было 
отмечено ранее, содержат универсальные признаки, которые также мо- 
гут иметь отношение к новой целевой задаче. Поэтому желательно со- 
хранить эти веса нетронутыми. Вместо этого лучше заставить сеть со- 
средоточиться на изучении специфичных для целевого набора данных 
признаков на последующих слоях. 

Варианты применения стратегий передачи знаний. Помимо 
принципов использования точной настройки необходимо также рас- 
смотреть возможные сценарии применения того или иного типа пере- 
дачи знаний. Соответствующий выбор зависит от нескольких факторов, 
но два наиболее важных из них — это объём целевого набора данных 
(малый или достаточно большой) и степень его схожести с исходным 
набором данных (рис. 3), к примеру, с точки зрения изображений и их 
классов. Принимая все это во внимание, можно выделить 4 основных 
варианта (рис. 4): 

Если целевой набор данных не сильно отличается по контексту от 
набора данных, на котором обучалась исходная модель, то в таком слу- 
чае нужно прибегнуть к использованию точной настройки, поскольку 
предварительно обученная на таком большом и разнообразном наборе 
данных, как ImageNet, нейронная сеть в своих начальных слоях фикси- 
рует общие детали, например, различные кривые и ребра, которые акту- 
альны и полезны для большинства задач классификации. Однако стоит 
учесть, что если целевой набор данных очень большой, то имеет смысл 
провести настройку весов во всей сети (или даже обучить сеть с нуля), 
а не только в конечных слоях. 

Несомненно, если целевой набор данных относится к очень специ- 
фической предметной области, например, медицинские изображения 
или китайские рукописные символы, и соответственно никакие предва- 
рительно обученные модели в таком случае не могут быть найдены, то- 
гда необходимо рассмотреть возможность обучения сети с нуля. 

Еще одна проблема заключается в том, что если целевой набор 
данных небольшой по размеру, то настройка предварительно обученной 
сети на таком наборе данных может привести к переобучению, особенно 
если последние несколько слоёв являются полносвязными, как в случае 
модели УСС-16. В таком случае можно попытаться использовать стан- 
дартные стратегии приращения данных (перемещение, вращение, пере- 
ворачивание и др.). Если в конечном итоге получится достаточно сильно 
«пополнить» данные, то применение точной настройки модели все еще 
может дать неплохие результаты. 

Если целевой набор данных действительно слишком мал, к при- 
меру, менее тысячи выборок, лучшим подходом является использование 
выходных данных одного из промежуточных слоёв, расположенного до 
полносвязных слоев, в качестве признаков для обучения линейного 
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классификатора (например, ЗУМ), то есть использование техники из- 


влечения признаков. 


Размер 
датасета 


Датасет большой, но 
сильно отличается от 
датасета, 
используемого для 
обучения исходной 
модели 


Датасет большой и 
схож с датасетом, 
используемым для 
обучения исходной 
модели 


Схожесть 
датасета 


Датасет малый по 
объему, и сильно 
отличается от 
датасета, 
используемого для 
обучения исходной 

модели 


Датасет малый по 
объему, но схож с 
датасетом, 
используемым для 
обучения исходной 

модели 


Рис. 3. Основные факторы целевого датасета 
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Рис. 4. Варианты использования готовой модели 
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Заключение 
В заключении хотелось бы отметить, что несмотря на все ограни- 


чения техника передачи знаний и, в частности, приемы извлечения при- 
знаков и точной настройки модели остаются очень перспективным 
направлением в области машинного обучения. Одной из важных задач, 
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которую предстоит решить является обеспечение возможности пере- 
дачи знаний между более разнообразными и отличающимися друг от 
друга задачами. 

Еще одна проблема заключается в осуществлении передачи в бо- 
лее сложных областях машинного обучения. Например, передача зна- 
ний в области обучения с подкреплением может стать гораздо труднее, 
поскольку исходные и целевые задачи становятся все более сложными. 
С практической точки зрения изучение данного вопроса, вероятно, ста- 
нет очень многообещающим направлением. Важно не ограничивать ис- 
следования по передаче знаний простыми доменами. 

Методы передачи знаний стали значительной областью изучений 
в машинном обучении. Они обладают многими идеологическими пре- 
имуществами, поскольку рассматриваются как компьютерный аналог 
важных аспектов человеческого обучения, а также практическими пре- 
имуществами, поскольку они способны сделать машинное обучение 
еще более эффективным. По мере увеличения вычислительной мощно- 
сти и применением исследователями машинного обучения для решения 
самых сложных задач, передача знаний будет становиться все более вос- 
требованной и широко применяемой техникой. 
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IMPROVING THE QUALITY OF MACHINE LEARNING 
MODELS IN THE PROBLEMS OF IMAGE CLASSIFICATION 
BASED ON THE APPROACHES OF FEATURE EXTRACTION 

AND FINE-TUNING THE MODEL 


This article indicates the main goal and features of the transfer learning 
method, and also describes the main strategies specific to this method. Further- 
more, the essence of methods such as fine-tuning the model and features extrac- 
tion is revealed, and an argumentation is given as to why convolutional neural 
networks should be used to extract features, given their architectural features. 
The rationale of the necessity for the exact tune-up of deep learning models is 
also provided. In addition, the article considers the basic techniques used to 
fine tune machine learning models. At the same time, key features of each ap- 
proach are noted. In conclusion, the article concluded the use of techniques for 
extracting attributes and fine-tuning the model is promising, also key points in 
the further development of these approaches are highlighed. 


Keywords: transfer learning, feature extraction, fine-tuning, image classifi- 
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РУБРИКА 3 


ПРИБОРОСТРОЕНИЕ И ЭЛЕКТРОНИКА 


САТЕСОКУ 3 


DEVICES AND ELECTRONICS 


Электронный журнал: наука, техника и образование (ISSN 2413-6220) 113 


УДК 621.382 
В.В. Андреев, А.Д. Перминова 


ИНТЕГРАЛЬНЫЙ ДАТЧИК ТЕМПЕРАТУРЫ 


Данная статья посвящена разработке интегрального однокри- 
стального полупроводникового датчика температуры, на выходе ко- 
торого вырабатывается напряжение пропорциональное абсолют- 
ной температуре и поиску путей повышения точности измерения 
температуры. Получены экспериментальные графики выходного 
напряжения, на основе которых скорректированы конструкция и 
топология делителя напряжения, что позволило повысить точность 
измерения температуры. 


Ключевые слова: Датчик температуры, интегральная микро- 
схема, температурный коэффициент, точность измерения. 


Введение. Для повышения надежности электронной аппаратуры 
обычно используют интегральные микросхемы, элементы которых об- 
ладают повышенной стойкостью к критическим воздействиям, таким 
как: ионизирующее излучение, сильные электрические поля, изменения 
температуры ит. д. [1-7]. Интегральные датчики температуры (ИДТ) яв- 
ляются неотъемлемой частью многих электронных устройств, так как 
обеспечивают требуемые характеристики аппаратуры в широком диапа- 
зоне температур, оптимальный тепловой режим ее элементов и защиту 
от перегрева. Отличительной особенностью ИДТ по сравнению с тради- 
ционными термодатчиками (термисторами, термопарами и др.) является 
сравнительная простота их использования - они не требуют линеариза- 
ции и компенсации холодного спая [8-9]. В настоящее время в электрон- 
ной аппаратуре широко используется ИДТ, представляющие собой ин- 
тегральную микросхему LM135-LM335 фирмы National Semiconductor 
[10]. Однако данная микросхема не всегда обеспечивает требуемую точ- 
ность измерения температуры и может иметь лишь ограниченное при- 
менение в аппаратуре специального назначения. 

Данная работа посвящена разработке интегрального однокри- 
стального полупроводникового датчика температуры, на выходе кото- 
рого вырабатывается напряжение пропорциональное абсолютной тем- 
пературе и поисков путей повышения точности измерения температуры. 

Конструкция интегрального датчика температуры. В данной 
работе рассматривается двухвыводная интегральная схема (ИС), на вы- 
ходе которой вырабатывается напряжение пропорциональное абсолют- 
ной температуре в К. 

Рабочий диапазон данного датчика температуры от —60 до 
+125 °С, с допуском + 2 °С (+ 20 MB по выходу). 


Принцип работы датчика основан на фундаментальном свойстве 
кремниевых транзисторов, определяемых формулой Эберса-Молла [11]: 
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TU s, 
kT 
Г, = Го 8» е , (1) 


где J, — ток эмиттера, Jọ — ток насыщения; [/, — напряжение «база- 


эмиттер» транзистора; q - заряд электрона; К — постоянная Больцмана; 
S, — площадь эмиттера транзистора; Т - температура. 

Используя формулу (1) выразим разность напряжений «база-эмит- 
тер» ( AU 6» =Ибы-И 62 ) двух транзисторов: 


ЕТ 
aus" Та За (2) 
а 5 52 


Данная зависимость линейна и используется в работе электрической 
схемы датчика температуры (рисунок 1) для транзисторов УТ16 и VTI7. 


RI 7 
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Рис. 1. Принципиальная электрическая схема 
интегрального датчика температуры 


Подставим значения постоянных переменных в уравнение (2), тогда 
для температур 0 °С и25 °С получим следующие теоретические значения 
AU 65 : 54,64 мВ при 0 °С (273,2 К); 59,64 мВ при 25 °С (298,2 К). 

Выходное напряжение датчика определяется выражением: 
RitRotRs 

Ко | 

Для соотношения площадей эмиттеров 1:10 температурный коэф- 

фициент AU бъ составляет 0,2 мВ/К. 


О вых - AU 65 (3) 


Коэффициент обратной связи может быть рассчитан по следую- 
щей формуле: 


к =——^?___—0,02. (4) 


ос 
К+ Ко + R3 


Исходя из формулы (4), имеем Kc, равный 0,02. 


ос? 
Чтобы вычислить температурный коэффициент датчика использу- 
ется формула: 
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1 
AU вий рода (5) 
АТ K 


oc 
Таким образом, выходное напряжение датчика, определяемое BbI- 
ражением (3), имеет температурный коэффициент равный 10 мВ/К. 
Экспериментальные результаты. На рис. 2 и 3 приведена теоре- 


тическая [/,,.(T) и экспериментальные зависимости для шести групп 
схем, состоящих из двадцати единиц каждая. Замеры выходных напряже- 


ний были проведены в четырех точках температуры: 298 К; 333 К; 
353 К;373 К. 


Зависимость Џвых(Т) 


3800 


— — 1ая группа : : | 
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Температура (Т.К 

Рис. 2. Зависимости выходных напряжений для групп схем 

в сравнении с теоретической зависимостью 
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Рис. 3. Значения экспериментальных и теоретической зависимостей 


в точке 333.2 К 


332.6 
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Реальные измерения выходного напряжения датчика температуры 
отличаются от теоретического значения. Но не выходят за пределы до- 
пустимой нормы (+10 MB). 

В таблице | приведены температурные коэффициенты (0.) для по- 
лученных замеров выходного напряжения. 


Таблица 1 
Температурные коэффициенты замеренных выходных напряжений 
Номер О вых (мВ), а 
группы | 25°С 60 °С 80 °С 90 °С d 
1 2979 3331 (10,1) 3532 (10,08) 3731 (10,04) 10,08 
2 2979 3329 (10,04) 3531 (10,1) 3730 (10,04) 10,07 
3 2984 3335 (10,07) 3537 (10,13) 3736 (10,04) 10,08 
4 2985 3338 (10,13) 3540 (10,13) 3740 (10,09) 10,12 
5 2975 3328 (10,13) 3529 (10,08) 3728 (10,04) 10,09 
6 2975 3329 (10,16) 3530 (10,08) 3730 (10,09) 10,12 


Температурные коэффициенты вычислялись по следующей формуле: 

a = (U wx (i) Оъых (1-1) )/ AT. (6) 

Например, температурный коэффициент для первой группы схем, 

при 60 °С, рассчитывается, используя формулу (6) следующим образом: 

a = (3331-2979)/(60-25) =10,1. 

Для повышения точности измерений можно скорректировать зна- 
чение К 

Корректировка коэффициента достигается путем изменения раз- 


мера резисторов, составляющих делитель напряжения в схеме, схема- 
тично показанной на рис. 4. 


Рис. 4. Схематичное изображение делителя напряжения 
интегрального датчика температуры 


Все резисторы, составляющие делитель, имеют один номинал со- 
противления равный К: 
l-A 
R=p 
Б+А 


+13), (7) 
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где b — ширина резистора; А — погрешность; р — удельное сопротивле- 


ние; [ — длина резистора. 
Суммарное сопротивление делителя определяется формулой: 


Р N 
n2 і=1 


где nı - количество резисторов подключенных последовательно; по — 
количество резисторов подключенных параллельно, № - суммарное ко- 
личество резисторов в делителе. 

С помощью формул (7) и (8) определяется коэффициент обратной 
связи: 


N 1 
Ур 2n+— 
i-l M 
= (9) 
К/по уп 
Реальные измерения выходного напряжения датчика температуры 
отличаются от теоретических значений. Проверим козффициент обрат- 
ной связи для реальных значений: 
АО [АО j* 2982/59,84 - 49,83. (10) 


бэ(реал 
Для данного К резисторы делителя напряжения интегрального 


датчика температуры корректируются. В результате резисторы имеют 
следующие номиналы и габариты: 
1. К, =16480 Ом (8 резисторов соединенных последовательно с ra- 


баритами 1.=126 мкм и b =10 мкм). 
2. R, = 686,66 Ом (3 резистора соединенных параллельно с габа- 


ритами 1.=126 мкми b =10 мкм). 
3. В. =17260 Ом (4 резистора с габаритами 1.=126 мкм и b=10 


мкм и 4 резистора с габаритами L =139 мкм и b =10 мкм соединенных 
последовательно). 

Таким образом, корректируя ширину и длину резисторов можно 
добиться оптимального значения К „при котором погрешность измере- 


ния выходного напряжения будет минимальна. 

Выводы. 

Разработана и проанализирована схема электрическая принципи- 
альная, лежащая в основе интегрального датчика температуры, на вы- 
ходе которого вырабатывается напряжение пропорциональное абсолют- 
ной температуре. Получены графики выходного напряжения для не- 
скольких групп микросхем, на основе которых проведено сравнения ре- 
альных значений с теоретическими. На основе экспериментальных ис- 
следований скорректированы конструкция и топология делителя напря- 
жения, что позволило повысить точность измерения температуры. 
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INTEGRATED TEMPERATURE SENSOR 


This article is devoted to the development of an integral single-chip semi- 
conductor temperature sensor, the output of which produces a voltage propor- 
tional to the absolute temperature and the search for ways to improve the ac- 
curacy of temperature measurement. Experimental graphs of the output voltage 
were obtained, on the basis of which the design and topology of the voltage 
divider were corrected, which allowed increasing the accuracy of temperature 
measurement. 


Keywords: Temperature sensor, integrated circuit, temperature coefficient, 
measurement accuracy. 
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УДК 541.49+544.342+546.492'131'141 
А.В. Шуваев 


СПЕКТРОФОТОМЕТРИЧЕСКОЕ ИССЛЕДОВАНИЕ 
РАВНОВЕСИЯ КОПРОПОРЦИОНИРОВАНИЯ 

B ХЛОРИДНОБРОМИДНОЙ СИСТЕМЕ РТУТИ (II) 
В ВОДНОМ РАСТВОРЕ 


Спектрофотометрическим методом изучено равновесие про- 
цесса Не С1, + HgBr, -2HgCIBr в водном растворе при ионной силе 


I 20,001 (НОО) u t 20+ ЕС. Результаты измерений анализирова- 


лись методами математической статистики; статистически зна- 
чимых различий между результатами не было обнаружено. Полу- 
чено значение константы равновесия 15К = 0,687 + 0,061, которое в 
пределах точности измерений совпадает с другими надежными ли- 
тературными данными. Отмечено, что оптические характери- 
стики хлориднобромидной системы ртути (11) оказались менее бла- 
гоприятными для измерения константы равновесия по сравнению с 
другими галогенидно-псевдогалогенидными системами ртути (11). 


Ключевые слова: диацидокомплексы ртути (I1), константы рав- 
новесий копропорционирования, планирование эксперимента, спек- 
трофотометрия, методы математической статистики. 


Введение. Процессы образования смешанных диацидокомплексов 
ртути (П) из соответствующих однородных: 

HgX, + НУ, -2HgXY (1) 
в различных растворителях представляют интерес для многих исследо- 
вателей [1-6]. На протяжении последних пятидесяти лет получено до- 
статочно большое количество данных о константах равновесий (1) [7]. 
Однако сопоставление результатов разных авторов для одной и той же си- 
стемы нередко обнаруживает их расхождения, существенно превышаю- 
щие уровень ошибок воспроизводимости. Критический анализ [8] приме- 
няемых методик исследования разными авторами показал, что в них BO3- 
можно содержание источников систематических погрешностей из-за 
необоснованности ряда предположений, которые даже при небольших 
погрешностях в первичных экспериментальных данных могут приво- 
дить к значительным искажениям конечных результатов. 

Информацию о характеристиках процесса (1) можно получить в 
рамках надежно себя зарекомендовавшего метода соответственных рас- 
творов [9] в виде линейной комбинации констант равновесий двухступенча- 
того замещения лигандов: 


HgX, +2У =X + НХУ +Y -2X -HgY,. (2) 
Впоследствии в рамках этого метода нами была изучена хлорид- 


нороданидная система ртути (П) в водном растворе [10]. Если в водных 
растворах изучение равновесий (2) не связано с проявлением каких-либо 
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сложностей, то при переходе к неводньгм растворителям они появляются в 
виде ограниченной растворимости и неполной ионной диссоциации со- 
лей, образования ионных пар, а также ион-дипольных взаимодействий. 
По этой причине более предпочтительны процессы лигандного копро- 
порционирования (1) в системах нейтральных, высокоустойчивых и хо- 
рошо растворимых в органических растворителях галогенидах и псев- 
догалогенидах ртути (II) (X, Y = А, Br, I, CN, $ СМ). 

Для изучения равновесий (1) бъл разработан спектрофотометриче- 
ский метод [8], основу которого составляет планирование эксперимента 
в виде поиска оптимальных условий измерения, обеспечивающих макси- 
мальную точность определения величины константы равновесия про- 
цесса (1) при заданном исходном уровне ошибок в первичных экспери- 
ментальных данных. Существенным отличием этого метода от аналогич- 
ных, применявшихся другими исследователями, является то, что в нем 
изначально отсутствуют какие-либо необоснованные постулаты и пред- 
положения. В дальнейшем для некоторых обозначенных выше систем ди- 
ацидокомплексов ртути (П) были получены данные в ряде растворителей: 
водном растворе [8, 11], бензоле, этаноле [8], диоксане [12] и газовой фазе 
[8]. Практическое использование применяемой нами методики исследо- 
вания [8] сопряжено с достаточно большим количеством вспомогатель- 
ных расчетов, необходимых для поиска оптимального плана экспери- 
мента. Впоследствии нами были найдены способы [13], позволяющие су- 
щественно облегчить и ускорить процедуру планирования эксперимента 
без потери точности определения констант равновесий (1). Поскольку 
еще не для всех галогенидно-псевдогалогенидных диацидосистем ртути 
(П) получены данные, было принято решение продолжить их исследова- 
ния. В данной работе приводятся результаты изучения хлориднобромид- 
ной системы в водном растворе. Вначале коротко изложим основы про- 
цедуры планирования эксперимента, подробно описанные в работе [13]. 

Методика исследования. Первоначально в широком интервале 
волновых чисел измеряются спектры растворов: а) только HgX,; б) 
только Ноу, ; B) НеХ, и HgY,B соотношении 1:1; г) Нех, и Ноу, при 
соотношении концентраций К = 10. Затем строятся графические зависи- 
мости т и а OT волнового числа V: 

2-=1- (£5, Ps): 


т = a = Ey / Ep2 (3) 
£59 T£ 


где 2,6, 5, — молярные коэффициенты зкстинкций Нох,, Ноу, и 


202 702, 71 
раствора при К =1. Выбираются оптимальные оптические характери- 
стики для исследований на основе поиска максимальных значений | т, | 
при условии, что величина а также по возможности была максимальной. 
Затем с помощью вспомогательной функции y решается задача поиска оп- 


тимальных концентрационных условий. Особенность такой функции за- 
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ключается в следующем. С одной стороны, в явном виде ее можно пред- 
ставить как функцию двух переменных: величины К - константы рав- 
новесия (1) и К — соотношения концентраций исходных реагентов. По- 
этому ее числовые значения можно предварительно рассчитать и пред- 
ставить их в виде набора графических зависимостей y — y(R) для ряда 


фиксированных К. С другой стороны, функция у достаточно легко 
определяется на основе спектрофотометрических измерений в виде 


у = Y(£1,£n,£45,£,,, А). Таким образом, используя функцию у, можно до- 
статочно легко и быстро, избегая многочисленных промежуточных рас- 
четов, осуществить переход от непосредственно полученных экспери- 
ментальных данных к конечному результату — величине К. При этом ос- 
новная задача состоит в том, чтобы минимизировать трансформацию 
уровня погрешностей первичных экспериментальных данных на конеч- 
ный результат. В работе [13] показано, что поиск оптимальных концен- 
трационных условий можно осуществить с помощью функции 

РОК) |= а,5+Х/|т|): 5, (4) 
где значения X и = табулированы для разных А, К. Для построения 
плана эксперимента необходима предварительная оценка величины К, 
ее можно осуществить, вычисляя значение Y, ,,, азатем и К, по резуль- 


татам оптических измерений при выбранных волновых числах. Область 
значений К, при которых функция f(R) достигает минимума, является 
наилучшей для дальнейших измерений. 

Восемь из десяти систем галогенидных и псевдогалогенидных ди- 
ацидокомплексов ртути (П) в водном растворе были изучены нами ранее 
18, 11]. В большинстве из них, за исключением цианидноиодидной и циа- 
ниднороданидной, оптические характеристики оказались благоприятными 
для исследования в рамках используемой нами методики. Иодиднороданид- 
ная система характеризуется крайне низкими значениями |т, |< 0,05 и 


Из нат, / € ig sc), «0,7, uro исключает возможность ее исследования. 


Осталась неизученной хлориднобромидная система, в данной работе 
приводятся результаты измерений в этой системе. 
Экспериментальные результаты и их обсуждение. В работе ис- 
пользовались дистиллированная вода; хлорная кислота марки "х.ч."; 
НоСЬ и HgBr, , синтезированные согласно методике [14]. Добавление в 


растворы хлорной кислоты с концентрацией 10? моль/л осуществля- 
лось, чтобы предотвратить гидролиз дигалогенидов ртути (11). Сначала 
из реактивов в калиброванной посуде готовили стандартные (исходные) 
растворы, а затем рабочие растворы. Максимальная концентрационная 
ошибка в приготовлении исследуемых растворов не превышала 0,3 %. Ре- 
гистрацию оптической плотности равновесных растворов осуществляли 
на спектрофотометре "Specord" при t = 20 + 1 °С. Во всех исследуемых рас- 
творах равновесие достигалось за время, не превышающее 1 мин. 
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Данные рисунка 1 использовали для выбора оптимальных волно- 
вых чисел. Модуль величины т; имеет максимальные значения при у = 


47 000 и 42 000 см“. Они небольшие - при 47 000 см“! всего около 0,16, 
а Ен.с, И &ньвг, ЗДесь приблизительно равны. При у = 42 000 см”! MaK- 


симум |m | еще меньше = 0,10, но более выгодно соотношение 
новь, / нас, - Для дальнейших измерений мы выбрали эти два значения 


волновых чисел. 


50 46 42 38 у тыс. см" 
Рис. 1. Зависимость m, и 10 ,,в,/&н,с, ОТ ВОЛНОВОГО числа: 


1- т = [Ж (Юс, F Dy, (Эм, + Dus, ) = f(v); 


2 E lg € ops / É несь = f(v) 


При выбранных волновых числах по результатам оптических H3- 
мерений при К- 10 вычисляли значение функции Ур, оно оказа- 


лось равным 0,542. Затем, используя графические зависимости у от 
је К при К=10 [13], определили предварительное значение KOH- 
станты равновесия (1) lg К =0,6( Hg —С1— Вг). На основании таблич- 
ных данных работы [13] рассчитали как составные части г, X, так и 
сами значения функции f(R) при двух волновых числах, эти данные 
приведены в таблице 1. 
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Таблица 1 
Расчетные значения функций z(R), X (К), f (К) 
при разных R для двух волновых чисел 
Параметр R 
у 5 10 Е: 26 30 40 
120,42; Ig È 2-06 É -40 
z(R) 9,0 6,0 5,2 4,9 4,6 4,4 
у=47 000 см; m, = 0,16; ва=0,02 
X(R) 2,18 | 335 | 453 | 5,79 | 824 | 10,78 
(15+ X/|m]|) | 15,12 | 22,44 | 29,81 | 37,68 | 53,01 | 68,88 
f(R) 136 | 135 155 185 | 244 | 303 
у = 42 000 см"; m, =-0,10; lg a = 1,05 
X(R) 1,4 1,7 1,9 2,1 25 3,0 
(15+ X/|m]|) | 15,5 | 185 | 20,5 | 22,5 | 27,5 | 31,5 
f(R) 139 | 111 107 110 126 139 


Можно видеть, что минимум функции f(R) достигается при 
К-5и 10 для v = 47 000 см"; Ё=10 и 15 jum v = 42 000 см". Для 
дальнейших измерений были выбраны именно эти значения. Следует 
также отметить, что абсолютные значения функции f (R) в области ми- 


нимума оказались выше в 2-3 раза, чем для ранее исследованных других 
систем диацидокомплексов ртути (П), для которых оптические характе- 
ристики оказались более благоприятными, чем в хлориднобромидной 
системе. В таблице 2 содержатся экспериментальные данные для одной 
отдельной серии измерений. 


Таблица 2 
Результаты измерений в серии 3 


Оптическая плотность 


Хо Состав раствора в кювете 0,1 см Количество независи- 
р-ра (ммоль/л) при у =47 000 см“ мых измерений 
1 Със, = 1,511 + 0,002 0,4650 + 0,0010 3 
2 Сив, = 1,512 + 0,002 0,4933 = 0,0017 3 
Сыс, = 1,511 + 0,002 
3 1,1100 + 0,0029 3 
С в, = 1,512 + 0,002 
Сива, = 1,511 + 0,002 
4 0,53402 = 0,0007 3 
Сев, = 0,1501 + 0,0002 


у =0, 575 + 0,026. В качестве ошибок – среднеквадратичные отклонения 


одного измерения 
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Всего было сделано 12 серий, в которых проводили от трех до ше- 
сти независимъх измерений для каждого раствора серии, для разнъх К, 
двух значений Сүн, И двух волновых чисел. Полная информация о по- 


лученных результатах содержится в таблице 3. 


Таблица 3 
Количественная характеристика результатов 
к NU Bon R , i 10 K ME ree ши 
1 47 1,5 5 | 0,596 23 3,77 | 1,37 | +1,09 11881 
2 -"- -"- 5 | 0,529 29 6,04 | 2,40 | -1,18 13924 
3 -" -" 10 | 0,575 26 5,22 | 1,30 | -0,36 1296 
4 -"- -"- 10 | 0,584 56 5,73 | 3,50 | -0,87 7569 
5 -"- 3,0 5 | 0,615 21 4,98 | 2,10 | -0,12 144 
6 -"- -"- 10 | 0,602 31 6,90 | 2,13 | -2,04 41616 
7 42 1,5 10 | 0,541 29 3,63 | 1,18 | +1,23 15129 
8 -"- -"- 10 | 0,565 23 4,70 | 1,16 | +0,16 256 
9 -"- -"- 15 | 0,549 36 4,66 | 1,34 | +0,20 400 
10 -"- -"- 15 | 0,513 28 3,87 | 0,93 | +0,99 9801 
11 -" 3,0 10 | 0,536 16 3,43 | 0,52 | +1,43 20449 
12 -"- -"- 15 | 0,565 30 5,43 | 1,50 | -0,57 3249 
5 Fe 5 к - средние квадратичные погрешности средних арифметических 
(F и К ) вотдельных сериях измерений 


На основе методов математической статистики [15] осуществ- 
ляли анализ полученных результатов. Для этого данные табл. 3 груп- 
пировали по признакам одинаковых Се, или V, или К. Для каждой 


группы производили расчет средних значений и среднеквадратичных 
ошибок (5 x) (табл. 4). Сопоставление значений среднеквадратичных 


ошибок осуществляли по Ё-критерию, а средних результатов - по /- 
критерию [15]. Можно видеть, что между группами с одинаковыми 


признаками — С у, К – не обнаруживается статистически значи- 


HgCl, ? 


MBIX различий. 
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Таблица 4 
Распределение величин К по группам с одинаковьгми признаками 


Са, =Ъ5 ммоль/л Са, =3,0 ммоль/л 


3,77; 6,04; 5,22; 5,73; EITE 
3.63: 4.70: 4.66: 3,87 4,98; 6,90; 3,43; 5,43 


К =4,70; S, =0,91 K =5,18; S, =1,42 
Е=2,43; ЕР-0,95, f, =3, f, =7)=4,35 
t=1,29; Р-0,95, f =10)=2,23 
v = 47 000 см“! у = 42 000 см“! 


3,77; 6,04; 5,22; 3,63; 4,70; 4,66; 
5,73; 4,98; 6,90 3,87; 3,43; 5,43 


К =5,44; S, 21,06 K' =4,29; S, 20,77 
Е=1,90; Р(Р=0,95, f, =5, f, =5)=5,05 
#=2,15; (Р=0,95, f =10)=2,23 
R-5 R -10 R-15 


. 604: 5,22; 5,73; 6,90; .4 97. 
3,77; 6,04; 4,89 3,63; 4,70; 3,43 4,66; 3,87; 5,43 


К 24,93; S, =1,14 К” = 4,94; 5, 21,31 K^ =4,65; $, 20,78 
F =S; / S? 2142; (Р -0,95, f, =5, f, =2)=19,30 
1= 0,01; (P 20,95, f =7)= 2,36 
Е = 52152 = 2,82; Е(Р = 0,95, f, =5, f, =2)=19,30 
1=0,41; t(P = 0,95, f =7)= 2,36 
Е= 52152 =2,14; Е(Р = 0,95, f, -2, f, =2) =19,00 
1= 0,35; t(P = 0,95, f =4)= 2,78 


Для всех серий независимых измерений результаты оказались 
приблизительно равноточны, поэтому в качестве наиболее вероятной 
оценки величины К приняли среднее арифметическое, а дисперсию 
среднего арифметического вычислили по Стьюденту, Р = 0,95: 


К +S- = 4,86+0,68. Соответственно логарифм константы равновесия 


(1) (X 2 CLY = Br) Ig K = 0,687 + 0,061. Этот результат находится B XO- 
рошем соответствии с надежными литературными данными 0,66 + 0,04 
[16] и 0,58 + 0,06 [17]. 


Несмотря на не совсем благоприятные оптические характеристики 
хлориднобромидной системы, за счет увеличения статистических дан- 
ных - числа серий независимых измерений, удалось получить результат 
с точностью, сопоставимой для других диацидосистем ртути (П) [8], ха- 


рактеризующихся лучшими для исследования оптическими характери- 
стиками. 
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SPECTROPHOTOMETRIC STUDY OF COPROPORTION 
EQUILIBRIUM IN CHLORIDE-BROMIDE SYSTEM 
OF MERCURY (II) IN AQUEOUS SOLUTION 


The equilibrium of the process HgCl, + HgBr, = 2HgCIBr in an aqueous 
solution at an ionic strength of I = 0,001 (НСІО,) and t = 20 + РС was studied 


spectrophotometrically. The results of the measurements were analyzed by 
methods of mathematical statistics, no statistically significant differences be- 
tween the results were found. The value of the equilibrium constant 
Тек = 0,687 + 0,061 is obtained; the result coincides with other reliable liter- 


ature data within the limits of measurement accuracy. The optical characteris- 
tics of the chloride-bromide system of mercury (II) is noted to less favorable for 
measuring the equilibrium constant compared to other mercury (II) hal- 
idepseudo-halogenide systems. 


Keywords: diacidocomplexes of mercury (II), equilibrium constants of the 
coproportion, experiment planning, spectrophotometry, mathematical statistics 
methods. 
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УДК 621.382 
С.А. Корнев, В.В. Андреев 


СИНТЕЗ И МОДЕЛИРОВАНИЕ МИКРОСХЕМЫ 
ДВУХКАНАЛЬНОГО ФОРМИРОВАТЕЛЯ 
ПРЯМОУГОЛЬНЫХ ИМПУЛЬСОВ 


Данная работа посвящена рассмотрению особенностей проекти- 
рования и моделирования цифровых микросхем, с использованием 
языков описания аппаратуры (HDL) и программных пакетов фирмы 
Cadence: Genus и Virtuoso. Показано, что использованный маршрут 
проектирования позволяет оптимизировать энергопотребление и 
площадь получаемой микросхемы, а также ускорить процесс разра- 
ботки, значительно снизив время, затрачиваемое на тестирование и 
верификацию устройства. 


Ключевые слова: Cadence Virtuoso, Cadence Genus, НОГ, Verilog, 
ModelSim, KMOII микросхема, цифровая схемотехника, языки onuca- 
ния аппаратуры, моделирование. 


Введение 

В современной электронике, для решения различных задач, все 
чаще применяются заказные интегральные микросхемы (ASIC), изго- 
товленные по КМОП технологии. К качеству и надежности таких мик- 
росхем и их элементов предъявляются повышенные требования [1-5]. 
Одним из путей повышения надежности и быстродействия микросхем 
является размещение на одном кристалле различных функциональных 
блоков [6-7]. Однако, в этом случае, может возрастать вероятность оши- 
бок, допускаемых при разработке изделия, а процесс верификации про- 
екта приобретает дополнительную значимость. 

В настоящее время для разработки и верификации цифровых схем, 
широко используются языки описания аппаратуры VHDL и Verilog [8, 
9]. Они позволяют произвести описание схемы на поведенческом или 
логическом (схемном) уровне, а также провести её тестирование. 

В данной работе рассмотрен процесс синтеза схемы двухканаль- 
ного формирователя прямоугольных импульсов, которая используется 
для управления (перемагничивания) фазовращателей в фазированных 
антенных решетках (ФАР). Одним из основных критериев оптимизации 
в процессе синтеза, является минимизация размеров данной схемы при 
её реализации на КМОП технологии с топологическими нормами 0,6 
мкм, предоставляемой компанией X-Fab. Сравнение размеров получен- 
ной в ходе синтеза схемы проводится с исходной реализацией на той же 
технологии, созданной без использования языков описания аппаратуры. 

Верификация и моделирование схемы на уровне регистровых пе- 
редач (RTL) проводились с использование программы Моде $ ит, co- 
зданной Mentor Graphics [10]. Для синтеза схемы и привязки её к вы- 
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бранной КМОП технологии использовалось ПО Cadence Genus, схемо- 
техническое и топологическое проектирование, а также моделирование 
схемы с использованием SPICE моделей BSIM3v3 производилось в Са- 
dence Virtuoso. 

Анализ исходной схемы формирователя прямоугольных HM- 
пульсов 

На рис. 1, представлена функциональная схема оригинального 
устройства, используемого для формирования прямоугольных управля- 
ющих импульсов. Данная схема является составной частью драйвера 
управления индуктивной нагрузкой и необходима для управления сило- 
выми К-МОП транзисторами, используемыми для формирования втека- 
ющих токов. В схеме присутствует два независимых канала, с двумя вы- 
ходами формирования прямоугольных импульсов в каждом. 


Канал 1 


Ае 


RST 


Рис. 1. Функциональная схема устройства 


К схеме предъявляются высокие требования по быстродействию: 
задержка распространения сигнала не должна превышать 50 нс, для каж- 
дого из управляющих выходов (ОПТА, OUTB, OUTC, OUTD) во всём 
температурном диапазоне работы микросхемы (-60 + 125 °С) и ёмкост- 
ной нагрузке в 1500 пФ, соответствующей ёмкости затвора используе- 
мого полевого транзистора. 

Рассмотрим функционирование устройства на примере работы от- 
дельного канала. Второй канал имеет полностью идентичную первому 


структуру блоков. 
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Блок 1 разрешает формирование прямоугольного импульса на вы- 
ходе схемы. Данный блок имеет активный низкий уровень разрешаю- 
щего сигнала, а также не позволяют произвести формирование прямо- 
угольного импульса одновременно на двух выходах одного канала. 

Блок сброса 2, позволяет сбросить все выходы драйверов отдель- 
ного канала (RST), либо двух каналов схемы одновременно (GRST). 
Сигналы сброса отдельных каналов имеют активный низкий уровень, а 
общий сброс - активный высокий уровень. 

Для формирования прямоугольных импульсов на выходе исполь- 
зуются блоки 3 и 4. При активном разрешающем сигнале, и отсутствии 
сигнала сброса, блок устанавливает высокий уровень на выходе, после 
детектирования фронта сигнала $їагї. При отсутствии сигнала, разреша- 
ющего формирование прямоугольного импульса и/или активном сиг- 
нале сброса, на выходе устанавливается низкий уровень. 

Блоки 5 и 6 информируют о состоянии выходов схемы: на выходе 
блоков будет высокий уровень, если хотя бы один из входов имеет высокий 
уровень. Также блок 5 можно перевести в высокоимпедансное состояние, 
подав низкий уровень на его вход разрешения (ЕХ АВ, EN. CD). 

Использованные при создании исходной схемы (разработанной 
без применения HDL) логические элементы и их количество представ- 
лены в таблице 1. Большое количество используемых инверторов объ- 
ясняется их использованием в детекторе фронта для формирования за- 
держки сигнала (блоки 3, 4). 


Таблица 1 
Элементы, использованные в исходной схеме 
Элемент Количество, шт. Площадь, мкм” 

И-НЕ, 4 вывода 4 1 104,00 
И-НЕ, 2 вывода 4 728,64 

ИЛИ-НЕ, 2 вывода 11 2 671,68 
HE 97 10 679,70 
Итого 116 15 184,02 


Согласно техническому заданию, микросхема должна обеспечи- 
вать время задержки распространения сигнала не более чем 50 нс, для 
выходов ОПТА, OUTB, OUTC, OUTD при ёмкостной нагрузке 1500 пФ; 
а также 50 нс, для выходов OUTAB, OUTCD, СОМ, при ёмкостной 
нагрузке 150 пФ. Значения параметров должны оставаться в заданных 
пределах во всём температурном диапазоне -60+ 125 °С. 

Синтез и верификация Verilog модели схемы 

На основе данных, полученных в ходе анализа функциональной 
схемы устройства, было разработано описание схемы на языке Verilog, 
с его использованием в процессе синтеза была получена схема, пред- 
ставленная на рис. 2. 
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outb 
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DFRROX1 №. C—outab 


Рис. 2. Схема устройства, полученная B результате синтеза 


Схема состоит из 4 D триггеров (DFRRQX 1), имеющих вход асин- 
хронного сброса с активным низким уровнем (RN), вход загрузки дан- 
ных по переднему фронту тактового сигнала (С) и информационный 
вход О подключённый к логической единице. 

Для формирования сигнала сброса О триггеров, в схеме использу- 
ются элементы ИЛИ-НЕ с 4 входами (МО412Х 1), два из которых HH- 
версны. Элементы ИЛИ-НЕ с 2 входами (NO2XI) и инвертирующий 
элемент c High Z состоянием (ITHXI) задают выходной сигнал на HH- 
формационных выходах OUTAB и OUTCD. Также для задания выход- 
ного сигнала СОМ используется элемент И-НЕ (NA2XI) с двумя вхо- 
дами, на входы которого поступают сигналы, формируемые элементами 
ИЛИ-НЕ (NO2X1). 

Для проведения верификации была создана тестовая среда, напи- 
санная на Verilog. В ходе тестирования (рис. 3) на входы разработанной 
модели подавались заранее определённые тестовые последовательности 
сигналов, при этом производилось сравнение выходных сигналов 
схемы, с ожидаемым идеальным состоянием выходов. 
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Входные сигналы 


Выходные сигналы 


Проверочные сигналы 


Рис. 3. Временная диаграмма сигналов тестируемой схемы 
в программе Моде $ ша. Входные сигналы: А — SET, выходные сигналы 
устройства: СОМ – OUTD, проверочная последовательность ожидаемых 
выходных сигналов: REF ОСТА - REF COM 


Ha puc. 3 представлена временная диаграмма сигналов TecTHpye- 
мой схемы. Входные комбинации сигналов считываются из файла — это 
позволяет не менять схему тестирования при изменении входных воз- 
действий и делает тест более универсальным. 

Анализ временных диаграмм, полученных в ходе тестирования, 
показал полное соответствие откликов разработанной и исходной схем 
на входные воздействия. 

Моделирование схемы на уровне регистровых передач не может 
дать полной информации о задержке распространения сигнала в схеме, 
так как не учитывает ёмкость нагрузки, разброс параметров технологи- 
ческого процесса и другие факторы, но в то же время позволяет прове- 
рить логику работы на наличие ошибок и исправить их без проведения 
ресурсоёмких и обладающих большей вычислительной сложностью де- 
тальных видов тестирования [11-12]. 

Следующим шагом после проведения верификации, является син- 
тез схемы на элементной базе предоставляемой выбранной технологией. 
В процессе синтеза производится оптимизация схемы по различным па- 
раметрам: энергопотреблению, занимаемой площади и задержке рас- 
пространения сигналов, также на схему накладываются временные 
ограничения [12-14]. 

Синтез схемы производился в ПО Cadence Genus, которая позво- 
ляет перевести Verilog описания схемы в список цепей (netlist) с цифро- 
выми библиотечными логическими элементами, предоставляемыми ис- 
пользуемой технологией. Использованные при синтезе элементы и их 
количество указаны в таблице 2. 
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Таблица 2 
Элементы использованные при синтезе схемы 
Элемент Количество, шт. | Площадь, мкм? 
D триггер 4 3 886,08 
Инвертор c High 7, состоянием 2 607,20 
ИЛИ-НЕ, 4 вывода (2 инверсные) 4 2 185,92 
ИЛИ-НЕ, 2 вывода 2 485,76 
И-НЕ, 2 вывода 1 182,16 
Итого 13 7 347,12 


После проведения синтеза, полученная схема была экспортиро- 
вана в САПР Cadence Virtuoso, для проведения более детального тести- 
рования и верификации, по сравнению с ModelSim, а также измерения 
задержки распространения сигнала. Cadence Virtuoso позволяет учесть 
при тестировании влияние температуры на параметры транзисторов, 
разброс параметров технологического процесса, а также другие пара- 
метры полупроводниковых устройств [15-16]. На этом этапе, к получен- 
ному цифровому ядру схемы были добавлены выходные драйверы и 
триггеры Шмитта на входы схемы. 

На Рис. 4 представлена временная диаграмма сигналов моделиру- 
емой схемы, полученная B Cadence Virtuoso. На выходы схемы подклю- 
чена ёмкостная нагрузка, согласно техническому заданию, это 1500 пФ 
для выходов ОПТА, OUTB, OUTC, OUTD и 150 пФ для выходов 
OUTAB, OUTCD, СОМ. Как видно из временной диаграммы, представ- 
ленной на рисунке, поведение схемы полностью идентично разработан- 
ной ранее Verilog модели. 
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Рис. 4. Временная диаграмма сигналов тестируемой схемы B Cadence Virtuoso. 
Входные сигналы: SET — СОЯ, выходные сигналы устройства: ОСТА - COM 
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Моделирование работы схемы проводилось при значениях темпе- 
ратуры, соответствующих границам рабочего диапазона (-60 и 125 °С) 
и температуре равной 25 °С. Также учитывалось влияние технологиче- 
ского процесса на разброс параметров элементов. Для этого было про- 
ведено Монте-Карло моделирование, где поведение схемы для каждой 
из представленных выше температур было промоделировано в 200 точ- 
ках. В качестве примера на рис. 5 и рис. 6 показаны гистограммы рас- 
пределений значений максимальной задержки распространения сигна- 
лов в полученной схеме при ёмкостной нагрузке в 150 и 1500 пФ соот- 
ветственно. На рисунках использованы следующие обозначения: р — 


среднее значение, min — минимальное значение, тах — максимальное 
значение, о — среднеквадратичное отклонение. 
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Рис. 5. Гистограмма распределения максимального значения задержки 
распространения сигнала при нагрузке 150 пФ, напряжении питания 5 В 
и температуре 25 “С. Параметры распределения: 

u — 13,09 нс, min — 12,40 nc, max – 13,96 нс, с — 0,29 пс 


em . С 


Рис. 6. Гистограмма распределения максимального значения задержки 
распространения сигнала при нагрузке 1500 пФ, напряжении питания 5 В 
и температуре 25 °С. Параметры распределения: 

u — 36,66 нс, тіп — 34,30 нс, тах – 39,43 nc, с - 0,92 пс 
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Данные о величине максимальной задержки распространения сиг- 
нала, полученные в ходе Монте-Карло моделирования при разных зна- 
чениях температуры и ёмкостной нагрузке в 150 пФ и 1500 пФ, пред- 
ставлены в таблицах Зи 4 соответственно. 


Таблица 3 
Параметры распределения значения максимальной задержки 


распространения сигнала (нагрузка 150 пФ, напряжение питания 5 В) 


Температура, Тип Мин., Макс., Ср. Ср. квадр. 
ме задержки нс нс знач., нс откл., пс 
Тон, 7,49 8,54 8,04 0,19 
Е Трн 9,64 10,69 10,11 0,19 
Bs Тон, 8,43 9,77 9,15 0,22 
Трн 12,40 13,96 13,09 0,29 
ios Тон, 9,56 11,22 10,47 0,26 
Тон 15,51 17,64 16,47 0,40 
Таблица 4 


Параметры распределения значения максимальной задержки 
распространения сигнала (нагрузка 1500 пФ, напряжение питания 5 В) 


Температура, Тип Мин., Макс., | Ср. знач., Ср. квадр. 
ще задержки нс нс нс откл., пс 

Tur, 19,58 21,43 20,47 0,33 

b Тон 26,50 29,99 26,10 0,63 

9s Тн, 21,86 24,36 23,14 0,42 

Тон 34,30 39,43 36,66 0,92 

е Тн. 24,87 | 28,16 26,61 0,55 

Тн 42,85 49,70 46,00 1,22 


На основе результатов моделирования, представленных в табли- 
цах выше, было определено, что задержка распространения сигнала не 
превышает предельные значения, указанные в технических условиях (50 
нс). В среднем изменение задержки распространения сигнала во всём 
температурном диапазоне работы составляет 3,5 нс для выходов с 
нагрузкой 150 пФ и 10 нс для нагрузки в 1500 пФ. 

Из полученных данных следует, что величина задержки распро- 
странения сигнала прямо пропорциональна температуре, а также что ве- 
личина задержки распространения сигнала при переходе выхода с низ- 
кого логического уровня в высокий Т, у, превышает задержку распро- 


странения сигнала при переходе выходного сигнала с высокого логиче- 
ского уровня в низкий T, . 
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Выводы 

В данной статье была разработана микросхема, реализующая циф- 
ровой формирователь прямоугольных импульсов, используемый как со- 
ставная часть драйвера управления индуктивной нагрузкой. Разработан- 
ная микросхема управляет силовыми М-МОП транзисторами, предна- 
значенными для формирования втекающих токов. 

Применение методики проектирования цифровых микросхем, с 
использованием языка описания аппаратуры Verilog, позволило опти- 
мизировать параметры схемы, в частности, уменьшить занимаемую 
площадь на кристалле и количество использованных элементов. Коли- 
чество использованных в схеме элементов, было уменьшено в 8,9 раза, 
а занимаемая площадь в 2 раза. 

Использование языка Verilog при проектировании позволило про- 
изводить полное тестирование поведения схемы, начиная с ранних эта- 
пов разработки, и значительно снизить время, затрачиваемое на тести- 
рование и верификацию устройства. Это было достигнуто путём пере- 
носа тестов с использованием SPICE моделей на этап схемотехниче- 
ского проектирования B Cadence Virtuoso, когда поведение схемы уже 
полностью отлажено. Моделирование времени задержки распростране- 
ния сигнала на схемотехническом этапе проектирования устройства по- 
казало, что полученная схема удовлетворяет ограничениям, указанным 
в техническом задании. 
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ПОРТАТИВНОЕ УСТРОЙСТВО ДЕТЕКТИРОВАНИЯ 
ГАММА-ИЗЛУЧЕНИЯ 


В данной статье представлено разработанное устройство де- 
тектирования гамма-излучения, используемое для мониторинга ра- 
диационного фона окружающей среды и обучения студентов. Рас- 
сматривается процесс проектирования устройства, выбор основ- 
ных компонентов устройства, подбор вспомогательных электрон- 
ных средств, выбор системной программной части и разработка 
прикладной программной части устройства. Выполнен краткий об- 
зор возможностей разработанного устройства. 


Ключевые слова: проектирование, разработка, спектрометрия, 
радиационная безопасность. 


На кафедре «Конструирование и производство РЭС» (МИРЭА - 
Российский технологический университет (РТУ МИРЭА)) совместно с 
предприятием Общество с ограниченной ответственностью научно-ис- 
следовательское производственное предприятие «Грин Стар Инстру- 
ментс» (ООО НИПП «ГСИ») разработано устройство детектирования 
гамма-излучения. Устройство регистрирует ионизирующее излучение и 
посредством USB порта передает обработанную информацию на персо- 
нальный компьютер или ноутбук для получения амплитудных спектров 
ионизирующего излучения. Далее полученная информация обрабатыва- 
ется специализированным программным обеспечением. 

Устройство является носимым, питание осуществляется через 
USB порт. За счет малого потребления энергии аккумуляторной батареи 
ноутбука хватает для длительной работы без дополнительной подза- 
рядки. В устройстве реализован активный, управляемый восстановитель 
базовой линии. Устройство позволяет отказаться от использования ка- 
ких-либо дополнительных средств при построении измерительных ком- 
плексов для регистрации ионизирующего гамма-излучения. 

Необходимость разработки такого устройства была продикто- 
вана тем, что имеющиеся на российском рынке изделия при сравни- 
тельно небольшой разнице в цене, сильно уступают по параметрам 
зарубежным аналогам. Стоимость устройства детектирования произ- 
водства фирмы ORTEC, состоящего из многоканального анализатора 
DIGIBASE и сцинтилляционного блока детектирования может дохо- 
дить до 9000 $, а стоимость аналога производства российской фирмы 
«ЭкоСфера» MKC-AT6101B составляет 7000 $. Устройство МКС- 
АТ6101В имеет повышенные массогабаритные характеристики и яв- 
ляется сложным в эксплуатации. 

По сравнению с существующими отечественными аналогами, 
разработанное устройство превосходит их по следующим основным 
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параметрам, а именно: знергетическое разрешение на энергии 662 
кэВ порядка 7 %, максимальная загрузка спектрометра 10° импульс/с, 
интегральная нелинейность - 0,05 %, вес - не более | кг. По сравнение 
с зарубежными аналогами устройство значительно превосходит их по 
критерию цена - качество. Стоимость разработанного устройства со- 
ставляет 3000 $. 

Наша жизнь тесно связана с использованием ионизирующих ис- 
точников, использование радиоактивного топлива на атомных злектро- 
станциях, использование в медицинских целях (рентгеновское обследо- 
вание, лучевая терапия). Вследствие этого появляется большое количе- 
ство отходов, хранение и утилизация которых иногда выходит из-под 
контроля. Завод полиметаллов, который сейчас принадлежит госкорпо- 
рации «Росатом», полвека назад складировал отходы урана на берегу 
Москвы-реки. На отрезке в полкилометра длиной опасные вещества 
просто засыпали песком. Сейчас радиационный фон там превышен со- 
всем немного, но, если немного копнуть землю показания приборов 
резко скачут вверх. Продукты питания способны накапливать в себе ра- 
диоизотопы, грибы (маслята, свинушки и пр.) способны аккумулиро- 
вать в себе Цезий-137. 

Вышеупомянутые факторы обуславливают актуальность разра- 
ботки устройства детектирования гамма-излучения. 

Для обследования территории с возможным загрязнением обычно 
используют дозиметры, построенные на счетчиках Гейгера. Эти при- 
боры позволяют оценить окружающий фон количественно, но для того, 
чтобы определить с каким конкретно элементом мы имеет дело нужен 
более сложный прибор. 

Гамма-излучение (гамма-лучи, у-лучи) — вид электромагнитного 
излучения, характеризующийся чрезвычайно малой длиной волны – ме- 
нее 2-10 1 м — и, вследствие этого, ярко выраженными корпускуляр- 
ными и слабо выраженными волновыми свойствами [1]. На данный мо- 
мент существуют несколько типов детекторов сцинтилляционных, газо- 
вых, полупроводниковых позволяющих регистрировать гамма-кванты. 

С учетом того, что данное устройство детектирования гамма-излу- 
чения ориентированно на широкою аудиторию целесообразно скон- 
струировать его по более дешевому типу. 

Использование фотоэлектронного умножителя (ФЭУ) совместно 
со сцинтилляционным кристаллом хорошо зарекомендовала себя в об- 
ласти регистрации радиационного загрязнения окружающей среды, ана- 
лизе радиоактивных материалов, спектрометрического исследования 
веществ. Счетчики производят как качественную регистрацию излуче- 
ния, так могут и измерять величину его энергии. 

Конструктивно, устройство детектирования гамма-излучения 
представляет собой блок РЭА цилиндрической формы, состоящий из 
трех функциональных модулей, объединенных внутри герметичного 
корпуса. Структурная схема устройства детектирования гамма-излуче- 
ния представлена на рис. 1. 
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Рис. 1. Структурная схема детектора гамма-излучения 


В устройствах часто используется модульная структура, так как 
она всегда считалось более оправданньгм решением [2]. 

Устройство детектирования гамма-излучения включает в себя: 

+ сцинтилляционный кристалл Nal (ТІ) диаметр 40 мм высота 40 мм; 

+ фотоэлектронный умножитель Hamamatsu R6231; 

* плата узла высоковольтного формирующую питание ФЭУ; 

Ф плата анализатора многоканального, на данной плате располо- 
жены основные элементы управления и обработки, такие как: 
микроконтроллер, ПЛИС, входная аналоговая часть и устрой- 
ство связи с портом USB компьютера. 

Как уже было сказано выше устройство детектирования гамма-из- 
лучения предназначено для регистрации ионизирующего излучения. 
Помимо регистрации активности устройство позволяет при помощи 
специализированного программного обеспечения позволяет определить 
тип найденного изотопа или нескольких изотопов. 

В устройстве детектирования гамма-излучения используется пре- 
образование заряда в цифровой код. При большом заряде с детектора и 
не слишком большой длине сигнального кабеля можно отказаться от 
применения ЗЧУ или использовать буферный линейный усилитель. Оп- 
тимизировать отношение сигнал/шум можно посредством изменения 
времени сбора заряда. Данный способ позволяет обеспечить хорошее 
энергетическое разрешение устройства детектирования. 
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Сигнал с детектора поступает на входную цепь устройства детекти- 
рования, которая адаптирует входной сигнал по форме для последующей 
обработки. Упрощенная схема входной цепи показана на рисунке 2. 


Рис. 2. Схема входной цепи 


Постоянная времени фильтра c, определяется по следующей 


формуле: 
та = Е(С+С,), 


где ty — постоянной времени фильтра, К — сопротивление, C - емкость, 


C, — емкость кабеля. 


При этом формируется сигнал, имеющий фронт нарастания, рав- 
ный по длительности токовому сигналу и спад, определяющийся посто- 
янной времени фильтра 1, . 


Далее сигнал подается на линейный усилитель, который обеспечи- 
вает «укладку» интересующего участка спектра в требуемый участок 
шкалы комплекса. Диапазон регулировки усиления 42 дБ - 8 значений 
через 6 дБ грубо с помощью ключа и 6 дБ с шагом 1/512 плавно (осу- 


ществляется программно). 

Усилитель оснащен восстановителем базовой линии на основе 
управляемого активного дифференциатора, и схемой управления (про- 
граммно) порогом обнаружения полезного сигнала. 

Преобразователь заряд-код собирает заряд на управляемом иде- 
альном интеграторе и с помощью устройства линейного разряда преоб- 
разует накопленный заряд в код. Время сбора заряда определяется дли- 
тельностью импульса одновибратора, регулируемой программно. 

Преобразование сигнала производится по последовательному 
принципу заряд-время-код, что обеспечивает получение высокой одно- 
родности ширины канала без использования дополнительных схем раз- 
равнивания. 

Далее преобразованный сигнал поступает на АЦП. В роли АЦП в 
данной схеме выступает микроконтроллер. Программное обеспечение 
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микроконтроллера обеспечивает работу обмена данньгми с персональ- 
ным компьютером, для последующей обработки полученных данных 
программой, установленной на компьютер. 

В приборе предусмотрена возможность подключения через к ПК че- 
рез USB порт на котором установлено ПО, позволяющее работать с полу- 
ченными от АЦП данными обрабатывать их и представлять пользователю 
в удобном для восприятия виде. ПО, используемое для обработки получен- 
ных данных разработано группой предприятий «ГринСтар». 

Источник высокого напряжения обеспечивает подачу высокого 
напряжения на цепи смещения детектора или для питания фотоэлек- 
тронного умножителя (ФЭУ). Величина высокого напряжения на вы- 
ходе источника задается программно с точностью 1/1024 от максималь- 


ного значения. Время нарастания и спада высокого напряжения также 
задается программно. 

Для реализации поставленных задач были подобраны компо- 
ненты, которые позволят легко спроектировать и изготовить печатную 
плату, имеющие малое знергопотребление, малую габариты и высокую 
надежность. При проектировании было учтено соотношение цена/каче- 
ство, чтобы устройство было доступным. 

Так при проектировании был выбран АКМ-контроллер от фирмы 
ST Microelectronics - STM32FIOIRET6. АКМ-контроллер на основе 
ядра ARM Соцех-МЗ характеризуется низкой ценой, широким набором 
периферии, высоким уровнем рабочих характеристик. Устройство, 
спроектированное на базе данного контроллера легко модернизировать, 
так как контроллеры фирмы ST Microelectronics взаимозаменяемы. 

Из недостатков стоит отметить сложность освоения этих микро- 
контроллеров по сравнению с Arduino, малое количество технической 
литературы, особенно на русском языке. 

Для программирования данного контроллера используются раз- 
личные среды программирования. Наиболее часто используется СооСох 
IDE (СоПОЕ). Данная среда разработки бесплатна, удобна в использова- 
нии, имеется много примеров, доступных для загрузки. Из недостатков 
стоит отметить возможность использования этой среды только под 
Windows, что в нашем случае было не критично. 

Микроконтроллер STM32F101 RET6 на базе ядра ARM Cortex- 
МЗ включает в себя: 32-битный микроконтроллер, выполненный по 
Гарвардской архитектуре; имеет несколько отдельных шин и 3-сту- 
пенчатый конвейер и более 10 регистров общего назначения, что поз- 
воляет выполнять операции параллельно, набор инструкций - Thumb- 
2 (совокупность 16 -и 32-битных команд, ориентированный на ком- 
пиляторы C / С + +). 

5ТМ32Е101 от фирмы ST Microelectronics был выбран не слу- 
чайно. Данная фирма одна из самых инновационных компаний в обла- 
сти полупроводников. Технологический портфель ST Microelectronics 
включает в себя самые передовые КМОП-технологии, а также дополни- 
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тельные технологические опции - встроенная память, аналоговые и си- 
ловые решения. ST Microelectronics имеет всемирную сеть дизайн цен- 
тров и заводов: микросхемы разрабатываются во многих странах Ев- 
ропы, в США, Индии и северных странах Африки. 

Благодаря наличию обширной линейки продукции и мощной тех- 
нической поддержки, компания ST Microelectronics является идеальным 
партнером для реализации комплексных проектов. 

Микроконтроллер имеет более 100 страниц официальной докумен- 
тации и переведенную на русский язык книгу по АКМ Cortex-MG [3]. Для 
освоения данного микроконтроллера выпущены различные платы от- 
ладки, подключаемые к персональному компьютеру через порт USB. 

Используя данный микроконтроллер можно отказаться от допол- 
нительных элементов, таких как ОЗУ, ПЗУ, АЦП и пр. Все это уже со- 
здано на базе этого микроконтроллера. 

Е1аѕһ-память в выбранном микроконтроллере составляет 512 
Кбайт, системная память (SRAM) – 48 Кбайт. Шестнадцать каналов 12- 
битного, достаточно быстрого (время преобразования | мкс), аналого- 
цифрового преобразователя способны обеспечить обработку информа- 
ции даже при больших загрузках входного сигнала. Цифро-аналоговый 
преобразователь отвечает за управление высоким напряжением, питаю- 
щим фотоэлектронный умножитель. 

Микроконтроллер содержит в себе температурный сенсор, обраба- 
тывающий информацию от внешних датчиков по распространенному 


каналу / ues Температурный режим работы данного микроконтроллера 
от —40 °С до 85 С, аустройство при максимальной нагрузке потребляет 


мощность не более 1 Вт, что в наших климатических условиях позволяет 
обойтись без охлаждения (даже пассивного) [4]. 

Отладка микропроцессора может осуществляться подключением 
по двум взаимозаменяемым интерфейсам SWD (Serial Wire Data) или 
JTAG (Joint Test Action Group). Связь с персональным компьютером или 
ноутбуком осуществляется с помощью встроенного универсального 
асинхронного приемопередатчика (Universal Asynchronous Receiver- 
Transmitter, UART). 

Основная схемотехническая обвязка (необходимая для работы 
устройства и не влияющая на его характеристики) формируется исходя 
из технической документации на компоненты и технических требова- 
ний к изделию. 

Питание устройства осуществляется от USB разъема, на который 
приходит 5 В, для работы большинства микросхем необходимо 3,3 В, а 
для некоторых и вовсе 1.8 В, поэтому необходимо поставить преобразо- 
вателей напряжения. Для преобразования 5 В в 3.3 В использован пре- 
образователь фирмы Recom, КПД которого составляет порядка 88 %, 
тепловыделение для данного устройства будет минимальным. Для даль- 
нейшего преобразования 3.3 В в 1.8 В был использован линейный пре- 
образователь фирмы Maxim МАХ8881. 
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Для подключения устройства к компьютеру необходимо исполь- 
зовать наиболее распространенный вид соединения. Шотландская ком- 
пания ЕТО] специализируется в области, связанной с шиной USB. Она 
разрабатывает, производит и осуществляет поддержку устройств и со- 
ответствующих программных драйверов для преобразования последо- 
вательной передачи данных no RS-232 или уровней TTL в сигналы 
шины USB, для того чтобы дать возможность современным компьюте- 
рам использовать простой классический К5-232 интерфейс, а разработ- 
чикам не приходилось вникать в сложности USB интерфейса [5]. Пре- 
образователь UART/USB FT-232R подходит для решения поставленной 
задачи. Все необходимые драйвера находятся в открытом доступе. И 
компьютер будет без дополнительных действий будет воспринимать 
подключенное устройство. 

В аналоговой части устройства детектирования гамма-излучения 
использованы стандартные методы обработки сигнала, которые вклю- 
чают в себя входной усилитель, восстановитель базовой линии, ФНЧ, 
ФВЧ и пр. Сигнал собирается на стробируемом интеграторе, собираю- 
щем полезный сигнал, время сбора задается через ПЛИС, и подается на 
вход микроконтроллера. 

Восстановитель базовой линии управляется программно с помо- 
шью ПЛИС. 

Для питания фотоэлектронного умножителя использована схема 
Ройера, которая позволяет «раскачать» напряжение до 1536 В, ток по- 
требления невысок (около | мА). Управление напряжением осуществ- 
ляется микропроцессором, программируемый шаг значений составляет 
— 1024. Далее напряжение распределяется по диодам фотоэлектронного 
умножителя. 

При выборе фотоэлектронного умножителя были рассмотрены не- 
сколько производителей. Российские производители сильно уступают 
по характеристикам зарубежных ФЭУ и проигрывают в соотношении 
цена-качество. Рассматривать российских производителей стоит в слу- 
чае, если того требует техническое задание. В данном случае такого 
ограничения не было. Выбор стоял между английской компанией ЕТ 
Enterprises и японской компанией Hamamatsu. Были выбраны ФЭУ 
фирмы Hamamatsu, данная фирма является лидером в производстве 
ФЭУ. Фотоэлектронные умножители фирмы Hamamatsu оснащены 
экранами, защищающими прибор от электромагнитных помех, что поз- 
воляет регистрировать даже слабые сигналы. 

Выбранный ФЭУ R6231 имеет диаметр 51 мм, эффективная об- 
ласть которого составляет 46 мм. Это обусловлено диаметром катодной 
сетки, которая регистрирует вспышки, происходящие в сцинтилляцион- 
ном кристалле. Это важно при выборе сцинтилляционного кристалла. 
Если кристалл будет больше диаметром, чем эффективная область реги- 
страции, то не получится собрать весь поток световых фотонов. 

При выборе сцинтилляционного кристалла рассматривались не- 
сколько типов кристаллов: Натрий йодистый легированный таллием 
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Nal (T1), йодид цезия легированный таллием CsI(TI) , бромида лантана 


легированный церием LaBr, (Се) (торговая марка BrilLanCe380) - фран- 


цузской компанией Saint-Gobain Crystals, йодида стронция легирован- 
ный европием SrL, (Би) [6]. 


Сравнив четыре типа кристаллов, безусловно выигрывает йодид 
стронция. Он не имеет собственной активности, как у бромида лантана 
и имеет лучшее разрешение среди всех четырех образцов - 2,9 Фо на 
энергии 662 кЭв по цезию-137. Но данный кристалл сложно купить в 
Российской Федерации, так как производит его всего одна компания, 
находящаяся на территории Соединенных Штатов Америки и финанси- 
руемая непосредственно правительством США. 

Второй образец бромид лантана имеет чрезмерно высокую цену 
для поставленной задачи. Йодид цезия уступает натрию йодистому по 
разрешению и имеет более высокую цену, поэтому было принято реше- 
ние использовать кристалл натрия йодистого легированного таллием. 
Данный кристалл обеспечит хорошее разрешение, порядка 6-8 %, что 
позволит легко обрабатывать полученные спектральные диаграммы. 

Выбирая размер кристалла, стоит учесть следующие факторы: 

* Большие кристаллы обеспечиваю быстрой набор статистики и 
малое время для усреднения результатов измерения, но имеют 
большой размер и вес, для их корректной работы требуется 
ФЭУ с большим диаметром катодной сетки. Для таких ФЭУ 
требуется больший ток питания, что уменьшает время автоном- 
ной работы. Цена данных кристаллов высока. 

* У малых и средних кристаллов время для набора статистики и 
усреднения результатов измерения выше, но данные кристаллы 
имеют малый вес, их можно использовать с меньшими ФЭУ, 
что дает длительную автономную работу. Цена таких кристал- 
лов невысока. 

Учитывая перечисленные факторы, необходимо брать небольшой 
кристалл. Как было сказано ранее диаметр кристалла не должен превы- 
шать диаметр эффективной области регистрации фотоэлектронного 
умножителя, отсюда следует что лучше взять готовый кристалл диамет- 
ром 40 мм и высотой 40 мм. 

Учитывая, что устройство должно быть максимально компактным, 
из разъемов имеет смысл распаять разъем для питания и передачи дан- 
ных - один USB type- В и один светодиод. 

Разъем для питания лучше брать самый распространенный, напри- 
мер, USB гуре- В, а электрическая схема устройства разработана с рас- 
четом того, что питающее напряжение будет равным 5 вольт. Преду- 
смотрена светодиодная индикация работы устройства и нештатных си- 
туаций. Его также можно подключить к GPIO. 

Имеет смысл оставить на плате площадки для SWD и JTAG под- 
ключений. Это нужно для прямой аварийной связи с устройством: SWD 
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позволит получить прямой доступ к микропроцессору, а ТТАС позволит 
напрямую отдавать команды процессору. 

Корпус устройства детектирования гамма-излучения должен 
иметь минимально возможный вес. Так как устройство является носи- 
мым целесообразно применение конструкторского решения в виде гер- 
метичного блока. Для выполнения этих условий корпус целесообразно 
изготовить из алюминиевого сплава Д16Т или аналогичного. Разъем для 
связи с ноутбуком или персональным компьютером используем фирмы 
Amphenol, данный разъем обеспечит защищенность ЇР-65. Для защиты 
самого разъема во время хранения устройства производитель предла- 
гает различные заглушки. 

Программная часть. Условно, программную часть можно разбить 
на несколько частей: системное ПО (загрузчик, драйвера) и прикладное 
ПО (с помощью этого ПО пользователь может видеть визуальное изоб- 
ражение спектра). 

Говоря о системном программном обеспечении, имеется в виду 
прошивка микроконтроллера и ПЛИС. Для написания программного 
обеспечения была использована среда CooCox ШЕ (ColIDE). Для озна- 
комления с микроконтроллером была приобретена отладочная плата. 

Обработка полученных спектров производится с помощью специ- 
ализированного ПО «eSBS», разработанного специалистами «Грин 
Стар». Данное ПО дает наглядное изображение получаемых спектров и 
обеспечивает возможность их динамического обновления (см. рис. 3). 


_ Рис. 3. Эмулятор анализатора eSBS 


Таким образом разработано недорогое, простое в изготовлении и 
эксплуатации устройство, позволяющее осуществлять мониторинг 
окружающей среды, что является одной из причин актуальности его раз- 
работки [7]. 

В данной разработке были использованы различные схемотехни- 
ческие и технологические решения для получения наилучших характе- 
ристик таких как: энергетическое разрешение при использовании деше- 
вого сцинтилляционного кристалла на энергии 662 кэВ порядка 7 %, за- 


грузка спектрометра до 10? импульс/с, что позволяет регистрировать 
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сигнал при высоких знергетических загрузках, интегральная нелиней- 
ность 0,05 Фо, небольшие габариты и малый вес. Низкое энергопотреб- 
ление было достигнуто за счет оптимизации схемы. 

Разработанное устройство детектирования гамма-излучения поз- 
воляет быстро и эффективно определять качественный и количествен- 
ный состав вещества, а также может сигнализировать об угрозе жизни 
человека. 
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PORTABLE OF A DEVICE FOR THE DETECTION 
OF GAMMA RADIATION 


This article presents a developed device for detecting gamma radiation used 
for monitoring the radiation background of the environment and teaching stu- 
dents. The process of device design, selection of the main components of the 
device, selection of auxiliary electronic means, selection of the system software 
part and development of the application software part of the device is consid- 
ered. A brief overview of the capabilities of the developed device is given. 
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УДК 004.925.8;744.4 
О.В. Сулина 


ПРОЕКТИРОВАНИЕ ПРЕСС-ФОРМ 
В CAD-CHCTEME KOMIIAC-3D 


Работа посвящена вопросам преподавания дисциплины «Проек- 
тирование пресс-форм» для подготовки специалистов по специаль- 
ности 15.05.01 «Проектирование технологических машин и комплек- 
сов». Обусловлена необходимость внедрения САПР систем для про- 
ектирования учебных пресс-форм. Предложено внедрить в про- 
грамму обучения систему автоматизированного проектирования 
«Пресс-формы 3D. КОМПАС-30». Рассмотрено применение воз- 
можностей программы для создания литьевой формы для детали - 
заглушка. Проанализированы возможности системы для проектиро- 
ваний пресс-форм, базы данных системы, степень автоматизации 
оформления основного комплекта конструкторской документации, 
соответствие оформления документации стандартам ЕСКД. Ин- 
терактивный режим проектирования пакета пресс-формы позво- 
лит свести к минимуму конструкторские ошибки и выбрать опти- 
мальную конструкцию формы обучающимися при проектировании 
первых пресс-форм, а также значительно сократить время разра- 
ботки одной учебной пресс-формы. 


Ключевые слова: пресс-форма, KOMIIAC-3D, ЗП-моделирование, 
полимер, термопластавтомат, учебный процесс, конструкторская 
документация. 


Термопластичные полимерные материалы на сегодняшний день 
используются практически во всех отраслях промышленности. Со- 
гласно данным [1] в 2018 году мировое производство пластмасс соста- 
вило 359 млн. тонн, 62 млн. тонн в Европе и 4 млн. тонн в России. Одним 
из основных способов изготовления изделий из пластмасс является ли- 
тье под давлением на инжекционно-литьевых машинах (термопластав- 
томатах/ТПА), этим способом изготавливается около 50 % всей пласти- 
ковой продукции. В России ежегодный прирост российского рынка 
ТПА превышает тысячу единиц [2]. 

Ключевым звеном, требующим значительных финансово-времен- 
ных затрат, в конструкторско-технологической подготовке производствен- 
ного процесса обработки полимеров литьем на ТПА является проектиро- 
вание и изготовление технологической оснастки — пресс-форм. При проек- 
тировании пресс-формы конструктор на основе исходных данных (мате- 
риал изделия, сложность поверхностей изделия, технологическая осна- 
щенность предприятия и др.) выбирает рабочую зону, рассчитывает схему 
расположения гнезд, их количество, проектирует литниковую систему, си- 
стему термостатирования и систему выталкивания полуфабрикатов из 
гнезд, причем конструкторское решение на прямую влияет на уровень эф- 
фективности производства и качества полимерных изделий. 
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Интенсивность работы проектировщиков сложной технологиче- 
ской оснастки зависит от слаженной работы специалистов, имеющих 
ценный опыт разработки и внедрения пресс-форм на производстве, и 
молодых инженеров, владеющих современными методами автоматиза- 
ции проектирования пресс-форм в САР-системах. Согласно исследова- 
ниям [3, 4] применение САПР при проектировании пресс-форм позво- 
ляет повысить скорость разработки техоснастки более чем в три раза и 
автоматизировать рутинную работу по деталированию чертежа общего 
вида, созданию сборочного чертежа, спецификации и пр. Поэтому, вы- 
пускники технических ВУЗов, способные участвовать в разработке 
пресс-форм (умеющие проводить типовые расчеты составляющих эле- 
ментов, расчеты на прочность, жесткость, устойчивость, долговечность 
и T. п.) с использованием современных информационных технологий 
широко востребованы современным производством. Такие профессио- 
нальные задачи готовятся решать обучающиеся специальности 15.05.01 
«Проектирование технологических маши и комплексов» (специализа- 
ции «Проектирование механообрабатывающих и инструментальных 
комплексов в машиностроении») при освоении дисциплин «Проектиро- 
вание пресс-форм» и «Проектирование и производство пресс-форм и 
штампов». 

Дисциплина «Проектирование пресс-форм» относится к дисци- 
плинам специализации базовой части учебного плана подготовки спе- 
циалистов и необходима для формирования у обучающихся профессио- 
нальных и профессионально-специализированных компетенций [5]. 
Дисциплина структурирована по следующим видам учебных занятий: 
лекции, практические занятия и самостоятельная работа. Обучение дис- 
циплине «Проектирование пресс-форм» должно строиться на интегри- 
ровании классических теоретических знаний по дисциплине (строение 
и свойства пластмасс, классификация полимеров, конструктивные раз- 
новидности пресс-форм и др.) для решения практических задач по про- 
ектированию различных типов конструкций пресс-форм в самостоя- 
тельных индивидуальных работах. 

Для автоматизации учебных проектно-конструкторских работ и 
актуализации образовательного процесса автор статьи предлагает внед- 
рить в программу обучения систему автоматизированного проектирова- 
ния «Пресс-формы 3D. KOMIIAC-3D». 

Система «Пресс-формы 3D. KOMIIAC-3D» обеспечивает: автома- 
тизированный анализ 30-модели детали и проектирование формообра- 
зующих элементов пресс-формы, моделирование литниковой системы, 
проектирование пакета пресс-формы в автоматическом или интерактив- 
ном режиме, контроль конструктивной допустимости деталей пресс- 
формы, автоматическое формирование в соответствии с ЕСКД основ- 
ного комплекта конструкторской документации, необходимой для изго- 
товления пресс-формы. Система содержит: базу данных оборудования, 
конструкций пресс-форм, параметрические библиотеки чертежей кон- 
структивных элементов пресс-форм [6]. Внедрение системы в учебный 
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процесс позволит повысить практическую составляющую процесса обу- 
чения за счет сокращения времени разработки одной учебной пресс- 
формы и, соответственно, создания благоприятных условий для изуче- 
ния обучающимися на практических занятиях в рамках объема дисци- 
плины специфики проектирования различных типов конструкций 
пресс-форм, систем выталкивания и охлаждения, а также конструкций 
крепления подвижных и неподвижных частей формы. 

Рассмотрим пример применения возможностей программы 
«Пресс-формы 3D. KOMIIAC-3D» для создания литьевой формы для де- 
тали — заглушка (рис. 1): 

1. По чертежу выполняется трехмерная модель детали (рис. 1). 


0001005001 aut" 


| а б 
Рис. 1. Чертеж (а) и трехмерная модель (6) детали заглушка 


2. Анализ исходной детали: автоматически создается технологиче- 
ская модель с локальной системой координат (ось направлена в сторону 
подвижной части и совпадает с нормалью плоскости раскрывания), про- 
водится анализ отрицательных или близким к нулевым уклонам, прово- 
диться анализ раскрываемости, система в автоматическом режиме рас- 
пределяет по адаптивному алгоритму поверхности на поверхности по- 
движной и неподвижной частей формы. Распределение поверхностей 
допустимо проводить и в ручном режиме, причем, при выделении по- 
верхности происходит селектирование связанных поверхностей и авто- 
матический перерасчет линии раскрывания. При распределении должно 
соблюдаться требование к положению изделия в форме: деталь после 
отливки должна оставаться в подвижной части пресс-формы. В соответ- 
ствии с проведённым разбиением поверхностей создаются полумат- 
рицы (технологические файлы создаются автоматически) (рис. 2). 
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а б 
Рис. 2. Полуматрицы заглушка-97.т34 (а), заглушка-98.т3а (6) 


Задается присоединительная точка впрыска, определяющая место 
подвода литниковой системы и направление впрыска материала (указы- 
вается направление, обратное направлению течению полимера в по- 
лость детали). Создается эскиз «Толкатели», где указаны форма и рас- 
положение толкателей (рис. 3). 


Рис. 3. Задание присоединительной точки впрыска, формы 
и расположения толкателей 


3. Моделирование литниковой системы. Выбор гнездности связан с 
учетом объема впрыска, усилия смыкания, пластикационной производи- 
тельности и размеров плит ТПА. Выбираем одноместную пресс-форму. 

Литниковая система должна обеспечить поступление расплавлен- 
ного полимера в формообразующую полость формы с минимальными 
потерями температуры и давления. Моделирование происходит путем 
настройки параметров: задаются параметры радиального, отводного и 
впускного каналов. Для детали заглушка необходим только централь- 
ный литник. Диаметр впускного канала принимаем 4,5 мм (на 0,5 мм 
больше диаметра сопла). Угол конуса 2°. 

4. Проектирование пакета пресс-формы. Задаются параметры 
пресс-формы: 

Ф особенности съема (толкателями - для детали заглушка); 

* центрирующие элементы (колонка-крепление); 

Ф плоскости раскрывания (одна); 

* ТПА (ДА 3130-125); 
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* параметры системы охлаждения — горизонтальное в подвижной 
и неподвижной матрице; 

Ф расположение детали — толщина (в подвижной матрице - 46 мм, 
в неподвижной матрице - 4 мм, длина рабочего хода - 48 мм). 

Определяется конструкция пакета: конструкция фланца в подвиж- 
ном и неподвижном основаниях (с посадкой, с креплением винтами), зада- 
ется смещение, определяющее величину межцентрового расстояния для 
одной пары колонок - крепления (0), определяется критерий однозначно- 
сти сборки (смещением), конструкция колонки-крепления (без лысок), па- 
раметры выборки в основании (присутствует — для уменьшения длины 
центрального литника в неподвижном основании), крепление механизма 
толкания (винтами), наличие защитных лапок (отсутствуют). 

Определяются параметры механизма толкания: размещение 
брусьев (горизонтальное), тип колонки (одним хвостовиком в подвиж- 
ное основание), наличие упора контртолкателя (отсутствует), крепление 
плит толкания (винтами), хвостовик (типовой). 

Определяются параметры впрыска: расположение бурта на литни- 
ковой втулке (посередине — для лучшей фиксации детали в неподвиж- 
ных матрице и основании), контакт с соплом (типовой). 

Задаются размеры пакета (таблица 1): 


Таблица 1 
Размеры пакета пресс-формы для детали заглушка 
основание неподвижное | 30 мм | матрица подвижная 46 мм 
плита неподвижная 0 плита подвижная 30 мм 
плита впрыска 0 плита прижимная 0 
матрица неподвижная 30 мм | основание подвижное | 30 мм 
плита съема 0 


Толщина элементов пакета выбирается согласно требованиям проч- 
ности, материалоемкости, надежности и долговечности пресс-формы. 

В результате проектирования создается математическая модель с 
числовым описанием основных параметров и эскиз пресс-формы с отоб- 
ражением характерных точек и элементов формы, доступные для осу- 
ществления редактирования. 

5. Расчет исполнительных размеров. При расчете исполнительных 
размеров поверхностей деталей формы (матриц, пуансонов и др.), 
оформляющих деталь, необходимо учитывать допуски на размеры де- 
тали и коэффициент усадки материала. Коэффициент усадки материала 
Полиэтилен 102 по ГОСТ 16337-77 1.5 96 — 2 % (принимаем 1,7 %), ко- 
эффициент учета отклонений (0,5). 

6. Создание основного комплекта конструкторской документации. 
В результате проектирования пакета пресс-формы система в автомати- 
ческом режиме создает ЗО-модель пресс-формы (заглушка-0.а3а), сбо- 
рочный чертеж пресс-формы (заглушка-0.сдуу) (рис. 5) и спецификацию 
(заглушка-0.зру), 30-модели (заглушка-1.т3 и др.) и чертежи деталей 
(заглушка-1.сА\ и др.), входящих в пакет пресс-формы. 
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Рис. 5. 3О-модель (а) и сборочный чертеж (6) пресс-формы 


Внедрение системы автоматизированного проектирования 
«Пресс-формы 3D. KOMIIAC-3D» в учебный процесс позволит значи- 
тельно увеличить скорость проектирования пресс-форм за счет исклю- 
чения рутинных операций и расчетов элементов формы, а также автома- 
тизированного создания основной конструкторской документации на 
изделие. Визуальное моделирование пакета позволит сосредоточиться 
обучающимся на принятии важных конструкторских решений и повы- 
сит качество учебных проектных работ. 
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MOLD DESIGN IN CAD-SYSTEM KOMPAS-3D 


The work is devoted to teaching the discipline «design of molds» for training 
experts in the specialty 15.05.01 «Design of technological machines and com- 
plexes». It is necessary to introduce CAD systems for designing educational 
molds. It is proposed to introduce a computer-aided design system «3D Molds. 
COMPAS-3D» into the educational program. The application of the program's 
capabilities for creating a mold for a plug part is considered. The possibilities 
of the system for designing molds, the system's database, the degree of automa- 
tion of the main set of design documentation, and the compliance of the docu- 
mentation with the unified system of design documentation standards are ana- 
lyzed. The interactive design mode of the mold package will minimize design 
errors and allow students to choose the optimal design of the mold when de- 
signing the first molds, as well as significantly reduce the development time of 
one engineering mold. 


Keywords: mold, KOMPAS-3D, 3D-modeling, polymer, thermoplastic auto- 
matic machine, educational process, design documentation. 
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